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Preface 


The following is a brief history of how this book came into ex¬ 
istence. In 1993-94, I developed an interest in the issues of elec¬ 
tronics packaging. By 1995, I could easily simulate an electronic 
system using state of the art computer programs and calculate its 
thermal and vibration characteristics. It became apparent to me 
however that without these sophisticated tools, I had no simple 
way to estimate the same characteristics and hence could not do 
back-of-the-envelope calculations. I noticed that there were plenty 
of good books and references on electronics packaging on the mar¬ 
ket but the majority seemed to make the assumption that the 
reader was already familiar with basic approaches and how to 
make rudimentary calculations. 

Later on, I discovered (much to my surprise) that there are not 
many engineers who have this set of tools. It was at that time that I 
embarked on developing a basic understanding of the engineering 
involved in electronics packaging and subsequently presenting it in 
this book. 

Herein, I have not tried to bring together the latest and most 
accurate techniques or to cover all aspects of electronics packag¬ 
ing. My goal has been to develop a book that can be read either in 
a week's time or over a few weekends and it would provide the ba¬ 
sics that an engineer, mechanical, biomedical or electrical, needs 
to keep in mind when designing a new system or troubleshooting a 
current one. Furthermore, this book will serve as a refresher 
course on an as-needed basis for program and engineering manag¬ 
ers as well as for quality assurance directors. 

This book is based on my seminar notes sponsored by the So¬ 
ciety of Automotive Engineers. 
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Introduction 


ISSUES IN ELECTRONICS PACKAGING DESIGN 

Let us assume that you have the responsibility of developing a 
new electronics system. Let us also assume that your budget al¬ 
lows you to bring a team of experts together. Where do you begin? 
Who do you hire? 

It does make sense to hire a team of electronics engineers to 
design the PCB’s and people to lay out the boards and maybe even 
those who will eventually manufacture them. Also, you have been 
advised that over-heating may be a problem, so you consider hiring 
a thermal engineer but one of your team member’s points out that 
he has a few tricks up his sleeve and it is better to spend the 
money elsewhere. 

In the last leg of your project you hire a junior sheet metal de¬ 
signer to develop your enclosure for you and you send the product 
to the market ahead of schedule. Everyone is happy, but .... 

In a few months, you have a problem. Your field units fail too 
often. The majority seems to have an overheating problem. There is 
a fan to cool the system, but it is not enough, you decide to add 
another one but to no avail. 

Well, your patience runs out and you decide to hire the ther¬ 
mal engineer after all. His initial reaction is that thermal consid¬ 
erations have not been built into the PCB design but after a few 
weeks he manages to find a solution, however, it is expensive and 
cumbersome. Well you have no other choice, you accept his rec¬ 
ommendations and all of the systems Eire retrofitted. 

Before you have a chance to take a sigh of relief, you have an¬ 
other problem facing you. The field units fail again but they seem 
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to have equally different reasons. Some fail at the PCB level while 
others on the surface of the enclosure and others for no apparent 
reason. 

What have you over-looked? 

What knowledge base do you need to develop to answer this 
question? This guide is developed precisely to answer this ques¬ 
tion. Our objectives are: 

• To develop a fundamental grasp of engineering issues in¬ 
volved in electronics packaging. 

• To develop the ability to define guidelines for system’s de¬ 
sign - when the design criteria and components are not fully 
known. 

• To identify reliability issues and concerns. 

• To develop the ability to conduct more complete analyses for 
the final design. 

TECHNICAL MANAGEMENT ISSUES 

Let us review the technical issues that require engineering 
management. These issues are briefly discussed next. 

Electronics Design 

An electronic engineer is generally concerned with designing 
the PCB to accomplish a particular task or choosing a commercial- 
off-the-shelf (COTS) board accomplishing the same tasks. In cer¬ 
tain applications, an Integrated Circuit (IC) or a hybrid must be 
designed specifically for particular tasks. Detailed discussion on 
this topic is beyond the scope of this book. 

Packaging/Enclosure Design 

There are four topics that I categorize under packaging and 
enclosure design and analysis. These are electromagnetic, thermal, 
mechanical, and thermomechanical analyses. We will not cover 
electromagnetics here, however, its importance can not be over¬ 
stated. Unfortunately, much of electromagnetic interference (EMI) 
or electromagnetic compliance (EMC) is done as an after event. 
Testing is done once the system is developed and often coupling 
and interactions are ignored. EMI is difficult to calculate exactly, 
however, back of the envelope estimates may be developed to en¬ 
sure higher end product compliance. Basically, the thermal analy- 
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sis is concerned with calculating the component critical tempera¬ 
tures. Mechanical analysis is concerned with the housing of the 
electronics (from component housing to PCB to enclosure and fi¬ 
nally to the rack) as well as the ability of this housing to maintain 
its integrity under various loading conditions such as shock and 
vibration. Thermomechanical management is concerned with the 
impact of thermal loads on the mechanical behavior of the system. 

In this work, we will set the foundation for thermal and me¬ 
chanical analyses of electronics packaging/enclosure design. 

Reliability 

While in my view thermal, mechanical, thermomechanical and 
EMI analyses are subsets of Reliability analysis, most engineers 
consider reliability calculations to cover areas such as mean time 
to failure (MTTF) or mean time between failures. This information 
helps us develop a better understanding of maintenance and repair 
scheduling as well as warranty repairs and merchandise returns 
due to failure. 
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Basic Heat Transfer: Conduction, 
Convection, and Radiation 


BASIC EQUATIONS AND CONCEPTS 

As electric current flows through electronic components, it 
generates heat. This heat generation is proportional to both the 
current as well as the resistance of the component. 

Once the heat is generated in a component and if it does not 
escape, its temperature begins to rise and it will continue to rise 
until the component melts and the current is disconnected. To pre¬ 
vent this temperature rise, heat must be removed to a region of 
lesser temperature. There are three mechanisms for removing heat: 
conduction, convection, and radiation. 

Conduction takes place in opaque solids, where, using a sim¬ 
ple analogy, heat is passed on from one molecule of the solid to the 
next. Mathematically, it is usually expressed as: 

KA 

0 = — (T’hot-T’cold) (2- 1 ) 

In this equation, Q is heat flow, T is temperature, K is 
thermal conductivity, A is cross-sectional area and L is the length 
heat travels from the hot section to the cold. 

Convection takes place in liquids and gases. The molecules in 
fluids are not as tightly spaced as solids; thus, heat packets move 
around as the fluid moves. Therefore, heat transfer is much easier 
than conduction. Mathematically, it is expressed as: 
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Q — hA(Tfo 0 1 T'eold ) 


(2.2) 


In this equation, Q is heat flow, T is Temperature as in equa¬ 
tion (2.1) but h is defined as the coefficient of heat transfer. A is 
cross-sectional area between the solid generating heat and the fluid 
carrying it away. 

Radiation takes place as direct transfer of heat from one re¬ 
gion to another. Similar to light, it does not require a medium to 
travel. It is expressed as: 

Q = osA(t£ qX -r c 4 0ld ) (2.3) 

Similarly, Q is heat flow, T is Temperature as in the other two 
equations. A , too, is the area of the region, however, s is emissivity 
- a surface discussed later - and a is a universal constant. 

We will talk about these equations in some detail and will 
learn how these equations will enable us to either evaluate the 
thermal performance of an existing system or set design criteria for 
new systems to be developed. We need to bear in mind that in gen¬ 
eral, these equations express physical concepts but do not produce 
“locally exact” solutions. For now, let me draw your attention to a 
few important points. First, there has to be a temperature differen¬ 
tial for heat to flow; next, heat rate depends on the cross sectional 
area; finally, while the relationship between heat flow and tem¬ 
perature difference is linear for conduction and convection, radia¬ 
tion-temperature relationship is extremely nonlinear. 

GENERAL EQUATIONS 

If we need to obtain a locally exact solution, we need to em¬ 
ploy a more general set of equations. These equations sire based on 
conservation of mass, conservation of momentum, conservation of 
energy and a constitutive relationship. This general form of equa¬ 
tions for fluid flow and heat transfer is as follows: 


P,t +{p“ih=0 

( puj),t + puk u j,k = ~Pj +(^ u k,k)j +(M«ij +u j,i)),i +ff (2.4) 


(pCp T ),t +P C p u k T ,k = -P u k,k + ( kT ,j),j + M»kjk) 2 + M( u i,j +u j,i> 


J+ 


P = P(P,T) 
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In these equations a comma denotes taking a derivative. Indi¬ 
ces t, i, j, and k denote time and spatial directions x, y, and z, re¬ 
spectively. 

Clearly, unlike the previous set of equations, these equations 
are not simple to solve nor can they be readily used to evaluate 
system performance or be used to set design criterion. Generally, it 
takes sophisticated computer hardware and software to solve these 
equations. The significance of these equations may be numerated 
as follows. 

1. They produce exact solutions of any thermal/flow prob¬ 
lems. 

2. These equations could be reduced to the simpler forms in¬ 
troduced earlier. 

3. They are used to develop a set of parameters that enables 
us to evaluate system parameters and design criterion 
above and beyond the information given to us by the previ¬ 
ous set of equations. These parameters are non- 
dimensional and can be used as a means of comparing 
various variables among systems that have different con¬ 
figurations such as size or heat generation rate. 


NONDIMENSIONAL GROUPS 

Most often results of engineering research and works in fluid 
flow and heat transfer are expressed in terms of nondimensional 
numbers. It is important to develop a good understanding of these 
nondimensional numbers. The set of interest to us is as follows: 

Nusselt Number 

The Nusselt number shows the relationship between a fluid’s 
capacity to convect heat versus its capacity to conduct heat. 


K 

Grashof Number 

The Grashof number provides a measure of buoyancy forces 
of a particular fluid. 
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Prandtl Number 

The Prandtl number shows the relationship between the ca¬ 
pacity of the fluid to store heat versus its conductive capacity. 



K 


Reynolds Number 

The Reynolds number gives a nondimensional measure for 
flow velocity. 

M 

We will revisit these equations and their significance in elec¬ 
tronics enclosure thermal evaluation later. 
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Conductive Cooling 


As it was mentioned earlier, conduction 
takes place in opaque solids, where, using a 
simple analogy, heat is passed on from one 
molecule of the solid to the next. 

Let us look at an example: 

Consider a layer of epoxy with a thermal 
conductivity of 0.15, a thickness of 0.01, and 
a cross-sectional area of 1. A heat source on 
the left-hand side generates a heat load of 
100. The surface temperature on the right- 
hand side is 75. What is the surface tem¬ 
perature on the left-hand side? For the sake 
of brevity, ignore the units. 


e = ^(r hot - 

^cold) 

or Q = ^(AT) 

in o 
'I o 
o 

II 

o 

o 

(AT) 

=> AT = 6.67 


T = 75+ 6.67 = 81.67 




-3 


II 


L = .01 

Cross Sectional Area * 1. 


Notice that this formula only gives the temperature at one 
point; namely, the left hand side. However, the temperature distri¬ 
bution in the epoxy is not known. This distribution can only be 
calculated by using other mathematical formulae. 
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THERMAL RESISTANCE 


Similar to electrical resistance to current flow, any given ma¬ 
terial also resists heat flow. This concept is very useful and can be 
developed to provide a systematic approach to solving heat flow 
problems. 

In electricity the relationship between the electric potential 
and resistance is defined as 

R 

ac 1 - id •—/S/S/VN—® 


where I is the electrical current. A similar relationship may also be 
developed for temperature, thermal resistance and heat flow. 


AT = QR if 



The previous example may now be solved using this approach. 
By using the concept of thermal resistance we obtain: 


R = — or R = - 0 01 - = 0.0667 
KA 0.15x1 

AT = QR => AT = (100) x (0.0667) = 6.67 
T = 75 + 6.67 = 81.67 


While a veiy simple problem was 
used to demonstrate the thermal resis¬ 
tance concept, this method can be applied 
to complicated problems with relative 
ease. 

Sample Problem and Calculations 

Consider this geometry of a typical 
chassis wall. Find the hot temperature if 
the wall temperature is maintained at 75 
°F; each opening is 5 x 1 in.; sheet metal 
is 0.050 in. thick aluminum (6061). The 
length is in inches, heat flow in BTU/hr 
and temperature in degrees Fahrenheit. 

Before tackling this problem, we 
need to know about thermal resistance 
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networks so that this and similar problems may be modeled prop¬ 
erly. 

Resistance Network 


Similar to flow of electricity through a network of various 
components, each having a different electric resistance, heat, too, 
may flow through different paths in parallel and/or in series, each 
having different thermal resistance. Thermal networks developed in 
this fashion provide a powerful tool to find an equivalent resistance 
for the entire network, hence allowing us to evaluate a temperature 
difference. 

Network Rules 



Since the elements of 
this network are either in 
series or in parallel, we 
first need to know how to 


Parallel 


find the equivalent resistance for each one. 

Series Rules 

When components are placed in series, the overall thermal re¬ 
sistance of such a network increases. 

^ total = R\ + #2 + ^3 + 

Parallel Rule 

When components are placed in parallel, the overall thermal 
resistance of such a network decreases. 



Sample Problem and Calculations 

Consider the chassis wall again. We need to find the hot tem¬ 
perature if the wall on the right hand side is maintained at a 75 °F 
temperature. The first step is to develop the representative net¬ 
work, then reduce it and finally find the equivalent resistance (Rj ). 
This process is depicted in Figure 3.1. 
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In Table 3.1, the length, area and resistance of each element 
is tabulated. Recall that R = L/KA and thermal conductivity for 
aluminum is 7.5 (Btu/(hr ft °F)) 




Figure 3.1. The Thermal Resistance Network of the Chassis Example 
Table 3.1. The Information Pertinent to The Chassis Example 



Now we need to find the equivalent resistance for the elements 
in series: 


111111 
-+-+- H -+-+-=> 

R 2 R 5 R 6 R 7 

111111 
6^67 09 + 09 + 09 + 8.89 + 444 


R 9 =1.21 


This enables us to replace the network representative with a sim¬ 
pler one in which the resistance elements are in series: 
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Rj — R\ + Rg + /?g => = 0.25 + 1.21 + 0.25 /fy —1.71 

A T = QR t => AT = 20x1.71 => AT = 34.2 °F 

r hot = 75 + 34.2 => r hot = 109.2 °F 

Two points must be noted here: 

1. No temperature variation in the vertical direction has been 
taken into account. 

2. We have only calculated the temperature at the high point. 
No other temperature information is known to us through 
this calculation. If critical components are placed inside, 
how do we know that we have not exceeded their operating 
temperature range? 

Before answering this question, we need to verify that we have 
a good solution here. Since this is a relatively simple problem, we 
can find a solution with a high degree of accuracy using finite ele¬ 
ment methods. 


Comparison with Exact Results 


Figure 3.2 may be considered to be the exact results for this 


problem using finite element analysis. The maximum temperature 
from this analysis is also 109.2°F. However, one will notice that the 
temperature distribution along the left edge is not uniform. The 
Resistance Network method has predicted that the temperature all 
along the hot side is uniform. 
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Figure 3.2. The “Exact” Solution Obtained From Finite Element Methods 
Assumptions 

The reason for this discrepancy is that in the previous tech¬ 
nique, it is assumed that heat flow is uniform along the direction of 
the thermal resistance. Effectively, this means that heat conduc¬ 
tion problem is one-dimensional. Clearly, this assumption does not 
hold true all the time as in the comers of this example problem. 
However, it has validity if used with caution. 

Temperature at Intermediate Points 

Recall the point made earlier in regards to placing critical 
components and the calculation of the internal temperature distri¬ 
bution. One needs to bear in mind that the relationship AT = QR 
not only holds true for the entire network but also for each element 
as well. Therefore interior temperatures may also be calculated. 
The only difference is that instead of the total resistance of the en¬ 
tire network, the proper resistance associated with the location 
must be used. Furthermore, keep in mind that Q is constant 
throughout the system and flows in the same direction. 

For example, temperature on the right side of the chassis 
openings is: 
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A T = QR % => AT = (2x3.41)x0.25 => AT = 1.705 

^"rightside = ^-705 => ^rightside = 76.705 

Similarly, the temperature on the left-hand side of the same 
openings is 

A T = Q(R s +R 9 ) => A7’ = (2x3.41)x(0.25+1.21) => AT = 9.96 

^leftside ~ 75+ 9.96 => ^leftside —84.96 


Exercise: 1C Temperature Determination 

One area of thermal modeling is the heat flow in between vari¬ 
ous layers of materials. An example of this configuration i s heat 
flow form a chip into its casing and its heat sink as shown in figure 
1 3.3. 1 As various surfaces come in contact with each other special 
considerations must be given to the interface and its numerical 
modeling. 

In this military application, all the heat is transferred via con¬ 
duction. As a result, spacers must be used to transfer the heat effi¬ 
ciently. In the selection of spacers, care must be exercised to 
choose compatible materials in their thermal expansion coeffi¬ 
cients. This topic will be discussed in Chapter 8. For now, the fol¬ 
lowing thermal coefficients may be used: 


Adhesive = 0.450 Btu/(hr ft °F) 
Silver = 280 Btu/(hr ft °F) 
Copper = 220 Btu/(hr ft °F) 
Insulation = 0.2 Btu/(hr ft °F) 


There are several issues here requiring us to exercise caution. 
First, the set of units presented is not consistent. Next, the thick¬ 
ness of copper in via holes is given in terms of its weight. Finally, 
the proper conduction area for the vias must be calculated. With 
these in mind, the numb er crunchin g is straightforward. The 


thermal network is shown in Figure 3.4 
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IC Case 


BLdhes,™, L= 008/12, A= 1/144, K=0.450 


Rspacer, L=.05/12, A= 1/144, K=280 


RadJiesive, L=.008/12, A= 1/144, K=0.450 


Rsovias, L=.100/12, A=9.756x10-3/144, K=220 


Rspreader, L=.05/12, A=l/144, K=280 


Rinsulator, L=.005/12, A= 1/144, K=0.200 


Radhesive, L=.003./12, A=l/144, K=0.450 
Metallic Core 


Figure 3.4. The Heat Flow Network 
Now the temperature of the IC may be calculated: 


A T = QR t => AT - (2x3.41)xl.37 AT = 9.34 °F 
T IC = 85 + 9.34 => T IC = 94.34 °F 


Notice that 2 watts of heat generation must be converted into 
Btu/hr. 

A strong feature of this approach is in its ease of modeling; 
making changes to the design and comparing results. For example, 
should we decide to add an adhesive to the bottom of vias and the 
silver spread er, we could easily observe the impact on the IC tem¬ 
perature. See Table 3.3. 


A T^QR t => AT = (2x3.41)x 1.58 AT — 10.77 °F 
T IC - 85 + 9.34 => T IC - 95.77 °F 
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Another feature is that the elements with the greatest resis¬ 
tance may readily be identified and hence optimized, should the 
need be. 


Table 3.3. Element Data for The Modified Heat Flow Problem 



MESZm 

Area 

Conductivity 

Resistance 

Adhesive 

.008/12 

1/144 

0.450 

0.213333 


BTiKW 

1/144 

280 

0.0021 

Adhesive 


1/144 

0.450 

0.213333 

50 Vias 

•1/12_ 


220 

0.5591 


Adhcsi ■ 

■ 008/141 

M M 



■o.2i; 

3338® 

E3HH9HI 


1/144 

280 

0.0021 

Insulator 


1/144 

.2 

0.3 

Adhesive 

.003/12 

1/144 

0.450 

0.08 




Total Resistance 

1.583299 


HEAT SPREADING 

In real life applications, as heat flows from one surface to an¬ 
other it spreads in a spherical fashion. In other words, it progres¬ 
sively covers a larger area. The formed cone angle depends on the 
thermal conductivity of the substrate material. By exploiting this 
physical phenomenon, significant issues in semiconductor pack¬ 
aging may be overcome. 

For example, CVD (Chemical Vapor Deposition) diamond has 
extremely high thermal conductivity (1500 W/mK and a spread 
angle (0) of nearly 80 degrees). This product could ideally be used 
as a heat sink for semiconductors or package heat-dissipation 
substrates. However, there are two drawbacks to this material; one 
is that due to its high price it is generally used for cooling of more 
expensive systems such as laser diodes, optical communications, 
or high-output semiconductor laser diodes. Another drawback is 
that diamond’s small coefficient thermal expansion (CTE) com¬ 
pared to larger CTE’s of other components which potentially cre¬ 
ates a reliability concern. As a solution to this problem, composite 
structure heat sinks are developed in which CVD diamond is 
formed either on a silicon substrate or combined with nickel in 
electrolyses plating. These products address the cost and CTE is¬ 
sues of the CVD diamond alone. 
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Figure 3.4. Heat Spreading Cone 

In the analysis, one has to exercise caution. Since heat 
spreads, a larger area must be considered. A good approximation is 
to use an average area of the top and bottom surfaces. To calculate 
this area, the spread angle must be known. The following two em¬ 
pirical formulae are available: 

/ „ \ 0.6 

0 = 90tanh{0.355fj } For Conductivity in W/(m C) (3.1) 

( rr \0-55 

0 = 90tanh{0.51[ J } For Conductivity in Btu/(hr ft F) (3.2) 

Note that K is substrate’s thermal conductivity. For a one¬ 
dimensional problem depicted in Figure 3.4, heat spreading is as¬ 
sumed in the x-y plate alone and the depth is t : 

x = y tan (0) 

Average Area = (d i + (d + 2x) i }/2 or 

Average Area = (d + x) (. 

Example 

A silicon chip measuring 0.07 in. x 0.07 in., 0.025 in. thick is 
mounted to an Alumina case of 0.025-in thickness with conductive 
epoxy 0.003 inches thick. Determine the chip temperature if the 
heat dissipation in 0.35 watts and ambient temperature is 75 °F. 
Thermal conductivity of Alumina is 17 Btu/(hr ft F) and of thermal 
epoxy is 1.25 Btu/(hr ft F). The First, step is to calculate the top and 
bottom area: 

Top Area = (0.07)(0.07) -.0049 in 2 

Bottom Area = (0.07 + x)(0.07 + x) 
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where x = y tan(0) and 0 is calculated from equation (10). 


0 = 90tanh{0,51 


} = 90tanh{0.51 


(3.14)(17) 

180 


} => 0=23° 


From here: 

x = .025 tan(23) 
x = .0106” 

Bottom Area = 0.0065 in 2 

Avereage Area = (.0065 + .0049)/2 Average Area = 0.0057 in 2 



Note that lengths must be converetd to feet and heat generation to 
BTU/hr. 


JUNCTION-TO-CASE RESISTANCE 

Many chip manufacturers provide a set of data commonly re¬ 
ferred to as junction-to-case resistance, which is defined as follows. 



Where j denotes the junction of the chip; c the outer casing; and 
qj the heat generated by the chip. 0j. c is the thermal resistance 

between the die and its outer casing. It may be used - with caution 
- to develop a system’s resistance network and hence calculate the 
die temperature without the exact knowledge of the material(s) 
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used in the chip encapsulation. Furthermore, this is a handy for¬ 
mula for evaluating the temperature of various chips in conjunc¬ 
tion with choosing a heat sink or a fan-heat sink combination. 


CONTACT INTERFACE RESISTANCE 

In many applications, components must be tied to together to 
develop the needed configuration. For example, PCB’s must be 
placed on edge guides or power supply units must be bolted to the 
chassis. This requires that two separate surfaces be joined at an 
interface allowing heat to flow across this interface. Because of 
surface nonregularities, the actual contact area is much smaller 
than what seems to be the contact area. Earlier, we discussed that 
the magnitude of heat flow area has a direct impact on the tem¬ 
perature difference; the larger the areas, the lower the temperature. 
Therefore, smaller contact area in the interface leads to higher 
temperature rises than expected. 

This problem is illustrated below. In Figure 3.5a, a close-up of 
the interface is depicted. Where the surfaces meet, the tempera¬ 
tures are lower (Figure 3.5b) and consequently, where there is a 
gap, the temperatures are higher. Similarly, in Figure 3.5c, the 
heat flux through the interface is shown. Clearly, the interface pre¬ 
sents a thermal barrier that needs to be addressed. 



a) b) c) 

Figure 3.5. Heat Flow Conditions at the Junction of Two Surfaces; 
a)Depiction of the Interface, b) Temperature Variation, c) Flux Variation 

This problem must be taken into account in thermal design or 
analysis. One remedy is to apply pressure. Another solution is to 
apply interface materials such as thermal grease. As pressure in¬ 
creases, the interface resistance decreases however, this decrease 
is bound by an asymptotic value. 

Modeling the Interface 

To model this in a network use the following relationship 
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AT = QR int where *int = , a 

^nt A 

In this relationship A is the apparent contact surface area and 
/zj nt is the interface resistance. Note that h mt increases with both 
temperature and pressure leading to an overall decrease in resis¬ 
tance. For instance, for A1 7075, h int varies from about 500 (Btu/hr 

ft 2 °F at zero pressure and 25 W/in 2 ) to nearly 5000 (Btu/hr ft 2 F 
at 400 psi pressure and 150 W/in 2 ). 

It is worth mentioning that while high altitude does have an 
impact on this resistance, it is not significant for most earth-bound 
applications. However, for space applications care must be exer¬ 
cised in devising appropriate interface pressure. 

Joint Resistance vs. Contact Pressure 



Figure 3.6. The Impact of Pressure on Interface Resistance. The pressure 
is in MPa, and the thermal resistance in Cm 2 -C/W. This picture is adapted 
from M. M. Yovanovich, J. R. Culham and P. Teertstra (1). 

Exercise - Calculate the Component Temperature 

Several heat generating components dissipating a total of 12 
watts are placed on the long end of an L-shaped Aluminum bracket 
as shown. This bracket is mounted on a chassis maintained at 100 
°F. Determine the component temperature. The bolt exerts 25 psi 
of pressure. As means of comparison, a finite element analysis re¬ 
veals that the temperature is 160.44 °F. 

There are three different ways to model this problem depend¬ 
ing on the thermal paths to be considered; namely, the heat flow 
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through the bracket and the heat flow through the interface. These 
will be examined here. All approaches, however, assume a 1-D 
conduction only heat flow. 



Bracket 

Segment 1: From heat source to the chassis; heat path length 
is measured from the center of heat source to the chassis; heat 
path area is the cross sectional area of the bracket. 

Length = L = 8 inches 
Area = 3 x 0.25 = 0.75 in 2 
K = 90 Btu/(hr ft °F) 

■^segmentl 


Segment 2: The portion of the bracket in contact with the 
chassis; heat path length is the bracket’s thickness; heat path area 
is assumed to cover only the area between the bolt and segment 1. 


T2 


90 x 


0.75/ 


1.42 


T44 
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This is justified on the basis that heat will not spread to the entire 
area under the bolt. 


Length = L = 0.25 in 
Area = 2 x 3 = 6 in 2 
K = 90 Btu/(hr ft °F) 


R 


segment2 


0.25/ 

/l 2 

90 x 6/ 
yUX /144 


0.0056 


^bracket ^segment! -^segment2 


1.426 


Interface 

Area = A = 2x3 = 6 in 2 


For A1 7075, hint varies from about 500 Btu/(hr ft 2 F) at zero 
pressure and 25 W/in 2 to nearly 5000 Btu/hr ft 2 F at 400 psi pres¬ 
sure and 150 W/in 2 ). Based on a linear interpolated of this data, 
hint may be found. 


hint = 781.25 Btu/(hr ft 2 °F) at 25 psi pressure, 
^interface = " =0.031 


781.25 x 6 / 


144 


Combination 

^ total = ^bracket ^interface = 1.426 + 0.31 = 1.457 
AT = QR = (12 watts x 3.41 (Btu/hr)/watts)( 1.457) 

AT = 59.6° F=>Component Temperature = 100+ 59.6 = 159.6° F 

A Second Approach 

For the sake of argument let us consider a slightly different 
thermal path and compare the results. In this approach, consider 
the length of the bracket to be 12 inches (=8+4) and the full area of 
the interface. 

Bracket 

Length = L = 8+4 = 12 in 
Area = A = 3 x 0.25 = 0.75 in 2 
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K = 90 Btu/(hr ft °F) 


R Bracket 




2.13 


Interface 

Area = A = 4 x 3 = 12 in 2 

hnt = 781.25 Btu/(hr-ft 2 -F) at 25 psi pressure. 

This data was linearly interpolated from the information given ear¬ 
lier that at 0 psi h = 500 and at 400 psi h=5000. 

^interface = - 1 10/ = 0.01536 

interlace 7gl _ 25x 12/^ 


Combination 

R-total ~ ^bracket ^ ^interface — 2.13 + 0.015 — 2.145 
AT - QR - (12 watts x 3.41 (Btu/hr)/watts)(2.145) 

AT = 87.77° F=>Component Temperature = 100 + 87.77 = 187.77° F 
Clearly, this is a very conservative solution. 

A Third Approach 

A third approach is similar to the first approach, however, to 
use the spread angle to estimate the effective length of the bend in 
the bracket. For a material with the thermal conductivity of 90 
Btu/(hr ft °F), the spread angle is approximately 53 degrees based 
on equation (3.2). 

Bracket 

Segment 1 

Length = L = 8 in 

Area = A= 3x 0.25 = 0.75 in 2 

K = 90 Btu/ (hr ft °F) 
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D _ 

^segment 1 


Segment 2 

Length = L = 0.25 in. 

Area = A = {0.25tan(53) =0.33} x 3 = 0.99 in 2 
K = 90 Btu/ (hr ft °F) 


/ 


T2 


90 x °%4 


■ = 1.42 


R 


segment2 — 




0.034 


^bracket ^segment 1 + ^segment2 1-454 


Interface 

Area = A = |0.25tan(53) =0.33 + 0.25} x 3 = 1.74 in 2 
bint = 781.25 Btu/(hr ft 2 °F) at 25 psi pressure. 


^interface 


1 


781.25 x 1-74, 


144 


0.106 


Combination 

Rtotal ~ ^bracket 4" ^interface ~ 1-454 + 0.106 — 1.56 
A T -QR- (12 watts x 3.41 (Btu/hr)/watts)(1.56) 

AT = 63.8°F=>ComponentTemperature = 100 + 63.8 = 163.8°F 


The finite element solution is 160.44°F and the solution by 
this approach is both conservative and yet close to the exact solu¬ 
tion. The component temperatures are tabulated in Table 3.4. 


Tab 


e 3.4. A Comparison of Component Temperatu res 


Approach 

Component 

Temperature 

Exact 

Solution 

1 

159.6 

160.4 

2 

173.28 

160.4 

3 

163.8 

160.4 
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It is obvious that the first approach while providing the small¬ 
est margin of error, it was not conservative. The second approach 
provided a very conservative solution while the third approach pro¬ 
vided a conservative and yet realistic solution. 

A Word on Edge Guides 

Board edge guides can be used to conduct heat away from the 
PCB. They also provide an added measure of vibration rigidity. The 
higher the contact pressure between the PCB and the guide, the 
lower the thermal resistance. Typical values range from 20°F 
in/watt for the U or G guides to 10 °F in/watt for the Z type guides 
to 4 °F in/watt to wedge clamps. 


2-D OR 3-D HEAT CONDUCTION 


It is possible to form a network of thermal resistance elements 
to cover an area (2-D) or a volume (3-D) and then write the equa¬ 
tions to balance the heat in and out of the system. This is, in a 
way, the basis of Finite Elements Analysis (FEA); however, from a 
practical point of view, for truly 2-D or 3-D problems, a commercial 
FEA package must be utilized to minimize both time and possibility 
of errors. 

In this approach, the resistance network is first developed and 
simplified. Then, the equivalent conductance (defined as the in¬ 
verse of resistance) network is drawn. Finally the heat flow is bal¬ 
anced at each node. To maintain consistency, it is assumed that at 
each node, all surrounding nodes have higher temperatures. The 
following problem provides an illustration. 

' The 


The enclosure contains 3 PCB’s as shown in figure 3.8. 


smaller board is 8.5 x 6.5 inches and the larger boards are 11.5 x 
6.5 inches. The gap between the inner walls and is 0.25” on all 
sides, except for the smaller board where the gap on one side is 
3.25 inches. Each board is 0.05 in. thick and the enclosure wall 
thickness is 0.10 inches. The top surface of the enclosure (on the 
smaller PCB side) is mounted to a surface that is virtually a ther¬ 
mal isolator. 

The voltage and current supplied to this system is 16 volts at 
20 amperes and it operates at 80% efficiency. The top PCB is 
smaller than the other two and dissipates % of heat generated by 
the other two boards. Furthermore, the top surface of the enclo¬ 
sure is insulated and the surface temperature is 75°F. 
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In developing this net work, we are in essence reducing a 
three dimensional problem into a one dimensional problem. One 
underlying assumption is that not only are the boards uniformely 
packed but also each component dissipates a similar level of heat. 
In many realistic problems, this is not the case but this approach 
could very easily put the designer in the right ball park. In this 
situation, to calculate Air 1 R and Air 1 L resistances, the volume 
above PCB 1 in divide d into two regions; the right region and the 
left region (figure 3.9). The heat flow distance (L) of each region is 
taken to be 6 in. while the cross-sectional area is assumed to be 
the one half of the area of the strip covering the gap between PCB 1 
and the surface above it ( = 2(12 + 7)(1.5)/2). Same logic is also used 
for the air gap between the boards. 



Figure 3.8. Three Dimensional View of The Enclosure 
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Figure 3.9. Initial Resistance Network 

This resistance network may be simplified as follows (figure 
3.10). 



Figure 3.10. Simplified Resistance Network 


This resistance network is then converted to a conductance 
network. Recall that conductance is defined as the inverse of 
resistance. 



Figure 3.11. The Conductance Network 
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outsite 



Now, we need to balance the heat equation at each node. To 
ensure that a consistent heat flow direction is followed, assume 
that the temperature of each node is lower than its surrounding 
nodes. 


Q\ +CR X (T 2 -T l ) + CR 3 (T a - 7\) = o 

02 + CR X (T X - T 2 ) + CR a (T 4 -T 2 ) + CR 2 (T 3 -T 2 ) = 0 

£?3 + cr 2 (t 2 - r 3 ) + cr 5 (t a -t 3 ) = o 

CR c (T s -T 4 ) + CR 3 {T x - T 4 ) + CR a (T 2 - T 4 ) + CR S (T 3 -T 4 ) = 0 


This leads to the following matrix equation: 


(CR, + CK 3 ) 
-CR, 

0 


-CR X 

( CR , + CRi ^ 

+ J 

-cr 2 


— 0? 3 — 


0 - C/?3 


-CR 2 

(CR 2 +CR 5 ) 

-cr 5 


-cr 4 

-cr 5 

' CR c + C«3 " 
v + + cr 5 j 


n 


Q\ 

Tl 

► — < 

22 

t 3 


2 3 

Ta 


CR C T S 


The solution to this equation will determine the temperature 
at each PCB as well as inside wall temperature (T 4 ). To develop a 
better feel for this approach let us asign the dimmensions and 
properties used in |Table 3.5| and obtain the following set of linear 
equations. 


' 0.9399 

-0.8495 

0.0 

-0.0904' 

T x 


43.64 

-0.8495 

2.2517 

-1.3239 

-0.0783 

t 2 


87.29 

0.0 

-1.3239 

4.0203 

-2.6963 

h 

* — i 

87.29 

- 0.0904 

- 0.0783 

-2.6963 

5.8740 

Ta. 


225.67 


The solution to this equation is tabulated in 


Table 3.6. 
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Table 3.5. The Dimensions and Properties Used 


Station 

Length 

Area 

K 

Resistance 

Conductance 

Standoff 1 

0.75 

0.05 

34 

5.294 

1.89E-01 

Standoff 2 

0.50 

0.05 

34 

3.529 

2.83E-01 

Standoff 3 

0.25 

0.05 

34 

1.765 

5.67E-01 







Air 1 R 

6.00 

28.5 

0.0153 

165.119 

6.06E-03 

Air 1 L 

6.00 

28.5 

0.0153 

165.119 

6.06E-03 

Air 1 




82.559 

1.21E-02 

Air 2-1 

0.25 

14.2 

0.0153 

13.760 

7.27E-02 

Air 2-4 

3.25 

14.2 

0.0153 

178.879 

5.59E-03 

Air 2-2 

0.75 

55.2 

0.0153 

10.647 

9.39E-02 

Air 2 




12.777 

7.83E-02 

Air 3-1 

0.25 

9.5 

0.0153 

20.640 

4.85E-02 

Air 3-4 

0.25 

9.5 

0.0153 

20.640 

4.85E-02 

Air 3-2 

0.50 

74.7 

0.0153 

5.246 

1.91E-01 

Air 3 




10.320 

9.69E-02 

Air 4-1 

0.25 

4.75 

0.0153 

41.280 

2.42E-02 

Air 4-4 

0.25 

4.75 

0.0153 

41.280 

2.42E-02 

Air 4-2 

0.25 

74.7 

0.0153 

2.623 

3.81E-01 

Air 4 




20.640 

4.85E-02 

Case 

0.10 

236 

0.0153 

0.332 

3.01E+00 

R1 




1.177 

8.49E-01 

R2 




0.755 

1.32E+00 

R3 




0.378 

2.65E+00 

RR1 




1.11E+01 

9.04E-02 

RR2 




12.777 

7.83E-02 

RR3 




3.71E-01 

2.70E+00 

R-Case 

0.1 

236 

0.0153 

0.332 

3.009 

CR1 





8.49 E-01 

CR2 





1.32E+00 

CR3 





9.04E-02 

CR4 





7.83E-02 

CR5 





2.70E+00 

CRC 

_ 




3.01E+00 


Tab le 3.6. The Solution to the Simultaneous Equati ons 


Location 

Power Frac¬ 
tion 

Q 

Temperature 

PCB 3 

0.4 

87.30 

218.11 

PCB 2 

0.4 

87.30 

164.29 

PCB 1 

0.2 

43.65 

137.83 

Wall 

Temp. 



147.53 


Ave. Temp. 

173.40 
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4 

Radiation Cooling 


Radiation heat transfer takes place when two bodies at a dis¬ 
tance exchange heat directly. It has an electromagnetic nature; 
hence, it does not need a medium. And, as the temperature differ¬ 
ence increases, radiation becomes more significant. In the past the 
impact of radiation was generally ignored when temperature levels 
were considered low enough. 

There are two categories of radiation: Wall to Wall Radiation 
and Participating Media. In Wall to Wall Radiation, exchange takes 
place between two surfaces or between a surface and its external 
environment. The medium does not have any influence on the pro¬ 
cess. In Participating Media Radiation, not only do the walls radiate 
to each other but they also exchange heat from the medium as 
well. Therefore, the medium influences the radiation process. 

In a majority of cases in electronics packaging, one need not 
be concerned with the medium participating in the heat exchange 
process. However, should the environment where the electronics 
system is placed contain excess moisture, gasses such CO 2 or par¬ 
ticulate such as soot, this type of radiation must be considered and 
its impact evaluated. Herein, we axe concerned only with the wall 
to wall radiation type. 

FACTORS INFLUENCING RADIATION 

Recall the general radiation equation 

Q = usA(T^-T^) 

More Specifically radiation heat exchange equation between 
two surfaces is written as 
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Qij = - £jAjF jjTj) 

o" is called the Stefan-Boltzman Constant ( 5.67xl0 8 W/(m 2 K 4 ) or 
0.1713x10‘ 8 Btu/(ft 2 h R 4 ). i and j refer to surface i and j. A, is the 
area of surface i. T is Temperature in absolute scale (R=°F+460 or 
K=°C+273). Fjj is the view factor between surfaces i and j. 

Here we have introduced two concepts: 

1. Emissivity ( s ), which is a function of surface properties, 
and 

2. View factor, which is a function of how two surfaces view 
each other. 

Surface Properties 

To develop a better feel for radiation, let us review some defi¬ 
nitions. 

In radiation, there is often a mention of blackbodies. A black- 
body is a standard with which other surfaces are compared. The 
level of radiation emitted or absorbed by a black body depends only 
on temperature and is independent of wavelength. 

Engineering materials including those used in Electronics 
Packaging do not exhibit blackbody characteristics; hence, the 
emitted or absorbed radiation levels are less than an equivalent 
blackbody. The ratio of the emitted energy of a “real” surface to an 
equivalent blackbody is called emittance and is denoted by z k , (A. is 
wavelength). Similarly, the ratio of absorbed energies is called ab- 
sorptance (ct>). Generally, both emittance and absorptance are 
functions of wavelength and temperature as well as surface prop¬ 
erties. For a blackbody emittance and absorptance is equal to 
unity. This also means that what is absorbed by a blackbody is 
emitted. 

A second, more practical idealization may be made. A gray 
body is defined as a body that emits and absorbs energy independ¬ 
ent of wavelength. Thus for a gray surface emittance (now called 
emissivity) is equal to absorptance (now called absorptivity) and 
they are less than unity (a - e < 1). They may still be functions of 
temperature. 

It is clear from this definition that a gray surface does not ab¬ 
sorb all of the incident energy. The remainder has to be either 
transmitted through the body or reflected away. Thus, by denoting 
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p as reflectance, and x as transmittance, one may easily conclude 
that a+p+x = 1. Most solid materials are opaque, so x=0. Therefore, 
a+p = 1 or e+p = 1. This relationship indicates that the higher the 
reflectivity of a gray surface, the lower its emissivity and absorp- 
tance become. 

This finding holds true for many engineering materials used in 
electronics packaging. The heat generated by these enclosures is 
generally in a wavelength spectrum where gray body approxima¬ 
tions are valid. It is interesting to note that the surface color will 
not have an impact in this spectrum and black and white surfaces 
emit and absorb equally. 

Solar radiation falls within a wavelength range where these as¬ 
sumptions do not hold and emissivity may not necessarily be equal 
to absorptivity. In this case the color of surface could have an im¬ 
pact. 

Another factor impacting the value of emissivity is surface 
polish. Surface polish increases p and thus lowers s. Non metallic 
surfaces tend to have larger emissivities than metallic surfaces. A 
polished metallic surface has low emissivity but as the same sur¬ 
face oxidizes, its emissivity increases. 

View Factor Calculations 

Recall the heat exchanged between two surfaces: 

Qij = °( E i A i F ij T ? ~ E j A j F ji T j) 


Fy is called the view factor. It is an indication of what percentage of 


one surface is seen by another surface. Unless one surface is com¬ 
pletely enclosed by another, the view factor between them is less 
than one. For most external radiation calculations, view factor may 
be assumed to be one. This presents a challenge for internal wall- 
to-wall radiation calculation. 

Calculation of radiative energy exchange between any two 
surfaces requires determining the geometrical configuration factor, 
or view factor, between the two surfaces. For two black bodies the 
view factor is defined as the fraction of the diffusely distributed ra¬ 
diant energy leaving one surface i that arrives at a second surface j: 


1 Cos f3 t Cos [3j dA i dAj 

F V = T J i -2- 

A i A, a, nr 


(4.1) 
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where A t and A , are the areas of surfaces i and j, respectively, and 

Pj and Pj are the angles between the position-dependent normal 

vectors to surfaces i and j and a line of length r connecting the 
points of evaluation of the normals (see | Figure 4.I.). 

The derivation of equation (4.1) can be found in Siegel and 
Howell [2]. The basic assumptions used in deriving equation (4.1) 
are that the two surfaces are diffusely emitting and reflecting and 
the two surfaces are isothermal. As a result of these assumptions, 
the view factor depends only on the geometry of the system. It is 
also important to note that for each surface I 

N 

Y^'ij = 1, where N is the number of surfaces. 

In general, equation (4.1) is not used directly; rather, different 
strategies are employed depending whether the model is 3-D, 2-D 
or axi-symmetric. Details of these algorithms and a complete dis¬ 
cussion of the entire view factor calculation can be found in Sha¬ 
piro [3]. 

The major complication with the view factor calculation is the 
possibility of partial blocking or "shadowing" between two surfaces 
by an intervening body. Three types of shadowing may exist be¬ 
tween two surfaces: total self-shadowing, partial self-shadowing 
and third surface shadowing. An in-depth discussion of this topic 
is beyond the scope of this book. In general, an accurate calcula¬ 
tion of view factors is a time and computation intensive task. 

The reader should not be dismayed by these complications. 
The view factors for known shapes and geometries are calculated 
and available in standard handbooks. 

It is often asked whether to include radiation exchange within 
an electronic enclosure. In general, internal radiation may be ig¬ 
nored unless, very specific conditions are encountered. These two 
points are illustrated in the examples below. For a more detailed 
discussion of this subject references (4) and (5) are recommended. 
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Figure 4.1. Surfaces View and Shadow Each Other 

EXAMPLES AND ILLUSTRATIONS 

To illustrate an application of wall to wall radiation exchange 
in Electronics Packaging, the solution to two example problems are 
presented in this section. The first problem consists of the flow of 
air past multiple heat-generating chips surrounded by walls of po¬ 
rous material. This problem involves conduction within the solid, 
porous and fluid materials, convection within the fluid, as well as 
radiation boundary conditions. This example illustrates the impact 
of radiation with the presence of conduction and convection (a 
subject that will be discussed in the next chapter). The second 
problem is flow over a heat generating step in a vertical channel 
and illustrates again the same phenomenon, however, radiation is 
allowed to escape to the environment. 

In these simulations, FIDAP, a Computational Fluid Dynamics 
software, was utilized and four node linear quadrilateral finite ele¬ 
ments have been employed. For more details see (6 and 7). 

Electronics Packaging Problem 

A two-cell vertical channel with heat generating chips was 
modeled; the bottom and top portions as well as the sides of the 
channel were comprised of porous materials. The left-hand channel 
had two heat generating steps and the right-hand channel in¬ 
cluded only one chip. The two channels were connected by a heat 
conducting material and there was convection to the environment 
from the outer walls. The inner walls were assu med to be gra y dif¬ 
fuse surfaces. The complete geometry is shown in | Figure 4.2J Table 
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4.1. summarizes the various material properties. For this simula¬ 
tion the full Navier-Stokes equations including buoyancy effects 
and the energy equation were solved. 

For radiation boundary condition specification purpo ses, 49 


radi ating surfaces were defined; these surfaces are shown in | Figure 


4.3. 


Note that the 27 surfaces in the middle of the computational 
domain were defined as bl ocking surfaces for view factor calcula¬ 
tion purposes. | Figure 4.4 I shows the resultant t emperature distri¬ 
bution for an emissivity of 0.8 while figure 4.5 shows the velocity 
field for the fluid as it passes through the system. To investigate 
the effects of radiation on the resultant temperature field three dif¬ 
ferent values of wall emissivity were employed [e =0.0, 0.1, 0.8). 
The chip em issivity remained constant (s =0.88) for all three cases. 


Figure 4.6 shows the temperature profile over the chips and the 


opposing wall in the left channel while Figure 4.7 depicts the tem¬ 
perature profile over the chip and the opposing wall in the other 
channel. Although the results indicate that for this particular sys¬ 
tem radiation does not have a significant effect on the resultant 
temperature distribution, the radiative interaction between the 
chip and the walls becomes more influential in reducing the tem¬ 
perature as the wall emissivity increases. 


Table 4.1. Material Properties of the Complex Problems 
Properties of air 

Density = 1.929e-5 
Viscosity = 4.689e-05 
Specific Heat = 1005. 

Conductivity = 6.655e-04 

Volume Expansion = .0033 at a reference Temp = 300. 

Conductivity of the solid and porous materials 
Porous base = 0.94 
Chip packages = .038 
Package leads = 3.8 
Central heat sink = 1.143 
Caging = 1.70 

Properties of the radiative surfaces 
emissivity = 0.88 

Stefan-Boltzman constant = 3.6577xl0- n 
Permeability of the different porous materials 
Porosity = 0.48, 

Permeability = 10'^ in the direction normal to the plane of the material 
Heat transfer coefficient at the outer walls 

h = 6.45E-03 at reference temperature = 300. 


§ All units are in watts, sec, inches 
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EFFECTS OF RADIATION IN 



TEMPERATURE 
CONTOUR PLOT 


LEGEND 
A - .3006E+03 
B- .3017E+03 
C- .3028E+Q3 
D- .3039E+03 
E- .3050E+03 
F- .3061E+03 
G- .3072E+03 
H - .3083E+03 
I - .3095E+03 
J- .3106E+03 
K- .3117E+03 
L- 312BE+03 
M - .3139E+03 
N - .3150E+03 
O- .3161 E+03 
P- .3173E+03 
Q- 3104E+O3 
R- .3195E+03 
S- .3206E+03 
T- .3217E+03 

MINIMUM 

.30000E+03 

MAXIMUM 

.32226E+03 


FIDAP 5.04 


Figure 4.4. Electronics Packaging Problem: temperature distribution 



Figure 4.5. Electronics Packaging Problem: velocity field 
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Figure 4.6. Electronics Packaging Problem: temperature along the two 

chips 



Figure 4.7. Electronics Packaging Problem: temperature along a line over 

the one chip 
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Flow in a Vertical Open-ended Channel 


In this problem, the effect of different combinations of modes 
of heat transfer is studied. The geometry is a vertical 2-D channel 
with a heat generating step and convection cooling on the outer 
wall opposing the chip. The conditions are such that flow with a 
uniform temperature enters the channel and passes over the step 
and then leaves the system. Figure 4.8 shows the geometry as well 
as the boundary 7 conditions and the input data. The solid was as¬ 
sumed to be 10 times as conductive as the fluid. A Prandtl number 
of 0.73 and a Reynolds number of 1.3698 was specified for the 
fluid. Four cases were studied: a) conduction only, b) conduction 
and radiation, c) conduction and convection, and finally d) all three 
modes com bined. In this analysis buoyancy effects were ignored. 
Figure 4.9 shows the temperature distribution for each case. It is 
clear that radiation and convection both reduce the overall tem¬ 
perature. In this problem, radiation has a dominant effect because 
heat could escape to the environment via the radiative mechanism. 



Figure 4.8. Geometry and Boundary Conditions 
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Figure 4.9. Open Ended Channel: Temperature Distribution For 
The Different Modes of Heat Transfer 

CABINET SURFACE TEMPERATURE 

An electronics cabinet dissipates 800 watts of energy and its 
physical dimensions are 12 x 18x26 inches. The surface is 
painted white. Assuming an ambient temperature of 50° C and 
considering conduction and radiation alone (no solar energy), we 
need to estimate the interior temperature. 

It is important to develop the ability to develop a “feel” for 
surface and interior temperatures of electronics enclosures with 
minimum information. The only information given about the cabi¬ 
net is that it dissipates 800 watts of energy and its physical dimen¬ 
sions are 12 xl8 x 26 inches. 

To make a back of envelope calculation, it is not necessary to 
know about how the components are mounted inside the enclo¬ 
sure. The fact is that the generated heat will escape; either to all of 
the surfaces or just a few depending on the materials used. 

A good assumption is that the heat is uniformly distributed to 
all of the external surfaces. Another assumption is that surface 
and interior temperatures are equal. For the sake of comparison, 
this case is solved along with the condition that heat is only trans¬ 
ferred to one of the largest surfaces. 
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1) 800 watts is distributed equally to all surfaces. Assume a 
sink temperature of 50°C (122°F) and no solar radiation. 
The generated heat should be set equal to the radiated 
heat (notice that no heat is conducted away from the 
cabinet): 

Generated Heat in Btu/hr = 800 x 3.41 = 2728 Btu/hr 

Radiation Heat Transfer Q = F e A cr (T su , f 4 - Tamb 4 ) where F 
is the view factor, e the emissivity, A, the area, a the 
Stefan-Boltzman (= 1.713 x 10' 9 ) constant and finally T 
temperature. The subscript surf stands for surface and 
amb for ambient. The surface is painted white so emissiv¬ 
ity is about 0.9 and since the cabinet is exchanging heat 
to the environment, the view factor F is 1. 

Ambient Absolute Temperature = 122 + 460 = 582°R 

Total area = 

2 (26 x 12 + 18 x 12 + 26 x 18)/144 = 13.83 ft 2 
2728 =(l)(0.9)(13.83)(1.713xl0' 9 )(T sur f 4 - 582 4 ) 

=> Tsurf = 702°R = 242°F = 116.6°C 

2) 800 watts is distributed only to one of the largest surfaces. 
All other assumptions remain the same. 

Total area = 26 x 18 / 144 = 3.25 ft 2 

2728 =(1)(0.9)(3.25)(1.713 xl0- 9 )(T SU rf 4 - 582 4 ) 

=> Tsurf = 901°R = 441°F = 227°C 

Clearly, every effort must be made to increase the heat 
transfer surface area. 
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5 

Fundamentals of Convective Cooling 


With the exception of space applications, there are hardly any 
practical thermal cooling problems that do not include convection. 
Convection involves heat transfer between the surface of a solid 
and its surrounding fluid. The rate at which heat is transferred is 
evaluated by 


Q — hA(T^ ol T'cold) 


This equation only looks simple! h - the heat transfer coeffi¬ 
cient - varies depending on the flow regime, i.e., laminar or turbu¬ 
lent, geometry and fluid properties. 

If we rewrite the convection equation slightly differently, we 
obtain: 


q = ^(T’wall - T «) 


q = Q/A is called the heat flux and is heat transfer rate per 
unit area. The units for heat transfer coefficient (h ) is (W/(m 2 C) or 
Btu/(hr ft 2 F)). T m is called bulk temperature arid is unaffected by 
the heat input from the solid. In electronics packaging problems, 
this assumption holds at some times and at other times, one must 
exercise caution in using this relationship. 

For a fluid to remove heat from a surface, heat must first be 
conducted into the fluid to be removed. Therefore, convection heat 
transfer not only depends on how fast the fluid flows but also on 
how well it conducts heat near the surface. The ratio between the 
fluid’s ability to conduct heat and then move it away is called Nus- 
selt number. It ties both conduction and convection together. 
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Recall that h c is the heat transfer coefficient, L is a charac¬ 
teristic length, and K y is the fluid conductivity. Effectively, this 

relationship is localized. However, it is also customary to have a 
similar relationship that is averaged over an area; 


where a bar over h or Nu indicates its average value. 

These relationships enable us to compare fluids under various 
flow regimes and conditions. For example, consider the impact of 
turbulence on Nusselt number for long ducts carrying gasses or 
liquids (4): 

For Laminar flow (Re<2100, Pr>0.7): 

Nu = 1.86(RePrD/L) a33 


For Turbulent flow (Re>6000, Pr>0.7): 
Nu = 0.023Re 0 ’ 8 Pr 033 


In these two relationships, Re is Reynolds number and Pr is 
Prandtl number, D is hydraulic diam eter and L is length. This re¬ 
lationship is graphically depicted in |figure 57Q It is quite clear 


from this figure that turbulent flow has a higher Nusselt number 
and hence a greater ability to remove heat. While this figure is ex¬ 
tended for the entire Reynolds number range, it should be kept in 
mind that the region in the 2000<Re<6000 is not governed by ei¬ 
ther of the two relationship and it remains - for this example - un¬ 
clear. 
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— 'Nu - Laminar D Nu - Turbulent 



Reynolds Number 

Figure 5.1 . A Graphical Depiction of Nusselt numbers for Laminar and 
Turbulent Flows. For this example, a length to diameter ratio of 20 and a 
Prandtl number of 0.75 was used. 


Flow Regimes, Types and Influences 


The value of heat transfer coefficient depends on the charac¬ 
teristics of the flow. Categorically, convection may be divided into 
two regimes; Laminar, and Turbulent. It may also be divided into 
two types; Free Convection, and Forced Convection. Compared to 
laminar flows, turbulent flows have a larger capacity to carry more 
heat. Similarly forced convection has a higher capability to remove 
heat. 


The greatest influences on convection are density, and pres¬ 
sure - temperature is related to both of these variables through an 
Equation of State or Constitutive relationship. An appreciable 
change in altitude will create a significant variation in either pres¬ 
sure or density that would warrant a recalculation. Many electronic 
packages are designed for sea level conditions and need to be de¬ 
rated when used in higher altitudes - utilized either on airborne 
system s or used in high altitude places. 

In I Figure 5.2| one may identify three areas of this flow; laminar 
where a coherent flow pattern exists, Transition region where this 
pattern begins to break down and finally turbulent where no co¬ 
herent flow pattern exists. 
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Figure 5.2. Different regions in a fluid flow adapted from Van Dyke, M.(8) 


Generally speaking, in natural convection, laminar flows have 
a Grashof number less than 10 7 , transition regions between 10 8 
and 10 10 . For values greater than 10 10 , the flow is considered to be 
turbulent. For a vertical heated wall, for lengths up to 18 to 20 
inches from the bottom of the wall, the flow remains laminar and 
beyond 24 inches, it is turbulent. 

For flows through ducts, laminar flows have Reynolds number 
less than 2000 and turbulent flows have Reynolds number greater 
than 6000, with the region in between being the transition region. 


FREE (OR NATURAL) CONVECTION 

One technique used in electronics cooling is free (or some 
times called natural) convection. By taking advantage of heat ris¬ 
ing, it is possible to design a system whereby cool air enters the 
enclosure form the bottom and warm air exits from the top. Note 
that the words “cool” and “warm” as opposed to “cold” and “hot” 
are used. When free convection is used to cool the system, tem¬ 
perature rises must be moderate. 

In free convection, fluid flow is developed under the influence 
of buoyancy forces. Depending on the orientation of heat source 
with respect to gravity, the flow field may be non-existent, simple 
or very complicated. It should be noted that flow might start as 
laminar but develop into turbulent. 
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To study the relationship between the buoyancy and the abil¬ 
ity of the fluid to remove heat, we need to reexamine the non- 
dimensional relationships. 


Nu 


hL 

K 


Gr = 


I?p 2 gpAT 

M 2 



Where £ is a characteristic length; p is the density of the fluid; 
g is gravity; /? is coefficient of fluid expansion; AT is the change in 
temperature; p is fluid viscosity; C P is fluid specific heat; h is coeffi¬ 
cient of heat transfer; K is fluid conductivity. 

Experiments (4) have shown that the following relationship 
may be developed: 

Nu = Constant Gr m Pr n 


or 


HL.c 

‘l?p 2 g/3AT 

m 

i 

_i 

K 

» 2 J 


K 


(5.1) 


This equation illustrates the relationship between the heat 
transfer coefficient, the ability of the fluid to absorb heat (C p ), its 

ability to conduct (K ) as well as its inherent buoyancy forces. Con¬ 
stants C, m and n are dictated by the flow regime and geometry. 
Most standard heat transfer textbooks provide this data for a vari¬ 
ety of geometries. In electronics packaging the vertical and hori¬ 
zontal walls (heated side facing up or d own) are quite frequently 


used to model the system’s enclosure, 
stants for these geometries. 


Table 5.1 


provides the con- 


Estimates of Heat Transfer Coefficient 

The equation for the Nusselt number (5.1) may be reduced to 
a specific equation for the heat transfer coefficient. This provides 
simplified means for back of the envelope engineering calculations. 
For exam ple, consider properties of air at 66°F as presented in | ta- 1 
ble 5.2. 
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Table 5.1. Empirical D a ta For Nusselt Number C o nstants 


Wall Orientation 

Laminar Flow, 
Gr<10 8 

Turbulent Flow, 
Gr>10 10 

C 

m 

n 

C 

m 

n 

Vertical 

0.55 

0.25 

0.25 

0.13 

0.33 

0.33 

Horizontal 
heat facing up 

0.14 

0.33 

0.33 

0.54 

0.25 

0.25 

Horizontal 
heat facing 

down 

0.27 

0.25 

0.25 

NA 

NA 

NA 


Table 5.2. Properties of Air at 66°F 


p (lbm/ft 3 ) 

0.076 

C P (BTU/lbm-F) 

0.24 

fj (lbm/ft-sec) 

1.225x10 5 

K( Btu/hr-ft-F) 

0.0147 


1.91xl0- 3 


Now, evaluate Prandtl, Grashof, and Nusselt numbers based on 
this property set. 


Pr = 


RTT J 

(0.24-)(1.225xl0~ 5 

lbmF 


lbm 

ft sec 


)(3600 


sec 

hr 


0.0147 


BTU 
hr ft F 


0.72 (dimensionless) 


I? (0.076^) 2 (32.2 —~y)(l -91x10" 3 -)AT 

Gr = -S-^. .-.. .= 2.3673 x 10 6 L 3 AT 

(1.225x 10 -5 _i^EL> 2 
ft sec 


0.0147 ' 


Equation (5.1) may be rewritten based on these values. 


hL 

0.0147 


C(2.3673x 10 6 1? AT ) m (0.72)" 
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From [Table_5jJ for laminar flows and a vertical wall, we obtain 
m=n=0.25. Thus the following relationship for the heat transfer 
coefficient h may be developed. 


h = 0.575C(—) 0 ’ 25 at 66°F (5.2) 

jL/ 

The value of L is generally determined by the flow path the over 
the heated surfaces. For example, for a flat vertical plate, L is the 
height of the plate; for a flat horizontal plate 

1 11 
L length width 

For turbulent flows, a similar relationship may be developed 
based on m=n=0.33. 


h = 1.863C(A7')°' 33 at 66°F 


(5.3) 


Notice that in this relationship h is independent of length L . This 
is consistent with McAdams findings (9). 

If we conduct a similar calculation for a variety of tempera¬ 
tures between 0°F to 300°F, we will soon discover that the form of 
this equation does not change appreciably. Therefore, the heat 
transfer coefficient for air (0 to 300°F) in free convection maybe cal¬ 
culated with some confidence from the following relationships. 

For laminar flows 


A = 0,56C(^) 025 
For turbulent flows 


h = 1.75C(Ar) 033 


Notice that AT is generally not known a priori and an iterative 
solution must be used. The value of C varies from 0.26 for horizon¬ 
tal plates with the heated side facing downward to 1.45 for small 
components such as resistors and wires. Its value is nearly 0.5 for 
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all other configurations of plates. Steinberg (5) provides a table for 
values of C for a variety of components. 

It is extremely important to remember that these are imperical 
formulae. Therefore the following dimensions must be used: 

h c is in Btu/(hr ft 2 F) 

AT is in degrees F 

L is in feet, a characteristic length 

Solution Procedure 

Our ultimate goal of doing any thermal analysis is to ensure 
that our critical component temperature remains below a safe 
threshold. So far, as we have talked about the heat transfer coeffi¬ 
cient value, we noticed that it depends on the surface temperature, 
which we were trying to find in the first place. It, therefore, be¬ 
comes obvious that we need to employ an iterative solution tech¬ 
nique to find both the heat transfer coefficient and the surface 
temperature. The very first step in accomplishing this task is to 
write the heat equation and balance the outgoing heat versus the 
incoming. Hereon, one of two approaches may be taken. 

Approach 1 

• Assume a temperature difference (AT) 

• Calculate h 

• Calculate the temperature difference based on the generated 

heat, h , and the thermal resistance. 

• Compare the assumed and calculated AT 

• Iterate if needed. 

Approach 2 

• Assume a h 

• Calculate the temperature difference based on the generated 

heat, h , and the thermal resistance. 

• Calculate h 

• Compare the assumed and calculated h 

• Iterate if needed. 

There are three considerations here. First, both approaches 
produce identical results, however, the second one may be more 
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practical to use. The reason is that the range of variations for h is 
by far narrower than that of temperature. In fact, for a moderate 
size enclosure, h is in the neighborhood of one (1). Typically, for 
external flows it has a range of 0.85 to 1.2 (Btu/[hr ft 2 °F]) and for 
internal flows, it varies between 0.75 and 0.85 (Btu/[hr ft 2 °F]). 

Second, we assume that we know the level of the dissipated 
heat. It is generally understood that the power is equal to voltage 
times the current. This power, however, is the dissipated heat 
when and only when there are no outputs from the system. For a 
system with input as well as output, the dissipated heat is the dif¬ 
ference between the power entering and the power leaving. A good 
rule of the thumb is that most systems are 15% to 25% efficient, 
meaning that the dissipated heat is between 15% to 25% of the 
total input power. 

Third, heat flow may take different paths. We need to know 
the level of heat flow through the particular path that is of interest 
to us. For instance, suppose that we have used the aforementioned 
approach to calculate the surface temperature of a sealed elec¬ 
tronic box. How do we calculate the interior air temperature? 
Should we use the total dissipated heat? The answers to these 
questions are based on the heat flow paths. 

The generated heat must take conduction, convection and 
possibly radiation paths to get to the surface. If we are interested in 
the interior air temperature, then, we must consider convection 
path and its heat flow alone. Again, a good rule of thumb is that 
only 25% of total heat pass through this path. This rule of thumb 
may be justified as follows. 

To establish a good natural convection flow, it has been rec¬ 
ommended (see Steinberg (5)) that board spacing be equal to or 
greater than 0.75 inches. Now, consider two parallel boards each 
with an area of A. Assuming that radiation is negligible, the heat 
paths between the two boards are as follows: 

Conduction: Rconduction = L/KA ATl — Qi Rconduction 

Convection: Rconvection = 1 / hA AT2 — Q2 Rconvection 

Total Heat = Heat through conduction + Heat through convection 
or Q=Qi+Q 2 . But ATi =AT 2 = temperature difference between the 
boards. By substituting L=(0.75 in/12) ft, h= 0.75 (Btu/[hr ft 2 °F]), 
and K=0.015 (Btu/(hr ft °F)) (Conductivity of Air), we obtain 
Qi=0.75 Q signifying that heat through conduction is 75% of total 
heat and Q2=0.25Q indicating that heat through convection is 25% 
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of total heat. An example at the end of this segment will clarify this 
point. 

High Altitudes 

As we ascend to higher altitudes, density of air drops. This 
drop could have a detrimental impact on air-cooled electronics 
packages. Recall that the heat transfer coefficient was derived from 
Nusselt number, which is the product of Grashof and Prandtl 
number. Also recall that Grashof number is a function of density. 
Therefore as density drops, the heat transfer coefficient and along 
with it the system’s ability to remove the heat is lowered. So, 


^sea level ^ 

L Ps ea level ^ sea level ^ 

m 

C p p 

— c 
K 

M 2 


K 


and 


^'altitude^ . r~ 

3 2 

^ Paltitude£altitude/^-^ 

m 

~C p p 

— L 

K 

P 2 


K 


By dividing these two equations, we obtain the following 
reationship: 


^altitude 
^'sea level 


Paltitude&altitude 

2 

Psea level £sea level 


Density is a direct function of pressure (Gas Law: P = pRT ) 

and m=0.25 for a large class of electronic systems. Furthermore, 
the change in gravitational constant is not very significant for the 
same class of problem. Thus, it can be shown that 


h 


altitude 


’'sea level. 


P altitude 
Psea level 


= h 


sea level . 


altitude-' sea level 


sea level 1 altitude 


(5.4) 


Pressure and Temperature maybe obtained from atmospheric 

data. 
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Board Spacing and Inlet-Outlet Openings 

Earlier, it was suggested that for a system cooled based on 
natural convection, the board spacing must be at least 0.75 inches. 
Adrian Bejan (10) provides the means of calculating an optimum 
gap spacing for natural flow and the associated maximum heat 
removal: 


D apt =23LRa- y \ 

(5.5) 

HW 1/? 
ft™ “0.45*47-— Ra'l 1 , 

(5.6) 


where: 

Ra L = 


^outlet ^inlet 


av 


H is the transverse length of the entire package and L is the 
length of the boards along the flow, W is the width of the stack, g 

is gravitational constant, a is thermal diffusivity, /? is the coeffi¬ 
cient of volumetric thermal expansion, p is viscosity and v is ki¬ 
nematic viscosity (= pj p where p is density). 

The opening area may be calculated based on the number of 
boards and their optimum gap size. The exit area should be de¬ 
signed as 1.5 times as large as the inlet area. 

Design Tips 

In designing enclosures using free convection one has to take the 
following points into account: 

1) The level of heat to be dissipated. 

a) Free convection is most effective in cooling relatively low 
component densities - thus moderate temperature rises. An 
approximate heat flux limit for free convection is 1 Watt/in 2 
for each board. 

2) The enclosure is sealed. 

a) Consider the PCB orientation with regards to gravity. Prop¬ 
erly oriented PCB’s allow for internal free convection to cir¬ 
culate the medium and lower the overall temperatures. 
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b) Also consider using fins. Again, the orientation of the fins 
(parallel to gravity) is crucial to developing proper flow pat¬ 
terns 

3) The enclosure is not sealed and there are airflow in-takes and 
outlets. 

a) When possible, ensure the vertical dimension is at least 
twice as long as the horizontal dimensions. This allows for a 
tower effect and creates an efficient flow environment. 

b) The in-take openings must be placed on the bottom of the 
enclosure while exhaust vents must be placed near the top 
for this configuration. 

c) In placing the components, locate the ones with the highest 
heat dissipation on the bottom of the enclosure near in take 
vents. 

4) Ensure an efficient conduction path. One way of doing so is by 
placing the highest-power-dissipating boards near the walls of 
the enclosure. This allows for radiation and possibly conduc¬ 
tion to play a role in heat removal as well. Also, place the 
boards at a minimum pitch of 0.75 inches. 

5) The efficiency of this cooling technique drops considerably by 
increasing altitudes. 

Cabinet Interior and Surface Temperature 

Recall the electronics cabi¬ 
net used previously dissipating 
800 watts of energy. Its physical 
dimensions are 12 xl8 x 26 
inches. The surface is painted 
white. Assuming an ambient 
temperature of 50° C and consid¬ 
ering conduction, convection and 
radiation as well as solar energy, 
we need to estimate the interior 
temperature. The system will be 
utilized in Orlando Florida 
(nearly sea level conditions) as 
well as Denver at 7000 feet. Do 
not take wind effects into ac¬ 
count. 

This method makes use of simple free convection and radia¬ 
tion assumptions to evaluate the interior and surface temperatures 
of the electronics cabinet considered previously. Briefly, it balances 
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the input heat energy and the output heat energy to provide an 
overall temperature. It is assumed that no wind is present. The ap¬ 
proach is; first to calculate the surface and next, to evaluate the 
interior temperature based on the surface temperature. 

Surface Temperature 

Input heat = Qelectronics Qsolar 

Output heat = Q CO nduction + Q convection Qradiation 

Where 


a) Qelectronics = 800 watts (given) 

b) Qsoiar = solar absorptivity x solar radiation x Area 

C) Qconduction - (K Aconduction)(Tsurface ” Tambient)/L 

d) Qconvection = h Aconvection (Tsurface ” Tambient) 

e) Qradiation = F £ A CJ (T^surface “ T^ambient) 


Heat Balance 

a) Qelectronics = 800 watts x 3.41 Btu/(hr-watt) = 2728 Btu/hr 

b) Consider that Sun is shinning on the largest surface alone 

A= (26 x 18)/(144) = 3.25 ft 2 

Qsoiar = solar absorptivity x solar radiation x Area 

Qsoiar = 0.2 x 355 x 3.25 = 230.75 Btu/hr - note that white 
paint has a solar absorptivity of about 0.2. Furthermore, only one 
of the largest surfaces (area=3.25 ft 2 ) is assumed to receive solar 
radiation (at noon assumed to be a maximum of 355 Btu/(hr-ft 2 ). 
Clearly, for a more realistic problem, the actual number of exposed 
surfaces must be used. Also, as solar radiation varies depending 
on the time of the day and the physical location of the unit. 

Total input heat = 2958.75 Btu/hr 


C) Qconduction - (K Aconduction) (Tsurface - Tambient) / L 


K = 90 Btu/(hr-ft-F) - Aluminum 6061 T6 

Aconduction = 2 (0.25 x 26)/144 = 0.09 ft 2 - There are two flanges 

on each side. 

L = 4/12 = 0.33 ft 
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Qconduction — (90 X 0.09) (Tsurface " Tambient) /(0.33) 
Qconduction = (24.5)(Tsurface - Tambient) 


d) Qconvection h Aconvection (Tsurfacc - Tambient) 


A = 2(26 x 12 + 26 x 18 + 12 x 18)/ 144 => A = 13.83 ft 2 

Qconvection = h (13.83)(T surface - T ambient) 


c) Qradiation F S A O (T 4 SU rface “ T 4 ambient) 


Qradiation = 1 x0.9x 13.83 X 1.713 X 1 0 -9 (T 4 SU rface - T 4 ambient) 
Qradiation = 21.32 X 10 -9 (T 4 surface “ T 4 ambient) 


It is assumed that the box radiates on all six sides. 

By summing the output heat and equating the algebraic sum to 
the total input heat, a nonlinear equation for solving T sur f a ce- Note 
that Tambient is 582 R. 

Input Heat = Output Heat 


2958.75 — Qconduction 4- Qconvection Qradiation 


2958.75 = (24.5)(Tsurface - 582) + h (13.83)(T surface - 582) + 

21.32 x 10 -9 (T 4 surfac e - 582 4 ) (5.7) 

To solve this equation, first h is estimated and this equation 
is solved for T SU rface- Then this temperature is used to get a better 
estimate for h . A good start for h is 1. 

For h= 1. 

=> T sur face=633.5 °R = 173.5 °F => AT =633.5 - 582 = 51.5 
^calculated = 0.29 (AT/L)°' 2 5 

where L, characteristic length, is =3HWD/(HW+HD+WD) 

L = 3 (26)(18)(12)/(26xl8 + 26x12 + 18x12) = 16.9 inches. 

L must be in feet thus: L = 16.9/12 = 1.4. 

^calculated = 0.29 (AT/L) 0 - 25 = 0.29 (51.5/ 1.4)0-25 

^ ^calculated = 0.7 1 


For h= 0.71 

—^ Tsurface = 637.1 °R= 177.1 °F => AT =637.1-582 = 55.1 
^calculated = 0.29 (AT/L) 0 - 25 = 0.29 (55.1/1.4)0 25 
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—i* ^calculated — 0.73 

For h = 0.73 

==> Tsurface=636.9 °R = 176.9 °F =i> AT =636.9-582 = 54.9 

This solution converges to the surface temperature of 176.9 °F 
and h=0.726. The choice of a characteristic length may be puzzling 
in some cases. In general, the temperature is not super sensitive to 
this variable. In here, had we chosen a value of 0.5 for L, the re¬ 
sults would not have been much different (174.3 °F with h=0.928). 

Calculation of Internal Temperature 

Earlier, the heat distribution between different paths was 
demonstrated and it was shown (in a general sense) that 75% of 
the generated heat escapes to the surface through conduction and 
only 25% of the generated heat is distributed inside of the cabin. 
Following the same heat balance principle: 

1) Input 

Qeiect.conv- = 0.25 x 800 x 3.41 = 682 Btu/hr 

2) Output 

Qoutput — h A (Tintemal - Tsurface) 

=i> 682 = h (13.830) (Tintemal - 176.9) 

Similar to the previous case, to solve this equation, first h is 
estimated and this equation is solved for Tintemal. Then this tem¬ 
perature is used to get a better estimate for h. A good start for h is 
0.75 for internal flows. 

For h = 0.75 

=> T surf ace=242.7 °F => AT =242.7-176.9 = 65.8 
hcalculated = 0.29 (AT/L) 0 -25 

where L, characteristic length is =3HWD/(HW+HD+WD) 

L = 3 (26)( 18)( 12)/(26x18 + 26x12 + 18x12) = 16.9 inches. 

L must be in feet thus: L = 16.9/12 = 1.4 
hcalculated = 0.29 (AT/L)°25 = 0.29 (65.8/ 1 .4)°-25 
^ hcalculated = 0.76 
For h= 0.76 

=^Tmtenai = 241.9 °F => AT =241.9-176.9 = 65 
hcalculated = 0.29 (AT/L)»-25 = 0 .29 (65/1.4)°-2S 
— r* hcalculated = 0.76 
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This solution has converged and the interior air temperature 
is about 242 °F. In order to calculate the component or board tem¬ 
perature, we can use the convection equation again with h = 0.76 
and Tbuik = Tintenai = 242. Suppose that one board of interest dissi¬ 
pates 100 watts (341 Btu/hr). The board temperature is therefore 
calculated as follows, assuming that the board is 17 x 25 inches. 

Q - hA(T hoarA - r interi0f ) 

341 - (0.76)(-~^^)(7' board - 242) 

152 = (7 board -242) 

Aboard = 242 + 152 = 394 

The reader should bear in mind that this is just the starting 
point in the design and not the final stage. Based on these num¬ 
bers, the designer then makes a decision on the viable approaches 
to thermal management. 

FIN DESIGN 

Fins improve the efficiency of convection cooling greatly. How¬ 
ever, there are a variety of fin designs and the designer still faces 
the issue of selecting the fins particularly for the external surfaces 
of the enclosure. In this segment, we will discuss a technique for a 
simple, first order calculation for the most basic shape of the fin, 
i.e., plate fins. 

Basic Procedure 

In this procedure, one must first evaluate the temperature 
distribution of the system without the presence of fins based on the 
techniques described previously and calculate a heat transfer coef¬ 
ficient. The next step is to increase the heat transfer coefficient to 
lower the critical/design temperature to the values set by the de¬ 
sign criterion. Then, based on the following formula, calculate the 
fin geometry. 

Original Area x Newly Calculated h = 

Fin Surface Area x Original h 

This formula is not an exact relationship but is useful for a 
quick and “back of the envelop” calculation. For a better explana- 
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tion of this approach, let us solve the previous example again and 
select a fin. 

RF Cabinet Free Convection Cooling 

The initial thermal characteristics of an RF cabinet was inves¬ 
tigated as a preliminary stage in the design cycle. Now, we need to 
select a heat sink to maintain a maximum interior temperature of 
82 °C (179.6 °F) assuming an ambient temperature of 50 °C 
(122°F). Recall that the RF cabinet dissipates 800 watts of energy 
and its physical dimensions are 12 x 18 x 26 inches. The surface of 
the enclosure is to be painted white. 

Analytical Approach 

For the sake of simplicity, the solution procedure will not be 
repeated here. Only equation 5.7 is rewritten here: 


2958.75 = (24.5)(Tsurface - 582) + heff (13.83)(T SU rface - 582) + 
21.32 x 10- 9 (T 1 2 * 4 surface " 582 4 ) 

The only difference is that h has been replaced with heff. The 
solution procedure is as before; however, there is no need to bal¬ 
ance heff against hcaicuiated- The results for a variety of hes are pre¬ 
sented in | table 5.31 Since, the design criterion is based on interior 
temperature, let us consider this portion of analysis. 


Calculation of Internal Temperature 

It will be assumed that 85% of the generated heat escapes to 
the surface through conduction. This may be achieved by mount¬ 
ing the heat generating elements directly onto the wall and use of 
copper straps on the inside. Thus, only 15% of the generated heat 
is distributed inside of the cabin. Following the same heat balance 
principle: 

1) Input 

Qeiect.conv- = 0.15 x 800 x 3.41 = 409.2 Btu/hr 

2) Output 

Qoutput = h A (Tintemal - T S urface) 


By equating 1 and 2, one obtains the following relationship: 


409.2 — h (13.83) (Tintemal - Tsurface) 
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The solution procedure is exactly the same as before and the 
values of h are balanced against the calculated h. Table 5.3 gives 
Tintemai and T sur f act . for a variety of heir values. 


Table 5.3. The heff, Internal and Surface Temperatures for a 122 °F 
_ Ambient Temperature _ 


heff 

Surface 

Temperature 

Interior 

Temperature 

Interior h 

1 

173.45 

216.72 

0.68 

2 

163.75 

207.02 

0.68 

3 

157.06 

200.33 

0.68 

4 

152.20 

195.47 

0.68 

5 

148.50 

191.77 

0.68 

6 

145.61 

188.88 

0.68 

7 

143.28 

186.55 

0.68 


Note that the value of Interior h has remained constant 
through the range of surface temperatures. Furthermore, increas¬ 
ing the exterior heat transfer coefficient has a diminishing return. 
And, even at the large value of he tr = 7, the interior temperatures 
are above the design criterion. What is a possible solution? Let us 
assume that we could increase the interior heat transfer coefficient 
by employing a fan, per se. How would this impact the design crite¬ 
rion? Table 5.4 presents the same results for Interior h = 1 and 2. 
Clearly, the design criterion could be achieved for Interior h = 2, 
and fiefr = 3. 


Table 5.4. The heff, Internal and Surface Temperatures for a 122 °F 


Ambient Tern] 

perature 

heB 

Interior h = 1 

Interior h = 2 

Surface 

Temp. 

Interior 

Temp. 

Surface 

Temp. 

Interior 

Temp. 

1 

173.45 

203.04 

173.45 

188.24 

2 

163.75 

193.34 

163.75 

178.54 

3 

157.06 

186.65 

157.06 

171.85 

4 

152.20 

181.79 

152.20 

166.99 

5 

148.50 

178.09 

148.50 

163.29 

6 

145.61 

175.20 

145.61 

160.40 

7 

143.28 

172.87 

143.28 

158.07 
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Fin Design 

It is important to be able to determine what fin size and 
spacing would produce the effective heat transfer coefficient calcu¬ 
lated above. Basically, for the same AT value: 

h Afin ~ heS Asides 


Suppose that the fin height 
is , the thickness is t, and the 
spacing is i 2 , tj is the fin 
efficiency and n is the number of 
fins and L is the total length to 
be covered by the heat sinks. 
The following relationship may 
be deducted: 



n = INT- 


h eff L , 
h(2( | +1 + / 2 )*1 


(5.8) 


Since, the fin geometry (i.e., l\, t, and l 2 ) is specified by the 

design engineer, these parameters must be subjected to the fol¬ 
lowing geometric constraint. 


L — 2 ^w(f + ^2) — L 


For this particular example, note that the natural convection 
coefficient of heat transfer, h, is 0.73 and the perimeter of the cabi¬ 
net is 42 inches. For the above values of hctr and a fin efficiency of 
0.9, the following fin size and spacing may be selected. 


Table 5.5. Proper Fin Spacing and Corresponding Thickness 


freff 

L 

(inches) 

'i 

(inches) 

t 

(inches) 

h 

(inches) 

n 

3 

42 

2 

0.225 

0.9 

37 

4 

42 

2.5 

0.2 

0.8 

42 

5 

42 

2.75 

0.125 

0.7 

50 
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Design Recommendations for This Example. 

It is quite possible for a 12 xl8 x26 in. box to dissipate 800 
watts of energy and maintain an 82 °C internal temperature. How¬ 
ever, one needs to exercise caution and careful analysis so that lo¬ 
calized temperature rises do not exceed the maximum limit. This 
requires that all high power heat sources be attached directly to 
the internal surface of the external heat sink. All high power heat 
sources must be mounted with copper straps. An internal low ve¬ 
locity fan must be incorporated so that localized hot spots do not 
develop. 

Fin Design Considerations 

A total of 37 fins with the following characteristics need to be 
utilized. 


The fin must be 2 inches tall. 

The fin must be 0.225 inches thick. 

The fin spacing must be 0.9 inches. 

A More Exact Procedure 

A more exact treatment of fin heat transfer is beyond the 
scope of this book. The primary reason is that on the one hand, the 
developed flow may become quite complex depending on the ge¬ 
ometry, and on the other hand, radiation may play a strong role in 
the heat transfer. The reader is encouraged to read reference 11. 

FORCED CONVECTION 

While there are a variety of forced-convection cooling tech¬ 
niques, these may be classified into three types; namely, Direct 
Flow, Cold Plates and Heat Pipes. In Direct Flow, forcing air or 
other coolants directly over components cools the PCB. Cold Plates 
cool PCBs are indirectly by either air or other fluids. Heat Pipes 
take advantage of phase change properties of certain materials to 
transfer heat from one place to another. 
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Air Flow Out 


Figure 5.3. In direct flow, the components are cooled directly. 


In Direct Flow cooling, the designer must be concerned with 
the following issues before choosing a fan or a pump: pressure 
drop/losses, flow paths, component spacing and distribution, and 
in case of air, effects of altitude. 



Coolant 


Figure 5.4. In Cold Plates, flow passes behind the PCBs. 

Cold plate design is on the basis that heat is transferred to 
coolant through a heat sink, therefore, it may be designed into the 
chassis. A cold plate design requires flow rate calculations to pre¬ 
vent either choking or lack of pressure. Furthermore, the designer 
must prevent coolant leakage. 

Heat pipes are designed to transfer heat from a hot spot to a 
location where it can easily be removed. They work best when heat 
is concentrated. They work on the basis of evaporation/ condensa¬ 
tion. Their performance (e.g., operating temperature) is based on 
wick as well as fluid media. We will not be concerned with heat 
pipe design or analysis issues in this book. 
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Figure 5.5. Heat pipes work best when heat source is concentrated in one 

area. 


DIRECT FLOW SYSTEM DESIGN 

To design an enclosure based on direct flow cooling principle 
is relatively simple. The first step is to determine the total heat dis¬ 
sipation and the second step is to establish the maximum allow¬ 
able coolant temperature rise at the exit. The third step is to cal¬ 
culate the flow rate that would maintain such a temperature rise 
for the given heat dissipation. 

These three steps are concerned with the entire system. In the 
forth step the designer must consider the heat distribution among 
various boards and the temperature rise of critical components; 
thus determining the required flow rate for each board. Based on 
this information, board and component spacing may be deter¬ 
mined. 

The last step is to develop the pressure drop curve for the 
system. The pressure drop curve also known as the Impedance 
curve enables us to find the fan or pump necessary for delivering 
the proper flow rate through our system. 

The Required Flow Rate 

Often, it is up to the design engineer to maintain a certain 
upper limit for the temperature rise. To do so, one has to evaluate 
the required flow rate. A simple one-dimensional relationship de¬ 
rived from the heat equation is as follows. 
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where m is the mass flow rate, Q is the dissipated heat, C p is 

the specific heat of the fluid, and AT is the temperature rise of the 
fluid form inlet to exit. 

We need to keep in mind that the end result of this procedure 
is to select the proper fan. Many fan catalogs provide fan charac¬ 
teristics based on volume flow rate such as cubic feet per minute 
(CFM). Furthermore, air is the coolant of choice for a large class of 
electronics enclosures, and so, it is only logical to rewrite this 
equation in these terms. 


pV = 


Q 


C p AT 


V = ■ 


pC p AT 


(5.10) 


Now consider the Equation of State for air (as a perfect gas) 

P a = pRT a , 

lb 

where P is expressed in ——, T in degrees Rankine, and R (1718 

ft 2 

ft lb 

- for air) is specific gas constant. Equation 5.10 may be fur- 

Sec R 

ther reduced to 

. QRT 0 

PoC p AT ' 

By substituting 0.24 and 1718, air’s specific heat and gas 
constant, respectively - as well as making sure that dimensions axe 
consistent - the following relationship may be developed. 

* 3.7 OT 

V - CFM (Cubic Feet per Minute) (5.11) 

Note that in this relationship, the following dimensions must be 
maintained. 
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Q is the power to be dissipated in Btu/hr, 

T a is the inlet temperature in degrees Rankine (°F+460), 

P is the inlet barometric pressure in lb/ft 1 2 3 

AT is the temperature rise across the equipment in degrees F. 

Sometimes barometric pressure is given in terms of inches of 
mercury (“Hg). In this case the flow rate equation may be rewritten 
as 


y 00524 QT a 

P.AT 


CFM . 


Now, the dimensions, which must be maintained, are as follows. 

Q is the power to be dissipated in Btu/hr, 

T a is the inlet temperature in degrees Rankine (°F+460), 

P is the inlet barometric pressure in “Hg 

AT is the temperature rise across the equipment in degrees F. 

Some industry standards such as ARINC 600 require that the 
flow rate be calculated per kilowatt of dissipated heat. For this rea¬ 
son, the same equation may be manipulated: 


y _ 178.8 QT 0 
P 0 AT 


CFM . 


Again, the dimensions, which must be maintained, are as follows. 

Q is the power to be dissipated in kW, 

T c is the inlet temperature in degrees Rankine (F+460), 

P is the inlet barometric pressure in “Hg 

AT is the temperature rise across the equipment in degrees F. 
Exercise 

An electronic box has the following components and maxi¬ 
mum potential heat dissipation distribution: 

1. PCB 1- uP Board, 26 watts 

2. PCB 2- Memory board, 26 watts 

3. PCB 3 - Network board, 19.5 watts 
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4. PCB 4 - Network board, 19.5 watts 

5. PCB 5 - Power Supply, 54 watts 

6. PCB 6 - El Network Board, 19.5 watts, quantity 2 

What is the needed airflow (in CFM) to have a inlet to exit 
temperature rise of 10 degrees C at Sea level condition? At Sea 
Level, standard temperature is 59°F and standard pressure is 29.9 
“Hg. 


T = 460 + 59 = 519°R 

Q = 26 + 26 +19.5 + 19.5 + 54 + 2 x 19.5 = 184 Watts 

Q = 184 (Watts) x 3.41 (Btu/hr)/(Watt) = 627.44 (Btu/hr) 

AT = 10°C = 18°F 

", (0.0524)(627.44)(519) 

(18)(29.9) 

V = 31.7 CFM 

So, the required flow rate to maintain a 10°C inlet to exit tem¬ 
perature rise is 31.7 cubic feet per minute. 

Board Spacing and Configurations 

The next step after calculating the required airflow is to spec¬ 
ify board spacing. For instance, what should the board spacing for 
the previous system be? To answer this question, we need to de¬ 
velop a better understanding of fluid flow and its relationship to 
geometry as well as changes in pressure or pressure losses. 

Flow Resistance 

Fluid flows from one point to another only and if only a pres¬ 
sure difference exists between these two points - much the same 
as heat flow between two points requires a temperature difference. 
And similarly, the path in between these points resists this flow. In 
the case of fluid flow, this resistance is due to friction and drag. 
The general governing equations are the Navier-Stokes equations, 
mentioned earlier, and an important subset of these is pipe flow 
equations. Piping systems and their associated flow equations 
closely approximate fluid flow through electronics enclosures. 

The relationship between the pressure drop and flow rate 
across a pipe of length L and an arbitrary cross section is 
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AP = N 


(5.12) 


where the resistance factor (or number of heads lost) is represented 

by 


N = f ~~~ , (5-13) 

U h 

f is the Darcy friction factor, L is pipe length and D h is hydraulic 
diameter. Darcy friction factor, a well studied subject for a variety 
of fluids and pipes, is generally inversely proportional to Reynolds 
Number. For PCB’s a typical value is about 0.02. For more details 
see Mark’s Standard Handbook for Mechanical Engineers (12). 

• 

Equation (5.12) may be rewritten in terms of flow rate (V )as 
follows: 


AF = N 


. 2 ^ 

pv 


(5.14) 


In this work, we are primarily concerned with air as the fluid, 
pressure units in inches of water ("H 2 O) and area expressed in 
terms of squared inches. To reflect these concerns, equation (5.14) 
may be reduced to 


crAP = £(pV) 2 "H 2 0, (5.15) 

where loss coefficient (or factor) is defined as follows. 

N 

& - 0 . 226 — 

A 2 


"H 2 Q 

(%-n ) 2 


(5.16) 


Note that V is the flow rate in cubic feet per minute (CFM), p 
is air density at fan inlet (lb/ft 3 ), a is the ratio of average density to 


Copyright © 2003 by Marcel Dekker, Inc. All Rights Reserved. 



standard density (0.0756 lb/ft 3 ) and A is the minimum flow area 
(in 2 ). 

Number of Heads Lost (N) 

It is obvious that evaluation of the temperature drop depends 
on knowing the number of heads lost. Good sources for this are 
references (5, 12, 13 and 14). Some typical values are shown in 
Table 5.6. 


Table 5.6. Velocity Head Losses in Electronics Enclosures 


Configuration 

Inlet 

90°tum 

PCB’s 

Outlet 

N 

1.0 

1.5 

0.5 - 2.5 

1.0 


The number of heads lost for a PCB, may be best evaluated 
based on the following formula if proper experimental data is not 
available. 


where typical values of /are between 0.01 and 0.02. L is the flow 
length and Dh is the hydraulic diameter based on open area. 

Flow Networks 

Earlier, we indicated that flow through piping systems might 
be adapted for flow through electronics packages. The elements of 
such a network are either in series or in parallel or both: 

Elements in Series: 

When the elements of a flow net work are set in series con¬ 
figuration, the total loss coefficient is the sum of the individual 
element losses and the total pressure loss is the sum of pressure 
losses in each segment. This may be expressed as follows. 

^ total = ^ 1 + ^ 2 + & 3 + • • • 

i 

m 

A^total = X + '' ‘ + 


Elements in Parallel: 

When the elements of a flow net work are set in a parallel 
configuration, the physical laws of fluid flow dictate that the pres- 
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sure losses be equal in all the branches. Based on this require¬ 
ment, the total loss coefficient may be developed as follows. 


total ' #7 #7 + #7 + 

^total ~ ^*1 ~ ^2 “ "' “ 

Flow Rate Distribution in Parallel Plates 

Now that we understand that pressure drop is directly related 
to flow rate and loss coefficient and indirectly to geometry, we can 
begin to develop a means of calculating board spacing based on 
available flow rate and required temperature rise. 

• 

Consider that the flow (at a rate of V) is first divided into two 
parallel channels and then joined again. It is clear that the sum of 
the flow rates in each channel is equal to the total flow rate. 




V = Vi + V 2 (5.17) 

Also, since the channels are parallel, the pressure drops in 
each channel are equal. 


AP, - &P 2 => itjfpVj ) 2 =I 2 (pV 2 ) 2 => 
^l(Vi) 2 =^ 2 (V2) 2 


Equation (5.17) along with (5.18) form a set of two equations 
for two unknowns. The solution to this set of equation is as follows. 


Vi = 


\- 


V 


and 


V 2 = 


1 


1 + 


*2. 

*1 


V 


Earlier, we raised the question of calculating the proper board 
spacing. Theoretically, one may first specify a board pitch, calcu¬ 
late the loss coefficient based on that value, and based on the total 
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required flow, determine the flow rate in each channel. Should this 
channel flow rate prove to be insufficient, board pitch may be al¬ 
tered. The reader may develop the relationships for three or more 
channels. 

Optimum Board Spacing and Heat flow 

Similar to the natural convection case, Adrian Bejan, et al, 
(10) provides the means of calculating an optimum gap spacing for 
forced flow and the associated maximum heat removal: 


^ t =2.7inf 

Q max =0MAT^nY 1 2 3 4 


where 

n L = 


^ API? ^ 
pa 


H is the transverse length of the entire package and L is the length 
of the boards along the flow, W is the width of the stack, a is ther¬ 
mal diffusivity, g is viscosity and v is kinematic viscosity (= 
g/p where p is density). AP is pressure drop across the flow. 

System’s Impedance Curve 

Earlier, it was mentioned that the last step before specifying a 
fan or a pump is to develop the impedance curve of the system. 
This is really a fancy way of saying to plot the pressure drop as a 
function of flow rate. In fact, we already know that the pressure 
drop is a quadratic function of flow rate as follows. 


cjAP - to t al (pV) 2 "H 2 0 

So, in a way, to develop the impedance curve of a system is to 
find .R tota |. The following are the steps required for AP Calcula¬ 
tions: 

1. Draw the flow network through the enclosure. 

2. Assign numbers to each section. 

3. Evaluate each segment’s number of heads lost ( N ). 

4. Evaluate the minimum flow area (A). 
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5. Calculate loss factor (it) for each segment. 

6. Calculate the total enclosure loss factor (it total ). 

7. Develop crAP table vs. V (flow). 


Sample Problem 


Consider the system in the previous exercise shown in figure 


5.6. The system dissipates a total of 184 Watts. We need to develop 


the system’s impedance in order to select a proper fan. The oper¬ 
ating conditions are between sea level and 7000 feet. Recall that 
the maximum potential heat dissipation distribution is as follows. 


1. PCB 1- uP Board, 26 watts 

2. PCB 2- Memory board, 26 watts 

3. PCB 3 - Network board, 19.5 watts 

4. PCB 4 - Network board, 19.5 watts 

5. PCB 5 - Power Supply, 54 watts 

6. PCB 6 - El Network Board, 19.5 watts, quantity 2 


Each board is 6.5 x 10 in. and minimum clearance for flow 
between boards is 0.25 inches. To develop the flow diagram, we 
need the geometry information for the following stations. 


Station 1: inlet contraction, area is 
nr 2 = (3.14)(1.5/2) 2 

Station 2: short pipe, area is 
nr 2 = (3.14)(1.5/2) 2 

Station 3: 90 degree turn, area is 
nr 1 = (3.14)(1.5/2) 2 

Station 4: expansion transition pipe, area is 

Average the two cross - sectional areas = (1.77 + 26.25)/2 

Station 5: Perforated surface, area is 
for each hole 7tr 2 = (3.14)(0.0.075) 2 

Station 6: expansion into enclosure, notice that air expands into 
area =(10)(4.8) 

Station 7: Flow through PCB’s, 

minimum flow passage area for each PCB is = (10)(0.25) 

Station 8: Outlet Perforated surface, 
for each hole nr 2 = (3.14)(0.15/2) 2 
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Table 5.7. Each Station’s Loss Factor 


Station 

Area (in 2 ) 

N 

.ft 


^ For Station 

1 

1.77 

Li 

7.24E-2 


0.0724 

2 



1.45E-3 


0.0014 

3 

1.77 


1.09E-1 


0.1086 

4 


■ 

1.73E-3 


0.0017 

5 

0.018 

■ 

1.05E+3 

. 1 . - V 1 

Station 450 5 

0.0052 

6 

48 

l 

9.81E-5 


0.0001 

7 

2.5 

l 

3.62E-2 

1 -V 1 - 

Station 6 7 

0.0010 

8 

0.018 

1.5 

1.05E+3 

- 1 . -- V 1 

V -ft Station 450^8 

0.0052 

System’s Total Loss Factor 

0.1955 


Now we can develop the impedance curve using a density of 
0.0765 lb/ft 3 leading to a =1. The impedance curve for this system 
at sea level is governed by the following relationship. 

AP = 0.1955(0.0765 V) 2 
AP = 1.14x10~ 3 (V) 2 "H 2 0 

At 7000 feet, the standard atmospheric conditions are as follows: 

T =493.73 °R 
P =1.6331xl0 3 lb/ft 2 
p =0.0620 lb/ft 3 


Now we can develop the impedance curve using a density of 0.062 
lb/ft 3 . This leads to a = 0.81046 
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0.81046 AP = 0.1955(0.062 V) 2 
AP = 0.927x10 _3 (V) 2 "H 2 0 

Figure (5.7) depicts the impedance cure of this electronic system at 
sea level and 7000 feet. 



0 10 20 30 40 50 


■ CFM | 

Figure 5.7. The System’s Impedance Curves at Sea Level and 7000 Feet 

It is noteworthy to consider that the number of heads lost for 
the PCB section was set to one (1). This resulted in a loss coeffi¬ 
cient of 0.1955 and a maximum sea level pressure equal to 2.86 
inches of water. Had we set that value to 2.5, the loss coefficient 
would have been 0.1970 and the maximum pressure, 2.88 inches 
of water. 

Fan Selection and Fan Laws 

From the point of view of this book, the issues involved in fan 
selection center around delivering proper flow rate so that the 
needed temperature rise may be maintained. However, there are 
other selection issues that need to be considered. These are 

1. Noise; mechanical and/or aerodynamic 

2. Reliability 

3. Power 

4. Size and Shape 

A given fan can only deliver one flow rate for a given system 
impedance. Earlier, we developed a technique for evaluating the 
required flow rate for a given temperature rise (AP). We also know 
how to develop the impedance curve for a system. The task of fan 
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selection is to find a fan with a performance curve such that it in¬ 
tersects with the impedance curve at the needed flow rate. Let us 
clarify this by looking at the performance curve in more detail. 


Fan Performance Curve 


Fan performance curve provides the pressures loss values for 
various delivered flow rates. The fan manufacturer generally pro¬ 
vides this data. The operating point of a fan is the intersec tion of 


performance and impedance curves as shown in Figure 5.8. 


Fan Laws 

Some times a fan evaluated for one set. of conditions has to be 
used under a different set of conditions. Fan laws allow us to ex¬ 
trapolate current data to the new operating environment. 


V A = V B 


2 

^ Size a 
S ize R 


" rpm a " 

^RPM b> 


Pressure 


Power A 


A = Pressure B 


Size. 


Size 


B 


'( RPM, 


l 


rpm b 


= Power B 


Size; 

V Size B j 


Y RPM a 

^RPM b 


Pa 
Pb j 


r 

Pa 

kPb 


(5.19) 

(5.20) 

(5.21) 


In these equations, Size is fan diameter, RPM is the revolution 
per minute, and Power is the required power to operate the fan. As 
an example, these equations indicate that all factors remaining the 
same, should we double the size of a fan, its flow rate would in¬ 
crease by a factor of 8, its pressure drop by a factor of 4 and power 
consumption by a factor of 32. 


Copyright © 2003 by Marcel Dekker, Inc. All Rights Reserved. 

























Figure 5.8. Various Fan Performance Curves Intersecting One Impedance 
Curve. The intersection determines the delivered flow rate. 

Fan Networks 

Many Systems Require fans to be placed either in parallel or 
serial configurations. The impact of these configurations on the 
flow needs to be understood. In a parallel configuration, the flow 
rate increases but the maximum pressure drop remains the same. 
However, in a serial configuration, the pressure drop increases but 
the maximum flow rate remains the same. In practice, fan trays, 
i.e., two or three fans in parallel, are used to increase the overall 
delivered flow rate to a system without increasing the size of a sin¬ 
gle fan. There is an added benefit as well. By having several vari¬ 
able speed fans in parallel, should one fails, the others may be de¬ 
signed to increase speed to compensate for the failed unit. Fans in 
series cure generally used to increase pressure such as a compres¬ 
sor. 
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Fans in Parallel 



Figure 5.9. Two Identical Parallel Fans 



Component Hot Spot 

The sole purpose of undertaking all the previous assumptions 
and calculations is to determine the temperature of the hottest 
spot in the system and the associated component. It will be shown 
later that reliability of a component is degraded as its temperature 
increases. The two major contributors to this hot spot temperature 
are: 


1. heat accumulated by the coolant, 

2. resistance of the film coefficient to cany heat away from 

the component. 

Based on these two factors, the hot spot temperature is: 
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^spot ^inlet ^^coolant 

The design engineer usually specifies Ar coo ] ant , the maximum 
temperature gain of the coolant. At times, a particular standard or 
code prescribes this value. A7fj] m , however, must be calculated. 
Recall that 


Q = h c AM mm . 


In practice, this equation is illusive and at times difficult to solve 
because /^varies depending on geometry, coolant type, and flow 
regime. Thus its exact value may not be readily determined. It is in 
fact the subject of much discussion and research. However, for de¬ 
sign purposes and in engineering practice, compromises and edu¬ 
cated assumptions must be made. 

A stack of plug-in PCB’s may be assumed to behave very close 
to ducts. Therefore, the following formula may be applied. 


h c =JC p G\ 


(CpM\ 

v * , 


-2/3 


where 


G is weight flow rate per unit area = 60 


( . > 

pv 

A total 
V J 


• 

V = flow rate (CFM) 

J = Colburn Factor 

C p = Specific Heat of Fluid 

H = Fluid Viscosity 
K = Thermal Conductivity 
p = Fluid Density 

total = total flow cross-sectional area (ft 2 ) 


(5.22) 


Steinberg (5) provides the following values for Colburn factor 


200 < Re < 1800 (Laminar Flow) Aspect Ratio >8, J - —r— , 

Re 0 ' 98 
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200 < Re < 1800 (Laminar Flow) Aspect Ratio ~ 1, J - 


2.7 


Re 


0.95 ’ 


10 4 < Re < 10 5 (Turbulent Flows) J = ^ 


Re 


0.2 


Where, 


D GD A 

Re =- and, 

D is the the hydraulic diameter = 


4 x Area 
perimeter 


Exercise 

Recall the Cabin Telecommunication Unit (CTU) in the previ¬ 
ous exercise. It potentially has the following heat dissipation distri¬ 
bution: 


1. PCB 1- uP Board, 26 watts 

2. PCB 2- Memory board, 26 watts 

3. PCB 3 - Network board, 19.5 watts 

4. PCB 4 - Network board, 19.5 watts 

5. PCB 5 - Power Supply, 54 watts 

6. PCB 6 - El Network Board, 19.5 watts, quantity 2 

Earlier, it was determined that the required flow rate to 
maintain a 10°C (18°F) inlet to exit temperature rise is 31.7 cubic 
feet per minute. Furthermore, we learned that each board is 6.5 x 
10 inches and minimum clearance for flow between boards is 0.25 
inches. In this example, there are seven parallel boards but 6 ducts 
for air passage. 

Now, we should determine the Hot Spot Temperature. Assume 
the following conditions for Air at 150 °F, 

C n = 0.24-^- , p = 0.050— , K = 0.017 BtU - and 
v lb °F Ahr hr ft °F 

lb 

the inlet air density p - 0.0765——. 

ft 3 


Recall that 
Q = A4A7 fllm 
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Heat flow ( Q ) and surface area ( A ) are known. To calculate ATfi im , 
we need to find h from equation 5.22. 


h 


= JC p G\ 


C p M 

K 


2 

3 


For this equation weight flow rate ( G ) and Colburn factor ( J ) must 
be calculated. 


G = 


P v 


Total Cross Sectional Area 0.25 x 10 


(0.0765X31.7) =m8 _lb. 


144 


x 6 


min fC 


change the time scale to hours 
lb 


G =1396.83- 


hr ft' 


The choice of the equation for Colburn factor ( J ) depends on 
the value of Reynolds number. 

GD 

Re =-where D is hydraulic diameter, 

P 

'0.25x10' 


D = ■ 


144 


10 + 0.25 
12 


D = 0.0406 ft 


Re = (1396.83X0.0406), 

0.05 

The aspect ratio = 10/0.25 = 40, thus based on the values of 
Re and aspect ratio, the following relationship for Colburn factor 
holds. 


J - ■ 


Re 


0.98 


= 0.00609, 


Now the heat transfer coefficient may be calculated. 
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h = (0.00609)(0.24)(1396.83)^ 


(0.24)(Q.05) 

0.017 


-2 

3 


2.57 


Btu 

hr ft 2 F 


The highest heat dissipation of a board (PCB 5) is 54 watts 
(=54 x 3.41 = 184.14 Btu/hr). It may be assumed that heat is dis¬ 
sipated on both side of this PCB. Thus we have 


A T. 


184.14 


film,PCB5 


- 6.5x10,. 
2x~-(2.57) 


= 79.36° F, 


144 


Component = 59 +18 + 79.36 = 156.36° F(or 69.0°C) 


The two El network PCB’s (19.5 watts) have been placed on 
either side of the enclosure, so heat dissipation is done via one 
surface only. 


ATf ,| m>PCB6 - 


66.495 


6.5x10 

144 


(2.57) 


57.32° F, 


Component = 59 +18 + 57.32 = 134.32° F(or 56.8°C) 


INDIRECT FLOW SYSTEM DESIGN 

In the previous segment, steps needed to analyze and design a 
cooling system for electronics equipment were out-lined and ex¬ 
plained. The basis of that approach is that the cooling medium; 
namely air, is in direct contact with the components. Hence, the 
term direct flow system was used. Electronic Equipment also take 
advantage of heat exchangers and cold plates. Since these devices 
separate the electronics components from the coolant, they may be 
termed as indirect flow systems. While from a design point of view, 
direct and in-direct flow systems with their unique design issues 
are essentially different, from a physical scientific point of view, 
they are very similar. They both involve understanding pressure 
losses, fan/pump specifications and mainly maintaining the criti¬ 
cal component temperature below a design value. Depending on 
the heat dissipation levels, the in-direct flow systems may employ 
fluids such as water, oils and in the case of some aircrafts even 
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fuel. Therefore, design consideration not only should include 
choking but also temperature below the flash point of the fluids 
used. The steps of designing a cold plate are very similar to the 
method outlined earlier and the reader is encouraged to consult 
other works (such as reference 14). 

The relationship for flow rate remains the same, however, one 
must make sure to use appropriate values of density and specific 
heat. 


V = 


Q 

pC p AT ’ 


As before, Q is the dissipated heat, C p is the specific heat of the 

fluid, and AT is the temperature rise of the fluid form inlet to exit. 
Similarly, the relationship between the pressure drop and flow rate 
across a pipe of length L and an arbitrary cross section is 


AP = N 


V 2A 


g 


where the resistance factor (or number of heads lost) is represented 
by 


/ is the Darcy friction factor, L is pipe length and D h is hydraulic 
diameter. Fluid flow along with Darcy friction factor in pipes and 
ducts is well for a variety of fluids and pipes. In fact, the values 
presented previously are based on these theories. For more details 
see Mark’s Standard Handbook for Mechanical Engineers (12). 

Resistance Network Representation 

In the forgoing discussion, only pure convection was consid¬ 
ered. In many practical problem, conduction and radiation may 
exist and their contribution must be considered. Earlier, we stud¬ 
ied how to develop thermal resistance networks for conduction 
problems. The same concept can be extended for convection as 
well. As before 
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A7 7 -QR , where R = — 
hA 


Note that the difference with conduction thermal resistance is 
that K/L- the ratio of thermal conductivity to thermal path length 
- is replaced with h - heat transfer coefficient. 
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6 

Basics of Shock and Vibration 


So far, we have been concerned with the thermal performance 
of electronics enclosures as a whole. There are other issues in their 
design that can be as (or even more) important. Mechanical and 
electromagnetic issues are two major design factors that must be 
considered. In general, once the enclosure is designed and the 
prototypes are developed, a few samples are tested in environ¬ 
mental chambers as well as vibration (shake) tables for compliance 
with various standards. Designs are, then, only modified to pass 
the given tests and the designers, in general, know nothing about 
the behavior of their system in the field and any relationship that 
their data has to either failure rates or repair/maintenance sched¬ 
uling. 

In this segment, we will review another factor involved in me¬ 
chanical analysis; namely, vibration management. The term man¬ 
agement is used because in the majority of times, these issues 
must be dealt with as opposed to be just “analyzed” or “designed.” 
Sources of vibration may be categorized as follows: 

> In Stationary Systems 

■ Unbalanced Loads 

> In Road Vehicles 

■ Rough Surfaces Of The Roads 

> In Sea Vehicles 

■ Fluid/Structure Interaction 

> In Air Vehicles 

■ Aerodynamic Loads 

The study of vibration is concerned with the oscillatory mo¬ 
tions of bodies and the forces associated with them. All bodies pos¬ 
sessing mass and elasticity are capable of vibration. 
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Vibration can be broadly characterized as linear or nonlinear, 
free or forced. In free vibration, the system will vibrate at one or 
more of its natural frequencies, which are properties of the dy¬ 
namic systems established by its mass and stiffness distribution. 
In forced vibration, if the excitation is oscillatory, the system is 
forced to vibrate at the excitation frequency. If this frequency coin¬ 
cides with one of the natural frequencies of the system, resonance 
occurs. Resonance may potentially lead to large oscillations caus¬ 
ing catastrophic failures. 

It is noteworthy to consider that all vibrating systems exhibit 
damping because friction and other factors dissipate energy. 

Harmonic Motion 

When an oscillatory motion repeats itself regularly in equal 
intervals of time, it is called periodic motion. The time x is called 
the period of the oscillation and its reciprocal f=l/x is called the 
frequency. The simplest form of periodic motion is harmonic mo¬ 
tion depicted as 


x = A Sin(2n —), 
x 

where A is amplitude of oscillation, r is period of oscillation 
and t is time. It is customary to write 

2n -w 
co = —- = irg 

x 


Angular velocity, co , is generally measured in radians per 
second, and is referred to as the circular frequency. Furthermore, 
x and / are the period and frequency of the harmonic motion, usu¬ 
ally measured in seconds and cycles per second respectively. 

Let us look at velocity and acceleration of a point whose dis¬ 
placement is governed by x: 


Displacement: 

Velocity: 

Accelation: 


x = ASin(cot) 


dx _ 

~dt ~ 

,2 

d x 
dt 


x, t - coACos(cot) 

2 

- x, tt - -co ASin(cot) 
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Thus by combining these equations, the differential equation 
describing the motion of a single degree of freedom system may be 
obtained: 

co 2 x + x, tt = 0. (6.1) 

Periodic Motion 

It is quite common for vibration of several different frequen¬ 
cies to exist at the same time. Consider a guitar string or a multi¬ 
ple-degree of freedom system. The vibration of each natural fre¬ 
quency contributes to the free vibration and one may see a re¬ 
sponse as follows. 



Figure 6.1 . Response of A Multiple-Degree of Freedom System 


It is quite possible to assemble these complex forms once the 
contributing frequencies and mode shapes are known. This ap¬ 
proach is called Inverse Fourier Transform. The discussion of this 
subject is beyond the scope of this talk. 

Vibration Terminology 

Peak Value: It generally indicates the maximum stress, which 
the vibrating part is undergoing. It also places a limit on the "rattle 
space" requirement. 

Average Value: The average value indicates a steady or, static 
value somewhat like DC level of an electrical current. It can be 
evaluated by: 
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x = lim — \x(t)dt 
T -> qo T q 

For example, the average value for a complete cycle of a sine 
wave [ A sin(t) ] is zero; whereas its average value for a half-cycle is: 

- A n 2A 

x = — f Sin(t)dt =- = 0.637A 

^0 ' T 

Mean Square Value: It is the average of the square values, in¬ 
tegrated over some time interval T. This term is generally associ¬ 
ated with the energy of the vibration. 

~T l T 2 

x = lim —\x{t)dt 
T — > oo T q 

For a simple harmonic system, this definition leads to 



Root Mean Square (rms): A common measure of vibration, it is 
defined as the square root of the mean square value. 


rms = 

Decibel: The decibel is a unit of measurement that is fre¬ 
quently used in vibration. It is defined in terms of a power ratio: 

P X 

db = 101og(—^-) = lOlog(-J-) 2 
P 2 x 2 

The second equation results from the fact that power is pro¬ 
portional to the square of the amplitude or voltage. 

db - 201og(—!-) 
x 2 
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Thus an amplifier with a voltage gain of 5 has a decibel gain of 14: 


db = 201og(5) = 14 

Octave: When the upper limit of a frequency range is twice its 
lower limit, the frequency span is said to be an octave. 

Free Vibration 

Often times, free vibration is demonstrated through the use of 
a spring-mass system as shown in Figure 6.2. Though this con¬ 
stitutes a very simple analogy, it leads to modeling simplifications 
for a large class of problems: 



A 

(A 

KD 


V 

7_ 


In Static Equilibrium 
KD=mg 




K(D+x) 


In Dynamic Equilibrium 
(Vibration) 

K(D+x)=mg 


Vme \7ma 

Figure 6.2. A Single Degree of Freedom System 

For the dynamic equilibrium, the forces may be balanced as 
follows. 


ma - mg - K(D + jc) or 
mx, tt - -kx 


( 6 . 2 ) 


Now we define co = — and obtain the following equation: 
m 
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x, tt +a> x = 0. 

Thus the motion is harmonic. The natural period of the oscillation 
and frequency are 

t - and / = —. . 

V K 2n\m 


First Application 

There is a point to all this! Consider that any physical system 
deflects due to an applied load. Therefore, theoretically, it is possi¬ 
ble to develop a force-deflection relationship for any physical sys¬ 
tem under applied loads. Now suppose that we limit our interest in 
only one location in the system and in one applied load. These ex¬ 
pressions are available for many well-defined engineering systems 
such as beams, plates and shells. See reference (12) as an exam¬ 
ple. This in effect is a reduction of a complicated system into a 
simple spring-mass model. 

As an example of this technique, let us calculate the natural 
frequency of a 10 pound block at the end of a cantilever beam as 
shown in Figure 6.3. Neglect the mass of the beam. 



Based on strength of material assumptions, the expression for 
the maximum deflection of the free end of a cantilever beam due to 
a concentrated load at that end is 
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where P is the load at the end of the beam, L is beam length, E is 
the Young’s modulus of elasticity and I is area moment of inertia. 
By rewriting this equation, we obtain 


" 3 El 

, I? 


\S = P. 


Notice this has the same format as a spring-mass force deflection 
relationship (i.e. Kx=F), Where 


K > ——, x > S and F > P(- mg) 

L 3 

Therefore, the expression for the natural frequency may be 
employed: 



Suppose that the beam is made of aluminum (E =10xl0 6 psi) 
and has a moment of inertial / = 0.03 in 4 . The equivalent mass for 
the 10 lb. force is 


10 lb-Sec 2 

m = -= 0.31 


32.2 


ft 


After substituting appropriate values, the natural frequency is cal¬ 
culated as 


/ = 6.5 or = 41 . 


This equation is valid only for a cantilever beam. For other 
types of beams appropriate force-deflection relationship must be 
used. Clearly, equations developed in this manner provide natural 
frequencies of equivalent spring-mass systems but the actual 
shape of the system (i.e. the vibration mode shape) is not readily 
available through this approach. 

A PCB behaves very similar to a spring and mass system for 
its fundamental mode of vibration. It is therefore, theoretically pos¬ 
sible that based on this approach, one would use the force- 
deflection relationships available in various handbooks to calculate 
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this fundamental frequency. Depending on the geometry and con¬ 
figurations, there are two ways to make this approximation: 
namely, based on Beam Deflection or Plate Deflection. 


An Example for Comparison 

Consider a large (12 x 6 x 0.08 in.) PCB (glass-epoxy with 
copper planes, E=3xl0 6 , v=0.18,) weighs 1.5 pounds. Find the 
natural frequency of the system for all sides hinges and the short 
sides clamped. 

The natural frequency of a plate of length a and width b (a>b) 
is 


Simply-Supported (hinged) Plate: f n = — 


Clamped-Simple Plate: f n = 


Va 2 + b 2 


M 

1/2/ 

Ip J 

V 


f-f 2 

sP, 

16 

2 + 2,2 + ,2 

a a b b 


8 


1/2 


The deflection at the center of the plate is 


wb 4 

Simply-Supported (hinged) Plate: y = 0.0838—— 

Et 3 

wb 4 

Clamped-Simple Plate: y = 0.0582——. 

Et 3 

Where E is Young's Modulus, t is plate thickness, w is weight per 
unit area, p is mass density per unit area, D is flexural rigidity and 

v is Poisson's Ratio. 


D = 


Et 3 

12(1 -v 2 ) ' 


The frequency of the equivalent spring mass system is calcu¬ 
lated from 




1/2 


The frequency of the plate is given in Table 6.1.| While, this 


example shows very good comparisons, in general it is better to 
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calculate the natural frequency using much more robust tech¬ 
niques. 


Table 6.1. A Comparison of Frequency Values 



Exact 

1 DOF 

% error 

SS frequency 

98.10 

94.08 

4.10 

CS frequency 

111.12 

112.88 

-1.59 


Second Application 

In many vibration problems, a system undergoes a given ac¬ 
celeration level (e.g. 2g’s). We may take advantage of this informa¬ 
tion to develop a means of calculating the maximum deflection that 
the vibrating part undergoes. Recall that for a simple harmonic vi¬ 
bration, we have 

x - X 0 Sin(co t) 

■y 

x, tt = -co X 0 Sin(a>t) 

X 0 is the maximum deflection of the mass, and x, tt is the ac¬ 
celeration. Clearly the maximum acceleration is 

x, tt = co 2 X 0 

Note that the negative sign is dropped because we are interested in 
the magnitude not the direction of acceleration; therefore, the 
maximum deflection is: 


X a = 


x, tt 


a) 


This equation in and of itself is not useful, but as mentioned 
earlier, in many cases, an acceleration level is specified . For exam¬ 
ple, it may be that the cantilever system in i figure 6.31 undergoes a 
2-g acceleration. Thus, the maximum dynamic displacement is 


= 


2(32.2) 


(0 


Recall that co = 41. Therefore, the maximum deflections is 
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X, 


2(32.2) _ q Q^g f eet or 0.456 inches. 

41 2 

In a later chapter, we will discuss how dynamic loads may be 
evaluated based on the maximum dynamic displacements. These 
loads are needed to calculate the dynamics stresses used to evalu¬ 
ate the enclosure vibration worthiness. For now, it will suffice to 
say that X 0 is equivalent to 5. One would obtain 

T.FJ y 

Dynamic Load P = ——-p- 
(o 2 L 3 

This will be discussed in more details in the Induced Stresses 
segment of this chapter. 

Mode Shapes 

A mode shape is the deformation of a system corresponding to 
a particular frequency. While, it is possible to calculate the fre¬ 
quency of a system, only the relative shape of the corresponding 
mode can be calculated. The exact magnitude is only obtained after 
the equation of motion is integrated (or solved for every time point). 
This is not a trivial task. To illustrated the reason behind this in¬ 
ability, consider equation (6.2) for a spring-mass system: 

mx, tt = -kx 

or 

mx, tt +kx =0 
Recall that 


2 

x = ASin(o)t)andx, (( -~o) ASin(cot) 

Now substitute these relationships into (6.2) 

-map" A Sin(co t) + kA Sin(o> t) = 0 
This relationship may be further reduced to 
(- mo ) 2 + k)ASin(cot) - 0 
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In here, product of three quantities; namely, (-mco 2 + k) , A and 
Sin(cot ) is zero. A is defined as the amplitude of vibration in a sin- 
gle-degree-of-freedom system (i.e., a spring-mass system) and thus 
corresponds to the mode shape. If A were to be set to zero, it 
would imply that the amplitude of vibration is zero. So A can not 
be zero and similarly Sin(cot) is non trivial. Therefore, (-mco 2 + k) 

must be zero leading to frequency calculation. Thus the exact mode 
shape of vibration may not be defined by any equations. 

This exercise illustrates how the frequency of a spring mass 
system may be obtained from a purely mathematical approach but 
this is obtained at the cost of not having any means of calculating 
A ; the vibration amplitude. 

Nowadays, most engineers employ finite element methods to 
determine the natural frequencies and mode shapes of complex 
systems. Generally, the mode shapes are stored in the same loca¬ 
tion as actual displacements. Thus, it is quite easy to imagine that 
actual displacements have been calculated. This is an oversight. 

The box below is an avionics communication computer sys¬ 
tem. The frequency of vibration is 60 Hz and the general shape is 
shown but the magnitude of displacement is not calculated. Notice 
how some of the PCBs have crossed each other. This obviously 
does not happen in real life but clapping of boards may indeed 
happen. The mode shape calculation enables us to pinpoint prob¬ 
lem areas. 



Figure 6.4. The Mode Shape of an Avionics Electronics Box 
Viscously Damped Vibration 

Damping acts as a force, which opposes motion. When a lin¬ 
ear system is excited, its response will depend on the type of exci¬ 
tation and the damping present in the system. For an undamped 
system, the amplitude of vibration is constant and does not change 
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with time, and its natural frequency is a function of stiffness and 
mass. Damping affects both the amplitude as well as the natural 
frequency. 

The governing equation is for a one-degree of freedom spring 
mass system is: 

+kx-F(t ) 

where represents the damping force. There are two classes 
of damping forces; one depends on friction called Columb damping 
and the other depends on velocity and is called viscous damping. 
Columb damping depends on surface properties and the level of 
applied forced normal to the surface and is difficult to quantify. 
Viscous damping, however, is expressed as 

fd ~ F x, t 

Thus, the equations of motion becomes: 

mx, tt +Cx, t +kx = F(t ) (6.3) 

Clearly, for a damped, free-vibration system, F(t) = 0 

Impact of Damping 

The solution to equation (6.3) when F{t) - 0 (i.e., no applied 
forces) is as follows 

x{t) = Ae at 

by substituting this solution into the governing equation, and ap¬ 
plying the boundary and initial conditions, coefficients A and a 
may be determined: 

a\ 2 - —~ ± Va where A = I — 

’ 2m V m 

On the one hand, A may be less than zero leading to imagi¬ 
nary values of a . Thus, the system is called under-damped ; the 
solution is oscillatory and vibration exists. On the other hand, 
A may be greater than zero leading to real values of a leading to a 
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non-oscillatoiy solution and vibration does not exist; the system is 
called over-damped. At the interface of these two conditions A = 0 
and this condition is referred to as critical damping-, however, vi¬ 
brations do not exist. This is depicted in Figure 6.5. 



Figure 6.5. Impact of Damping on Vibration 

Vibration only exists when the system is under-damped. How¬ 
ever, the amplitude of vibration as well as the system’s natural fre¬ 
quency are no longer constant. For a spring-mass system, the am¬ 
plitude of the oscillations is described as 


Xa = 


2 [ (v 0 +X 0 fo)) 


, 1/2 


co D 




x Q and v 0 are the initial position and velocity, respectively, co is 

C 


the undamped natural frequency and £ = 


2 m co 


is defined as the 


damping ratio. a> D the natural frequency of the damped system: 


co D 

Note that the natural damping present in most materials 
hardly exceeds 15 percent, i.e., £ = 0.15 and oj d =0.99 co . Therefore 
in evaluating natural frequencies and mode shapes - so long as 
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artificial damping does not exist - the impact of damping is ig¬ 
nored. 

Forced Vibration 

The general governing equation for a single degree of freedom 
system with damping is as follows: 

Mx, tt +Cx, t + Kx = F(t) 


If F(t) is harmonic, it may be expressed as F{t) = P 0 Sin{mt ), thus 
Mx, tl +Cx, t +Kx = P 0 Sin(tut ). 

The solution to this equation is a combination of general and par¬ 
ticular solutions: 


x{t) = {ACosm D t + BSinco D t)e ^ mt + 

.—— Sin(rnf - 0 ) 

A /(l-r 2 ) 2 +(2r^) 2 


(6.5) 


where 


X 0 = —, tan <9 = 

U jy > 


m 


K 


1 - r 


r - 


0) 


Equation (6.5) has the following characteristics: 


1. There is a start-up "transient segment" which diminishes 
exponentially. The decay rate depends linearly on both the fre¬ 
quency as well as damping. 

2. Under resonance conditions, i.e., the forcing frequency 
equal to the natural frequency, r =1, the magnitude of vibration is 


finite and approaches 



Engineering Applications 

From a design point of view, we are interested in a few guide¬ 
lines to assess the behavior of our system, particularly, as the ex¬ 
citation frequency nears the natural frequency, we need to know 
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Will the displacements be magnified? 

Will the excitation force be magnified? 

Dynamic Magnification Factor 

Dynamic magnification factor is a measure of dynamic dis¬ 
placements of a system compared to input or static values. It is an 
indicator of how displacements are magnified in a system vibra¬ 
tion. It is the ratio of steady state response amplitude to the static 
response amplitude. For a spring-mass system. 


D 


m 



_ 1 _ 

V(l-r 2 ) 2 +(2£-) 2 


( 6 . 6 ) 


where r = — is the frequency ratio and £ =-is the damping 

cd 2 Mco 

ratio. Note that D m depends only on g and r : 

1. For r « 1, D m s 1, i.e. X =X 0 , in other words, the prob¬ 
lem can be solved as a static problem. 

2. For r > 1.45, D m < 1. Thus when the frequency of the ex¬ 
ternal load is about 45% greater than the natural fre¬ 
quency, the amplitude of the steady-state response is less 
than the static displacement. In other words, the impact of 
dynamic forces are less severe than the static forces. 

0. For r= 1, D m = 1/(21;), D m is only a function of £. For a 

typical value of £ , say £ =0.01 (1% damping), D m - 50. Now 

assume that {a> was not calculated accurately, and co*-0.9 
cd , i. e., 10% error in calculation. The erroneous frequency 
leads to a magnification factor of nearly 5 instead of 50. 
Clearly, large errors in displacement may be introduced if 
the natural frequencies are not calculated accurately. 

Transmissibility 

Transmissibility is a measure of the magnification of input 
forces throughout the system It is defined as the ratio of transmit¬ 
ted force through the springs and damper to the amplitude of the 
applied force. 

For a single spring-mass system, 
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(6.7) 


T - F T- V 1 + ^) 2 

F o yj(\-r 2 ) 2 +(2 £r) 2 
Note that T m depends only on E, and r : 

1. For r « 1, T m = 1, i.e. the problem can be solved as a static 
problem. 

2. For r > 1.50, T m < 1. Thus when the frequency of the ex¬ 
ternal load is about 50% greater than the natural fre¬ 
quency, the amplitude of the transmitted force is less than 
the applied value. 

3. For r= 1, T m is only a function of £ . 



For a typical value of £, say £ = 0.01 (1% damping), T m = 50. 
Similar to Dynamic magnification factor a 10% error in frequency 
calculation leads to significantly large errors in the transmitted 
forces and resulting induced stresses. 



Figure 6.6. A Graphical Representation of Dynamic Magnification Factor 
or Transmissibility 
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Typical Transmissibility in Electronics Enclosures 

Steinberg (15) suggests that for an electronics system 

T= Ajfn 


f n is the fundamental frequency and the factor A varies de¬ 
pending on the particular configuration. 


1) Mounted brackets A=2 or 3 

2) Spring edge guides, A=0.5. 

3) PCB’s, A depends on the natural frequency: 

i) For fundamental frequency between 200 to 300 Hz, A= 1 

ii) For frequencies above 400 Hz, A=2. 

iii) For frequencies below 100 Hz, A=0.5. 


Sample Problem 

Recall the cantilever beam shown in 


figure 6.3 


Now assume 


that the 10 pound block is replaced with an oscillatory force of 10 
pounds with the following frequencies: m = 30, 40, 50 and 60. As¬ 
suming a damping ratio (£) of 1%, we like to calculate dynamic 


magnification factors and displacements as well as transmissibil- 
ities for this system. 



Figure 6.7. A Cantilever Beam Under Vibratory Loads 

Dynamic magnification factor and transmissibility are defined 
by equation 6.6 and 6.7 respectively, 


D 


m 



_i_ 

V(l-r 2 ) 2 +(2^) 2 
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, Ft Vl + (2^) 2 

F ° Vo - r 2 ) 2 + ( 2£-) 2 


where r - — and £ =-— . Static deflection was discussed earlier as 

co 2 Mat 

PI? 

= -, E = 10xl0 6 psi and I = 0.03 in 4 . So, 

3 El K 

S 0 = — = —• ()X( - 1 ^ ----= 0.0192 inches. 

3EI 3x10x10^x0.03 

Previously we calculated the natural frequency of this system to be 
co =41, and £ =0.01: 

for q = 30/41 = 0.73, thus 


F>m = -f-—. '" 2 - 14 and 

V(l-0.73 2 ) 2 + (2x0.01x0.73) 2 

S l - D m S a = 2.14 x .0192 = 0.041”. 


T = 

1 m 


y[\ + (2x0.73x0.01)" 


V< 


1-0.73 2 ) 2 + (2x0.01x0.73) 2 


2.14 


Similarly, 


r 

F>m 

<5 

T 

* m 

40/41 = 0.98 

22.6 

0.430” 

22.6 

50/41 = 1.22 


0.039” 

2.04 

60/41 = 1.46 


0.017” 

0.88 


Near resonance, this system will experience a theoretical dis¬ 
placement of 0.43 inches compared to 0.019 inches static dis¬ 
placement. In real life, the system may fall apart long before mag¬ 
nification factors of 22 are realized. Clearly, these vibrations need 
to be managed through selecting a proper isolator. 
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Vibration Isolation 

The function of an isolator is to reduce the magnitude of mo¬ 
tion transmitted from a vibrating foundation to the equipment or to 
reduce the magnitude of force transmitted from the equipment to 
its foundation. There are three types of isolation; namely, Fre¬ 
quency Isolation, Velocity Isolation, and Force Isolation. Further¬ 
more, there are four types of isolators; Rigidly Connected Viscous 
Damper, Rigidly Connected Coulomb Damper, Elastically Con¬ 
nected Viscous Damper, and Elastically Connected Coulomb 
Damper. The materials used in these isolators are as follows: 

> Metal springs used where large static deflections are re¬ 
quired or environmental conditions do not allow elasto¬ 
mers. 

> Elastomers providing highest levels of energy storage and 
dissipation with an added benefit of being molded to cus¬ 
tom shapes. 

> Plastic isolators being similar to elastomers but providing 
more rigidity. 

For more information the reader is referred to Shock and Vi¬ 
bration Handbook (16). The steps to select an isolator are as fol¬ 
lows. 

1. Select your desired magnification factor - in other words, 
what is the highest allowable displacements. 

2. Select the damping ratio. 

3. Calculate the spring stiffness. 

Example 1 

Find an appropriate spring to contain the magnification fac¬ 
tors below 0.5 for the cantilever beam shown in | figure 6.7.| Use 
30% damping. 

Notice that g =0.3 and that D m <%. Thus, 


^/(l-r 2 ) 2 +[2(0.3)r] 2 2 


which leads to (l-r 2 J +[2(0.3)r] 2 >4. The solution to this inequality 
is r >1.65. 
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This value is indicative of the new frequency ratio of the sys¬ 
tem with isolators. Recall that 



With g =32.2 ft/sec 2 and M =10/32.2 =0.31 slugs => co„ = 1.796^ 


1.65 =-===> K=182.2 lb./ft 

1.796VK 


This is the total stiffness. The stiffness of the cantilever beam is 


^bearn = — = 3 ( 1 Qx1q 6 X 0 - 03 ) = 520.83 lb/in = 6249.96 Ib/ft 
I? 12 3 


For Serial Springs: 


1 1 1 


K K K 

total beam isolators 


From here ^j so i ators = 187.67 lb./ft 


Example 2 

A 3.0 lb. electronic system can withstand a 5g sinusoidal vi¬ 
bration. The base will be loaded to 12g and frequency of 62 Hz. De¬ 
sign an isolation system with a damping factor of 20% for the 
mountings. 

This problem is concerned with forces and their magnification 
and/or decay in the system. Therefore, the transmissibility must 
be taken into calculations. Note that 20% damping means £ =0.2. 


l + [2/0.2>] 2 


5g 


}-r 2 J +[2/0.2 )rf l2g y(l-r 2 f +0.16r 2 


1 + 0.167-' 


_5_ 

12 


This expression leads to r 4 -2.762r 2 -4.76 = 0. The solution to this 
equation is r=1.99. Thus for r>1.99 the response of the system will 
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be less than 5g. Note that r is related to the spring stiffness 
through natural frequency. 

This value is indicative of the new frequency ration of the 
system with isolators. Recall that 


r - 


y- and /« =Y 

fn 



With g =386 in/sec 2 and M =3/386 =0.31 => f n = 1.805 Vk 


- —;=> 1.99 K< 297.8 lb./in 

1.805VK 

Note that the softer the spring, i.e., more isolation, the larger are 
the displacements 

Maximum PCB Deflection 

Earlier, it was demonstrated that the primary natural fre¬ 
quency of a plate may be calculated from the maximum deflection 
and vise versa, having the natural frequency enables one to calcu¬ 
late the maximum deflection. Earlier it was discussed that 


X =—or X a =— —., inches, (6.8) 

<y 2 (2 nf n ) 2 

where g=386 in/sec 2 is the gravitational constant and G is the 
level of acceleration experienced by the system. Clearly, G-G in T 
and T is the system transmissibility. Furthermore, we had previ¬ 
ously suggested that T =^/7T 5 thus equation (6.8) becomes 


x = 9.8G ( „ 


(fn) 


1.5 


(6.9) 


Alternatively, we may have a displacement limitation and 
therefore be interested in knowing the associated frequency for a 
given displacement: 

2 


fn 


9 .SG in 

x 0 


( 6 . 10 ) 
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Maximum Desired PCB Deflection 

Steinberg (15) suggests that if the peak single-amplitude dis¬ 
placement of the PCB is limited to X max , the component can 
achieve a fatigue life of about 10 million stress reversals in a sinu¬ 
soidal vibration environment. For sinusoidal vibrations 

^max-^^-ches (6.11) 


The desired corresponding frequency is given by f d . 


fd~ 


' 9.SG in ChrJEy3 

0.000225 

v 


( 6 . 12 ) 



Where b 

B is the length of the PCB edge parallel to component (inches) 
L is the length of the electronic component (inches) 
h is the thickness of the PCB (inches) 

C = 1.0 for Standard DIP 

= 1.26 for DIP with side-brazed lead wires 
= 1.26 for Pin Grid Array (PGA) with 2 parallel rows of wires 
= 1.0 for PGA with wires around the perimeter 
= 2.25 for leadless ceramic chip carriers 
= 1.0 for leaded chip carriers where the lead length is about 
standard DIP 

r =1.0 when the component is at the center of PCB 

= 0.707 when the component is @ 14 point X and !4 point Y 
on a PCB supported on 4 sides. 

= 0.5 when the component is @ l A point X and % point Y on 
a PCB supported on 4 sides. 
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Applications to Electronics Enclosures 

One has to bear in mind that vibration design and analysis 
does not end with PCB’s. In fact, the supporting structure(s) must 
be evaluated as well. The structure of the enclosure is made up of 
engineering components such Beams, frames and plates; each re¬ 
sponding to the vibration {or shock) excitations. Their response 
also impacts the well being of the entire electronics system. An¬ 
other aspect is in the mechanical design of the electronics package 
in of fasteners & clamps as well as components. Finally, one has to 
have a good understanding of the response of the PCB to make 
proper use of ribs so that a redesign does not become necessary. 

Figure (6.8) shows a corner of the supporting structure of a 
PCB containing relatively heavy transformers. This bracket fails the 
test conditions because of stresses that develop at the comer. A 
redesign is necessary for the system to remain functional under 
the test conditions. 



Figure 6.8. The Stresses Developing at the Comer of Supporting Structure 

Designing for Shock Environment 

The physics of shock environment is essentially different than 
that of vibration. Generally speaking, vibration is a steady state 
phenomenon whereas shock is transient; therefore, the absorbed 
energy must be released over a longer period of time than required 
by vibration isolators. Isolation for each has different requirements 
and if both must be considered, then compromises must be made. 

There are two basic types of shock environments; one is due 
to a sudden motion of a support or foundation, and the other is 
due to applied or generated forces by the equipment and experi¬ 
enced by the support. In an electronics packaging situations, the 
primary concern is with the former. 

There are three approaches to modeling a shock environment; 
namely, pulse shock, velocity shock and shock response spectrum. 
In a pulse shock model, accelerations or displacement shapes are 
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well known - generally a sinusoidal or square wave form. While 
analysis is quit easy using this approach, it is unrealistic. In a ve¬ 
locity shock, the system undergoes a sudden change of velocity. A 
common test is the drop test; where, equipment is dropped from a 
predetermined height and is expected to survive. The last tech¬ 
nique is called shock response spectrum where the structural re¬ 
sponse to motion as opposed to the motion itself is evaluated. It is 
the most accurate approach but at the same time, it is the most 
difficult as well. Practical applications of this techniques require 
use of numerical modeling and will not be discussed here. 

The steps for pulse and velocity shocks are demonstrated by 
the following two examples. 

A Pulse Shock Isolation Example 

The PCB in Figure (6.9) has a natural frequency of 64 Hz and 
must pass a 10G sinusoidal vibration test and a 50G shock with a 
half sine pulse of 0.005 second. Select a set of isolators that will 
permit these constraints. 


i 

8 0 

k 

40 Pin 

i 

50 



r 

E-3 


c\ . 

L=2.0 

0.005 

n 

M ► 


B=12.0 


Figure (6.9). A PCB Subjected to Shock and Vibration 
System without Isolators 

First study the system’s response in vibration without isola¬ 
tors: 


X 


e 


9-8 G in T 
f 2 

J n 



X e - 0.0957 inches 


The maximum allowable displacement is 
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0.00022 B 


X r 


C hr -Ji 


B -12 
C = r -1 
h = 0.08 
L = 2.0 


X mSiX = 0.0233 inches 


So there is a great chance that the system will fail under the vibra¬ 
tion condition. 

System with Isolators 

By using a trial and error method, we could assume an isola¬ 
tor frequency and then evaluate the system first for shock response 
and second for vibration environment. 

Design 1 : Assume a set of isolators with a resonant frequency 
of 5 Hz to provide vibration isolation for the 64 Hz PCB. For shock 
environment the input factors are: 

Gin = 50, 

/„ =5, and the shock pulse frequency is /„ =---= 100 

" p 2x0.005 


K _ isolator frequency _ 5 _ ^ 
pulse frequency 100 

A = 2R . From here A =0.10 


(9.8)(50X0.1) 


X = 1.96 inches 


This dynamic displacement value is too large for this applica¬ 
tion. A more appropriate value should be in the range of 0.3 or 0.4 
inches. 


Design 2 : Assume a set of isolators with a resonant frequency 
of 20 Hz. 


isolator frequency = . From here A- 0.4 

pulse frequency 100 
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=> X = 0.49 inches 


Y _ (9-8)(50)(0.4) 


This value may still be too large. 

Design 3 : Assume a set of isolators with a resonant frequency 
of 30 Hz. 


R = — = 0.3, 
100 


X = 


(9.8)(50)(0.6) 

30 2 


A = 2R . From here A= 0.6 

X = 0.33 inches 


This value is acceptable. We need to check for vibration. For 
small levels of damping, assume £,=0. Then the dynamic Magnifica¬ 
tion factor is 


where r = 

1 - r 2 /natural 

Notice here that D m is the dynamic magnification factor of the iso¬ 
lators and not of the PCB. 

r = 64/30 = 2.1333 => D m = - 0.2816 

A negative sign only means that the system response is out of 
phase with the input). 

Xj = X e D m => Xj = (0.0957)(0.2816) => X ,■ = 0.0269 inches 

This is larger than the allowable displacement and therefore 
this set of isolators are not acceptable. 

Design 4 : Assume a set of isolators with a resonant frequency 
of 25 Hz. 

R=-^- = 0.25 A = 2R . From here A = 0.5 

100 
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x = ( 9.8)(5 0 )(0.5 ) x = 0 392 i nc hes 

25 2 

Now let us check for vibration. 

r = 64/25 = 2.56 D m = - 0.18 


X,- = X e D m => X t = (0.0957)(0.18) => X t = 0.0172 inches 

This value is below the allowable limit and therefore the de¬ 
sign is acceptable. 

A Velocity Shock Isolation Example 

Consider a 20 pound electronics box which has a fragility level 
of 12 G’s. In other words, it can not withstand more than 12 G’s. 
We like to find the maximum height it may be freely dropped with¬ 
out any damage to the system. The natural frequency of the system 
is 125 Hz. The conditions are such that there is no rebound. 

There are two points for consideration here. First is that the 
system may be assumed to behave similar to a spring-mass system 
at the moment of impact. Thus by equating the kinetic energy of 
the system with the energy required to deflect the spring a distance 
A , we obtain the maximum deflection. 


= v Ie 1 = _ 

V k Inf/ 


(6.13) 


where f n is the natural frequency of the system (without isolators) 
or of the isolators if they exist. V is the change in velocity and is 
calculated by 

V^cjlgh. (6.14) 

g is gravitational constant, h is the drop height and C is the coeffi¬ 
cient of rebound. It is equal to 2 for a full rebound and 1 for no re¬ 
bound at all. 

It may be shown that the transmitted shock acceleration is 
given by the following relationship. 
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M n v 

g 


(6.15) 


By combining equations 6.14 and 6.15 and rearranging the 
terms, an expression for the height may be developed: 


h = 


8 n 2 C 2 f n 


Now, we may substitute the values relevant to this problem: 


h = 


(12 2 )(386) 


(8)(3.14 2 )(1 2 )(125) 


h - 5.64 inches. 

Therefore, if the electronic box is dropped from any heights above 
5.64 in., it will be damaged. Suppose that this system had to sur¬ 
vive a drop form a height of 12 inches. What type of isolators 
should we choose? To do this, we need to calculate the change in 
velocity: 


V = V2(386)(12) = 96.25 


in 

Sec 


From equation 6.15, calculate the natural frequency of the isola¬ 
tors. 


^mav g 


max 
J n 

f n ~ 


2xV 

(12)(386) 


2tt(96.25) 
f n = 7.66 Hz 

The dynamic deflection of the isolator is: 
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96.25 
2^(7.66) 


1.999 in. 


The displacement of this box will be approximately 2 inches. To 
calculate the stiffness, 





K=4n 2 f Z M 
K = 4(3.14 2 )(7.66 2 )f^ 


in 

This is the stiffness of the dynamic system to ensure survival form 
a 12 in. drop test. 

Maximum Desired PCB Deflection 

Similar to vibration, Steinberg (15) suggests that if the peak 
single-amplitude displacement of the PCB is limited to X mdiX , the 
component can achieve a fatigue life of about 10 million stress re¬ 
versals in a sinusoidal vibration environment. For a shock condi¬ 
tion, 


* - ti.mniB inches (6. 16 ) 

max ChrJZ 

The desired corresponding frequency is given by /^ 
f 9.8<7 ■ AChr-JT 

r _ _ Irl _ 

Jd ~ 0.001325 

V 

Where A is shock amplification factor (Gout/Gin) typically is 



(6.17) 
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Isolator Frequency 

A — L A , A = —--—-—- . 

Pulse Frequency 

C,h,r and L are defined by equation 6.11. 

Equipment Design 

Aside from the PCB, the electromechanical engineer must be 
mindful of the design aspects of other parts of the packaging. 
These include the chassis design, which is generally formed from 
sheet metal or cast from a light metal and machined; the cabinet 
design and its structural strength as well as its shock and vibra¬ 
tion issues; and, the vibration of wires and cables as well as their 
interaction with other components that may lead to failure. 

Methods of construction of the chassis or the cabinet play an 
important role in the survival of the system as a whole in vibration 
and shock environments. For example, the friction between bolted 
or riveted joints dissipates energy and as a result represents a 
great level of damping in the system. Therefore, these types of 
structures are superior to welded structure in vibration environ¬ 
ments. Furthermore, in welded joints, stress concentration is a 
common defect leading to low fatigue life. Reducing the number of 
welded joints as well as their heat treatment to reduce the residual 
stresses may eliminate this problem. 

As for wires and cables, it is prudent to tie wires that extend 
in the same direction, support the harness length, minimize the 
lead length, and clamp the wire near the termination to a struc¬ 
ture. 

Induced Stresses 

In calculating the natural frequencies and mode shapes, the 
applied loads (either forces or acceleration, e.g., gravity “G” loads) 
do not play a role. However, when we speak of a system’s failure in 
shock or vibration, we are making a reference to the deflection (and 
by necessity to stress) at a critical point in the system. These val¬ 
ues may not be evaluated without knowledge of the loads applied 
on system. 

Fortunately, calculating the system’s response under dynamic 
loads is not difficult at all. The principles of Strength of Material 
still holds as one realizes that the components used in an elec¬ 
tronic enclosure are beams, plates and frames. The question re¬ 
mains; how should one treat the “dynamic” loads? The steps to 
calculate the induced stresses are as follows: 
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1. Calculate the natural frequencies. 

2. For a given load, calculate transmissibility of the system. 

3. Calculate the dynamic load (a product of transmissibility, 
stress concentration factor and static loads). 

4. Calculate the stresses from the dynamic loads. 

While an engineer must develop a solid foundation of the 
analysis tools that he/she uses, the techniques described here 
would be best utilized for developing a “sense” for the system. Now 
a days, finite element analysis (FEA) tools can be readily used to 
calculate the loads and stresses on the electronics enclosures. 

Sample Problem 

As an example, let us consider the beam problem that we 
used previously to calculate the dynamic magnification factors. 
Now we need to understand the stress levels so an appropriate 
material may be selected. 


10 Sin( TO t) 


- > 

12.00 

Figure 6.10. A Cantilever Beam Under a Vibratory Load 

Based on Strength of Materials theories, the stresses in a 
beam are given by the following relationship. 

_ Me 

°"max _ ~j~ ■ 

Where M is the bending moment, I is moment of inertia, and c is 
the distance from the outer fibers to the neutral axis. For more 
details, the reader is encouraged to read any textbook on Strength 
of Materials. The bending moment is evaluated from the Shear & 
Moment Diagram, which may be developed by using a free body 
diagram. Note that the dynamic bending moment is the static value 
multiplied by transmissibility. 
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Recall that for this system, <y=41, £=.01 and c=0.5. From a 
free body diagram (not shown here), we obtain 

M = 120 lb-in. 

Therefore a static stress value is 


_ Me 
^ static - j 


120x0.5 

0.03 


= 2000 psi 


Previously, we calculated for n = 0.73, T m =2.13; therefore, 


M= 120 x 2.13 = 256.8 lb-in leading to 


CT 1 = 


2 56.8x0,5 
0.03 


= 4280 psi. 


Similarly, 


r 

T 

1 m 

M 

a 

40/41 = 0.98 

22.6 

2715.9 

45,264 

50/41 = 1.22 

2.04 

245.4 

4090 

60/41 = 1.46 

0.88 

106.0 

1767 


Note that moment is in lb-in and stresses are in psi. Further¬ 
more, static bending stress is 2000 psi, while dynamic stress may 
rise to as much as 45,264 psi. Clearly, at dynamic loads, failure 
may occur! 

Vibration Failure Modes 

Now that we have mentioned “failure,” let us review a few 
modes of failure. In a shock or vibration environment, a system’s 
failure may depend on several factors: 

• Deflections causing chatter - This includes the failures of 
components, which have moving parts such as relays. 

• The stresses beyond yield. 

• The stresses beyond fatigue levels. 

• Sliding of cables, and connectors. 
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Therefore, a system’s failure should not be only defined as 
when it “falls apart” but when the design criterion are not met any 
longer. 

Random Vibration 

Earlier, we mentioned that in a periodic motion, displace¬ 
ments follow a predictable pattern. However, in certain vibration 
cases, the amplitude of displacements or accelerations does not 
follow a particular pattern. Should we plot these variations a func¬ 
tion of time, we would develop a graph very similar to Figure 6.11. 
This figure may be considered to be a typical time history curve for 
random vibration. 



Figure 6.11 . A Typical Time History in Random Vibration 

Application of Fourier Series enables us to change the 
mathematical form of any curve into a series of sinusoidal curves 
each with its own frequency and amplitude. This concept is shown 
graphically in Figure 6.12. 



An interesting lesson of this transformation is that all the fre¬ 
quencies within a given bandwidth are present at all times. There- 
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fore, in a random vibration environment, all resonance frequencies 
are also present. As a result, failures may be produced that can not 
be duplicated in a harmonic vibration sweep. 

In a way, one may say that modeling random vibration is an 
attempt to model nature. In the discussion on Free and Forced Vi¬ 
bration, it was pointed out that these types of vibrations are peri¬ 
odic. However, Random Vibration is different from these two in that 
it is non-periodic. While the solution to periodic vibration is deter¬ 
ministic, probabilistic approaches must be used to find the dis¬ 
placements and/or stresses in a random vibration excitation. This 
constitutes the most important difference between sinusoidal and 
random vibrations. In other words, we can guess what maximum 
accelerations, forces or displacements the system will undergo but 
we can never say with certainty at which time. 

Superficially, the same failure modes exist; however, in ran¬ 
dom vibration, all of the natural frequencies are excited at the 
same time. This leads to an overall amplification of displacements 
and/or stresses. Failures are therefore more severe. 

Solution Techniques 

The steps to solve random vibration problems are relatively 
straightforward. These are first, knowing the excitation spectrum; 
in other words, the input root mean square of disturbance (e.g. 
power spectral density (PSD)). Second, determining the equipment 
response, e.g., calculating the response power spectral density 
( G out ) and the response stresses values. Finally, calculating the 
probability of failure (e.g., IS, 2S, and 3S probabilities). 

Excitation Spectrum 

As it was pointed out in random vibration, the exact location, 
velocity and/or acceleration may not be known. However, a root 
mean square value may be determined: 

2 

S.D.-Limit (6.18) 
Af ->0 Af 

where q> may be the root mean square of acceleration, velocity or 
displacement Spectral density may be classified as 

1. Power Spectral Density (PSD) 

2. Velocity Spectral Density (VSD) 

3. Displacement Spectral Density (DSD) 
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Power Spectral Density is commonly used and is available 
from various handbooks and /or standards. 

/-*2 

PSD = Limit— (6-19) 

Af->0 A J 

where G is the root mean square of acceleration. 

Notice that since PSD and G (acceleration in gravity units) are 
related, a root mean square acceleration may be calculated from 
the PSD curves. More specifically, G rms is the square root of the 
area under the PSD curve. Note that since G rms is the area under 
the PSD curve, various PSD curves may indeed produce the same 
Grms ‘ 

Let us explore the Power Spectral Density Curve in some more 
detail. It is plotted on a log-log coordinate system and is a function 
of frequency. In other words, the X axis shows the frequency and 
the Y axis shows the power spectral density (PSD) which has the 
units of G 2 /Hz. The simplest PSD curve is a straight horizontal line 
known as white-n oise as show n in Figure 6.13. A more typical PSD 
curve is shown in |Figure 6.14 


Log 
Power 
Sp ectral 
Density 


0.20 


White Noise 
Curve 


10 


2000 


Log Frequency (Hz) 


Figure 6.13. A White Noise PSD Curve 

The square root of the area under these curves provides a root 
mea square of the acceleration loads (G rms ). For example, to cal¬ 
culate the G rms for the white noise input PSD as shown in figure 
6.13, the following steps need to be taken: 

G rms = 4PSD (f 2 -/,)= V(0.20)(2000-10) 
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G rms =19.95 

However, one should bear in mind that this is a log-log scale 
and not be fooled by “straight lines” and usage of the triangle rule - 
except if the lines are horizontal. 



Figure 6.14. A Typical PSD Curve 

To calculate the area under the sloping lines, special attention 
must be given to the basic area relations using Log equations. For 
example, using the point slope equation from analytic geometry: 

Ln PSD = mLnf + Ln n , 

m is the slope of the line and Ln n is the intercept on the PSD axis. 
This equation may be rewritten as 

PSD - mf n 

The area under the curve can then be determined by integration. 

2 2 2 

A rea = J ydx = J PSDdf = J mf n df (6.20) 

1 l 1 
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Probability Distribution Functions 

Earlier, it was mentioned that random vibration is probabilis¬ 
tic. In order to predict the probable acceleration levels, it is neces¬ 
sary to use the probability distribution function. While providing 
details is beyond the scope of this work, for a Gaussian (or normal) 
distribution curve, there are three bands of interest; namely, IS, 
2S, and 3S. S stands for standard deviation. 


1. A IS probability means that there is a 68% possibility that 
the accelerations will reach the G rms levels. 

2. A 2S probability means that there is a 27% possibility that 
the accelerations will reach the 2 G rms levels. 

3. A 3S probability means that there is a 4% (possibility that 
the accelerations will reach the 3 G rnts levels. 


The rms acceleration of the white noise depicted in Ifigure 6.13| was 
shown to be 19.3. This means that there is 68% (IS) probability 
that the acceleration levels will reach 19.3, or there is 27% (2S) 
probability that the acceleration levels will reach 38.6, or there is 
4% (3S) probability that the acceleration levels will reach 57.9. 


Equipment Response 

It may be shown that the response of a single degree of free¬ 
dom system to a white noise random vibration input can be deter¬ 
mined from equation 6.12. 


G out PSDT„ 


( 6 . 21 ) 


Where: PSD = Spectral density at resonance = G 2 /Hz 

co n = Natural Circular Frequency of System 
T m = Transmissibility at resonant frequency 


It should be mentioned that the displacement of a PCB in vi¬ 
bration behaves similar to a single degree of freedom system and 
therefore this approach may be utilized. Furthermore, the above 
equation is only va lid for a PSD input that is flat close to reso¬ 
nance, as shown in Figure 6.15 However, the errors of deviating 
from this rule are quite small in most practical cases. 

A frequently asked question is how a spring mass system (a 
single degree of freedom) would vibrate when subjected to a (flat 
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spectrum) random vibration environment. The answer is that while 
the system cam only vibrate at its natural or resonant frequency, it 
is the amplitude of vibration which varies. Steinberg (15) suggests 
using the term “the number of positive zero crossings” as opposed 
to “resonant frequency.” This is the number of times the displace¬ 
ment trace crosses the zero axis with a positive slope and For a 
single degree of freedom system, this is exactly equal to the reso¬ 
nant frequency. 



Figure 6.15. Proper PSD and Transmissibility Values For Equipment Re¬ 
sponse Calculations 

In most practical cases, the magnification factor must be 
evaluated from available test data and in its absence must be esti¬ 
mated. The general rule is that low displacements and low strains 
lead to high transmissibilities and visa versa. Therefore, high fre¬ 
quencies, low input G’s or any other factor such as high pressure 
interfaces between ribs, supports, connectors and conformal coat¬ 
ings tend to lower transmissibility. Earlier, some guidelines were 
provided for periodic vibration environments. The same rules hold 
here as well. 

Dynamic Loads 

Dynamic loads (forces and moments) may now be evaluated 
from the output acceleration. One has to exercise caution in using 
stress concentration factors if there are holes or notches in the 
system and to take proper units into account. In the British Sys¬ 
tem, dynamic Force = Pounds x G oul . Whereas, in other systems, 
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G out must be converted into an appropriate units. Based on these 
values of forces and moments, dynamic stresses may be developed. 

Probability of Failure 

There are three levels of stresses corresponding to the three 
levels of the probable G out levels. As an example assume that 

IS stress is 4624 psi. Then, 2Sand 3S stresses are 9248 psi and 
13872 psi, respectively. If the material is aluminum - with a yield 
stress of 12,000 psi, then there is a 4% chance that the part will 
undergo yielding. 

Maximum Desired Deflection 

Similar to vibration, Steinberg (15) suggests a maximum de¬ 
sired PCB deflection for at least a 10 million vibration cycles: 

X - Q Q002 2j8 inches (5.22) 

ChryfL 


The corresponding resonant frequency is 

- a 0.8 


fd~ 


29 AChr 


(f)(PSD)L 


0 . 00022 B 


(6.23) 


B, C, h, r, and L are defined for equation 6.11. 

Exercise 

Calculate the Stresses for the beam shown in 


Figure 6.11 us¬ 


ing the previous white noise PSD shown in i Figure 6.13.1 Determine 
the probability of failure for 6061-T6 Aluminum and Stainless 
Steel. Assume a 10% damping in the system due to coatings, wir¬ 
ing and so on. The natural circular frequency is 41. For Aluminum 
the yield stress is about 12 ksi and for steel it is 40 ksi. 

The calculation of stresses is as follows: 


^out ^/%jPSD T m ! 2 


T =l/2£ at resonance =-= 5 

m 2 0.1 
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•v/41 • 0.25 • 5 


out 


3.58 


F = /Z>-G ou t =10-3.58 = 35.8 lbs 


M= F-g = 35.8-12 = 429.6 lbs-in 


(7 = 


KMc 

I 


429.6 0.5 
0.03 


= 7160 psi 


K is the stress concentration factor. Since there are no holes 
or notches in the system, K = 1. 

a =7160 psi is the IS stress level compared to a static value 
of 2000 psi. The 2S and 3S are 14320 (= 2 x 7160) psi and 21480 
(= 3 x 7160) psi, respectively. Clearly, there is a 27 % chance that 
aluminum will fail and thus it is not a good choice. 

Multiple Degree of Freedom Systems 

When a multiple degree of freedom system is subjected to a 
random vibration input, the motion of each mass in the system will 
influence the motion of every other mass in the system, as shown 
in Figure 6.15. 



Figure 6.15. Response of Multiple Degrees of Freedom Systems. 


Theoretically, the response of the multiple degree of freedom 
system to a random vibration input can be determined using an 
integration technique. The IS G rms response is obtained by inte- 
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gration under the area of the random input curve and transmissi- 
bility output curve, with the use of equation 6.16. 


< 624 > 

For a more detailed treatment of this subject see Steinberg (15). A 
note of caution: in this book, we have approached solving engi¬ 
neering design issues using “back of the envelop” techniques. We 
need to keep in mind that any time we analyze multi-degree of 
freedom system without the use of numerical techniques and com¬ 
puters, we run the risk of making errors. The chance of making 
errors increases as with the increase in the degrees of freedom. 

Numerical techniques, in general, and finite element analysis 
(FEA) in particular have become the tool of choice for many engi¬ 
neers and researches to analyze and simulate various configura¬ 
tions. These tools, however, must not be used blindly. To develop 
an understanding of FEA, the next chapter is given to this topic. 
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7 

Introduction to Finite Element Method 


In the previous chapters, we developed tools for developing an 
understanding of heat transfer and shock and vibration issues. 
While these tools may suffice for solving design issues, they are not 
appropriate tools for conducting an in-depth analysis. For realistic 
problems of today, a numerical method - and mainly finite element 
method - is utilized. Nowadays, it is possible to use CAD software 
to “design” a system and then by the click of a button, solve for 
stresses and strains in various components without the slightest 
notion of any degree of accuracy. 

Since a CAD model is the closest to the “real thing,” an inter¬ 
esting question arises: what is the difference between a Real-World 
system, an Engineering model, and a FEA model? 

As engineers, we transform real-world systems into engineer¬ 
ing models by making “engineering” assumptions. This enables us 
to simplify the physics of the system to a level where realistic solu¬ 
tions may be sought. By further transforming these models into 
idealized mathematical models, FEA models are developed and are 
characterized by nodes, elements, and boundary conditions. 

The purpose of this chapter is to provide a basic understand¬ 
ing of the Finite Element Methods (FEM) and its proper applica¬ 
tions in the context of stress analysis, vibration and heat transfer. 

The basic idea in finite element method is to find the solution 
of a complicated problem by replacing it with a simpler one. The 
earliest application of FEM may be considered to be the ancient 
technique of finding the lower and upper bounds to the circumfer¬ 
ence of a circle. However, the FEM as known today was presented 
by Turner and his colleagues in 1956 for the analysis of aircraft 
structures. 
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An Example of Finite Element Formulation 

The following is an illustration of finite element formulation 
for a bar with two different cross sectional areas clamped at one 
end, and loaded with a lateral force at the free end as shown in 
Figure 7.1. Our goal is to calculate the displacements of this bar at 
various locations particularly where the two cross sections meet 
and at the end where the load is applied. Note that in this formula¬ 
tion displacements are calculated first. They are called the primary 
variables. If stress or strain field is needed, it must be developed 
from the calculated displacements. 



Figure 7.1. Two Bars Loaded at One End 


Consider a generic line element of length 1 (as in Figure 7.2) 
and assume that the displacement varies linearly from one end to 
the other - each end is called a node. Then, a simple finite element 
model of this system is as in Figure 7.2. 



Figure 7.2. A Simple Line Element of Length 1 
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Figure 7.3. The Finite Element Model of the Original System 

Now, we assumed that the displacement varies linearly from 
node to the other. In other words: 

u(x) - a + bx (7.1) 

Let us apply the boundary condition: 

atx- 0 m(0) = Uj 

at x = 1 u(l) - uj 

Thus one concludes that 
a - Uj 

Uj - Uj > 

By substituting these relationships into equation (7.1), we de¬ 
velop the elemental equation for displacement based on nodal val¬ 
ues: 

u(x) = ^1 - ~~j u i + 0 < x < l (7.2) 

In general this equation may be written as 
n 

U(X)= X (PjUj 

i = 1 

where (p l ’s are called the shape or trial functions. In this case, 
n =2, and 
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The goal of this formulation is to find the displacements. 
Thus, from a stress-strain relationship, we can drive a force- 
displacement relationship as follows: 


o = Ef , 


where o is stress, E is young’s modulus of elasticity and e is 
strain. Multiply each side by the cross sectional area, and we ob¬ 
tain 


A o = EAf 

Note that the right hand side is equivalent to Force and strain 
is the derivative of displacement with respect to length; Thus, 

r r* du 

F = EA— . 
dx 

For the each element, we balance the force at each node. For the 
first element: 


(f), -<f), 
-(f), <f), 


(7.3) 


For the second element: 


(fh ~(fh 

~(f ) 2 (f ) 2 


-f 2 

F3 


(7.4) 
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Notice that in a general formulation, the stiffness and the 
length of the two elements are different. Now combine the two ele¬ 
ments 



(7.5) 


In matrix notation, this equation is written as 


[K]{u} = {F} 


In general, if [K] - defined as the characteristic or stiffness 
matrix - depends on {11} the problems is nonlinear and if [K] does 
not depend on { 11 } the problem is said to be linear. The size of the 
matrix [K] for practical problems easily exceeds several million en¬ 
tries. Thus, a great deal of research has been dedicated and con¬ 
tinues to be dedicated to solving these equations most efficiently. 

To obtain the displacements, first, we need to apply the 
boundary conditions to equation 7.5 (i.e. zero displacement at node 
1 and force equal to F at node 3) and specify the values of element 
length as well as stiffness. Then, we can solve for « 2 and «3 • 

Formulation of Characteristic Matrix and Load Vector 

In the previous example, we developed the finite element ma¬ 
trices based on our understanding of the physics of the problem. In 
general, there are three approaches to this formulation: 

1. Direct Approach 

2. Variational (Energy) Approach 

3. Weighted Residual Approach 

Direct Approach 

In this approach, as in the previous example, the basic phys¬ 
ics is utilized in order to form the elementa l matrices. As an exam¬ 
ple consider the bar problem in Figure 7,4, 
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Figure 7.4. Bar with Different Cross-Sectional Areas 

From physics, we know that force and displacement have the 
following relationship. 

[K]{u} = {F} 


or 



To calculate the [K] matrix, we need to take advantage of in¬ 
fluence function or coefficient. The stiffness influence coefficient is 
defined as the force needed at node i (in the direction of ujj to pro¬ 
duce a unit displacement at node j (i.e., uj= 1) while all other nodes 
are restrained. Therefore, 

Ku=Fi. 

However from strength of material 
F = a A , and 

K h = F i = EeA = EA — 
u i ( 

where t is element length and At is the change in length. Recall 
that 
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Thus 


At - u ■ - u ■ =1-0 = 1. 
J 1 


K ii 


EA 

1. 


Note that for calculating Ky a negative force must be used because 
of the direction of the coordinate system: 



The matrix equation may now be assembled easily. 


E A 

( ~ h 

,EA 


.EA, EA, 


(“)l ( 7~)l + (~T~)2 


£ 

0 


-03 


0>3 


"1 

« 2 

“3 


f 2 -f 2 

f,-f 3 

F 4 


Variational (Energy) Approach 

In this approach, one takes advantage of the fact that nature 
would take the path of least resistance, i.e., minimum energy ex¬ 
penditure to change from one state into another. Thus by express¬ 
ing the energies involved and then by minimizing it, not only the 
equations of "motion” are defined but also the boundary conditions 
are determined. The advantages of variational formulation is as 
follows: 
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1. Generally the energy functional (7 ) has a clear physical 
meaning. 

2. The functional contains lower order derivatives of the field 
variables. 

3. Some times both upper and lower bounds may be found. 

4. Complicated boundary conditions - either natural or free 
are satisfied. Only geometric - or forced - boundary condi¬ 
tions need to be imposed. 

For our example: 

I = Strain Energy - work done by external forces 


/ = in*-w 

k H 




Af — a e dr 
0 2 


and a = Ee , so fl. = 




EA; 2 


i 

V 

o 


dr 


but 


£ - 


u j~ u i 


t 


substitute back in the Strain Energy relationship and obtain 
E A ( 9 9 

El /, = - U~ +U; - 2U;U J 

k 21 \ 1 J 1 J 


or in matrix form 


EA 
2 1 


[ 

' 1 -f 


U; U • 



1 J 

-1 1 

hi 


The work done by external forces is the sum of applied forces 
in the direction of displacements: 

W = YF-u- 

p ^ ii 

i 


The functional may now be formulated as follows. 
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-1 1 


'J)k 


~(Z f i u ih 


This energy functional must be minimized with respect to dis¬ 
placements; thus, 


dl 

d “i 


= 0 


for i = 1,2, and 3 


This differentiation leads to the following set of elemental equations 
which may be assembled for the overall equation. 


EA 

' 1 -f 


I 

[<->* 

-1 1 

k 

t ‘j 


Weighted Residual Approach 

Another approach to develop the elemental equations is the 
Weighted Residual technique. The first step of this approach is 
based on the idea that in many practical problems, the governing 
differential may be obtained much easier than the energy func¬ 
tional. The next step is to find an approximate solution - such as 
finite element approximation - and minimize the error. This step is 
expressed mathematically as 


Ku = F or Ku-F = 0 
If u is an approximation to u then 


K u- F = R 


Needless to say, the best approximation provides the lowest 
error, i.e., lowest R. To achieve this, the method of Weighted Re¬ 
sidual Approach dictates that we multiply the error by a "weight" 
function, integrate the product over the domain of interest, and set 
the outcome equal to zero. 

This technique is best suited for areas of physics such as 
fluid dynamics, where development of the differential equations is 
much easier than the energy functional. 
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Finite Element Formulation of Vibration Equations 

Recall that the finite element matrix equation has the follow¬ 
ing general form: 

[K]{u} - {F} 

This is has the same form as the force deflection relationship 
in a spring mass system. One may say that to develop the equa¬ 
tions for a finite element analysis, a continuous (engineering) sys¬ 
tem is broken into a series of interconnected spring-mass-damper 
systems. 

It is important to consider that any time a continuous system 
is broken into a spring-mass-damper system, some level of error is 
introduced. The designer/analyst must have a handle on this er¬ 
ror. As an example, assume that a 10% error in natural frequency 
calculation has taken place. This error leads to a dynamic magnifi¬ 
cation factor of 5 instead of 50. Clearly, large errors in displace¬ 
ment may be introduced if the natural frequencies are not calcu¬ 
lated accurately. 

There are two classes of dynamic response problems: 

1. Wave propagation type where the variables change with 
time but no periodic behavior is identified. 

2. Vibration types where a periodic/harmonic response is ex¬ 
pected. 

Wave Propagation Type 

Let us look at the vibration and numerical solution of this 
problem. The equation of motion for a spring-mass system is: 


Mu, tt +Cu, t +Ku = f(t) (7.6) 

In this equation,/(t) is a general time-dependent forcing 
function. Now, time derivatives may be replaced by an approxima¬ 
tion similar to the following. 

u, t = Ut+At+U ‘ (7.7) 

At 

u - u ‘ +2 + u t (7 3 ) 

Ht At 2 
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Equations (7.7) and (7.8) point out that to calculate the veloc¬ 
ity of at time t , the location at time t + At must be known. Fur¬ 
thermore, to find the acceleration at the same time, the location at 
time t + 2At must also be known. By substituting this approxima¬ 
tion in equation (7.6), we obtain: 


M 

9 u t+2At 
At 2 


[At 


2 M 
At 2 


\ u t+At + 


M 

A^ 


C_ 

At 


+ K 


u t = f(t ) 


At time t = 0, the initial conditions, i.e., the location of the 
mass as well as it velocity must be known. Then, by taking advan¬ 
tage of equation (7.7), the location of the mass at time At may be 
calculated. The calculation may begin in earnest at time 2 At by 


u 2At - 


A r 


f.. .. 1 

' C 2M ^ 


( M C J 


S' 

II 


|»A t ~ 

—--+ K 



At 2 j 


,A t 2 A/ ; 

/ 


Once x 2 ^ t is evaluates, x 3A , is evaluated as follows. 


h 3A t 



/(AO¬ 


's. 


r C_ 
k A t 


2 M 
At 1 j 


U 2A t 




\ 


This approach is called time marching and may be continued 
to evaluate the displacement at any time. However, there is an in¬ 
herent inaccuracy in this approximation and as time progresses, 
the errors accumulate and eventually the approximate solution be¬ 
comes completely erroneous. The task of finding proper time ap¬ 
proximations to alleviate this problem rests on the shoulders of 
mathematicians and some sophisticated strategies have been de¬ 
veloped. A discussion of these strategies is beyond the scope of this 
work. 

Despite the inherent inaccuracies, this approach is applicable 
to any dynamic system provided that the mass, damping, and 
stiffness matrices are available. Furthermore, f(t) may be any 
function of time including periodic functions. However, this is not a 
good approach to vibration problems because 

1. The time step required to model several periods accurately 
will be prohibitively small for any realistic problem 
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2. Natural frequencies may not be readily calculated. The time 
response is a combination of all frequencies involved and 
discrete values may not be extracted. 

Vibration Type 

Another mathematical approach for solving time dependent 
problems is to utilize Fourier Transform and construct a time re¬ 
sponse curve using frequencies and mode shapes. To do that, the 
first step is to calculate the natural frequencies and mode shapes. 
Earlier, it was shown that damping has little effect on the value of 
natural frequencies, therefore, only the undamped equations will 
be considered. 

+ [K]{U} = 0, 

Assume displacements as follows: 

{U} = {X} e " ia>t 
Thus, we obtain 

[ [K] - to 2 [M] ] {x} = 0 

This means that the determinant of [[K] - co 2 [M]] must be zero. 
Therefore values of co must be found to satisfy this condition. 

Methods of Solving This Equation 
Two general types of methods are available: 

1. Transformation methods such as Jacobi, Givens and 
Householder; these methods should be considered when all 
of the frequencies are needed. 

2. Iterative methods such as power methods; these methods 
must be considered when few frequencies and mode shapes 
are required. 

It should be pointed out that there are two possible ways of 
formulating the mass matrix M. The first is called Consistent Mass 
Matrix. This matrix is fully populated. It is also possible to di¬ 
agonalize this matrix, (i.e., only the terms on the diagonals are 
nonzero). In effect, it is assumed that the mass of the systems is 
lumped at the nodes. There are advantages and disadvantages in 
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using either of the methods. This will be illustrated in the following 
example. 

Example: Find the natural frequencies of longitudinal vibra¬ 
tion of the unconstrained stepped bar shown below: 



Figure 7.5. Longitudinal Vibration; The First Rod has a cross-sectional 
area twice the second one. 


Suppose that we only use two elements: 


[*](D =<^>1 




" 1 

-f 

_ 4EA 

' 1 

-f 

-1 

1 

L 

-1 

1 / 

' 1 

-1' 

2EA 

' 1 

-f 

: -1 

1 

L 

-1 

1 


Consistent Mass Matrix - To use this formulation, one must 
use the energy approach to develop the mass matrix. The result is 
as follows. 


[^]<b=(^)i 

[M\ 2 )=(^t)2 


2 r 

pAL 

2 r 

1 2 

6 

1 2J 

"2 f 

pAL 

'2 r 

1 2 

r 12 

1 2 


where p is density. The assembled matrices are: 
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[M] = 


pAL 

12 


2 -2 
-2 3 

0 -1 
[4 2 0 
2 6 1 
0 1 2 


0 

-] 

1 


The determinant of [[K] - (O 2 [M] ] must be zero. 

pL 2 6) 2 

Define p 2 - —- and form the determinant: 

24 E 


" 2EA "j 

V L J 


2(1-2 J3 2 ) 
-2(1 + J3 2 ) 
0 


“ 2(1 + p 2 ) 
3(1-2 p 2 ) 
~{1 + P 2 ) 


0 


-(l + y9 2 ) 
O-2/0 2 ) 


= 0 . 


This equation may be expanded to provide the following equation. 
18/? 2 (1 - 2/? 2 )(/? 2 - 2 ) = 0 . 


Notice that there are only three degrees of freedom in this 
system, i.e., ui, U 2 , and U 3 . Therefore, there are only three possible 
solutions to this equation: 

p 2 -0 => a 2 = 0 => co\ - 0 


A zero frequency corresponds to a rigid body motion; i.e., the 
entire system moves together. This is caused because no boundary 
conditions have been applied. The other two modes of vibrations 
have the following frequencies. 


P 2 


2 


P 2 = 2 


2 

_ 12E 

, ^ E 

=> ^2 : 

'y 

=> = 3 - 46 a 2 


pL 2 

] P L 2 

2 

_ 48E 

1 E 

=> 0)3 = 


=> a? 3 = 6.92 


PL 2 

Vp l 
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Lumped Mass Matrix - To use this formulation, one must 
“lump” the mass at each node. For this example we obtain: 


mm =(-^)i 




'1 o' 

_ pAL 

'1 o' 

0 1 

2 

0 1 

" 



'2 r 

_ ML 

1 0 

1 2 

4 

0 1 


This leads to the assembled mass matrix as following. 
pAL 


[A/] = 


2 0 0 
0 3 0 
0 0 1 


Define / 3 2 = pL 02 
8E 


and form the determinant: 


m 


2(i-fin 


-2 

3(1 ~fi 2 ) 


•1 


1 (1-J3 2 ) 


0. 


This equation may be expanded to provide the following equation. 

6/3 2 {l~ fi 2 )(fi 2 - 2) = 0. 

This leads to the following frequencies: 

fi 2 = 0 => a ) 2 - 0 fti] = 0 a rigid body motion 


>2 __ 2 8 E 


fi = 1 COj - 


fi 2 =2 


pL 2 
2 16E 

"3 =-f 

pL 2 


a> 2 = 2.83 

r 

"3= 4 , 


pL 2 


pL 


Impact of Boundary Conditions 

Now include the effect of boundary conditions. To do so, we 
need to go back to the assembled equation [ [K] - to 2 [M]l and set 
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the rows and columns associated with the zero displacements 
equal to zero. This in effect reduces the size of the matrix by the 
number of degrees of freedom set to zero. For this example, the 
first row and column must be set to zero and the resulting equa¬ 
tions for the stepped bar are as follows. 

Consistent formulation 

3(l-2/? 2 ) -(1 -/? 2 ) =q 2 PLV 

-(l-/? 2 ) (1-2 p 2 ) ’ 24E 

This equation leads to 
<y, = 1.96 
(o 2 = 5.16 

Lumped formulation 

-1 (1 -ft 1 ) 8E 

This equation leads to 

co x =1.84 

co 3 = 3.55 

Comparison with Exact Solution 

In this example, we assumed that the bar has two different 
cross-sectional areas. How would the results change if the two 
links were the same size (as the smaller one)? Carrying the equa¬ 
tions leads to the following frequencies for the consistent mass 
matrix. 
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p 1 = 0 => co] = 0 


a>l = 0 


p 2 co f = -~7 => &>2 =3.46 


p\f 


p 2 =2 


2 16E 

®3 = — 




&>3 = 6.92 






for the lumped matrix 


P = 0 => co 2 = 0 => <y| = 0 

&) 2 = 2.83 


p 2 =\ 


co] 


8E 


pV 


pL l 


P 2 =2 


2 16E 

"3 =—y 
pL 2 


co 3 - 4 




These results indicate that as long as the boundary conditions 
have not been applied, (in other words, the bar has free-free 
boundaries), the cross-sectional has no impact on the frequencies. 

Once the boundaiy conditions have been applied, the results 
may be compare with the exact solution. Table 7.1. shows the re¬ 
sults. 


Table 7.1. Comparison of Results 


Solution Type 

(01 

(02 

Consistent 

L6t JJ 

5 - 63 U 

Lumped 

Pi 

3 55 iJ 

Exact 

Pi 

4 - 71 JJ 


It was clearly shown - by way of the example - that the lower 
the frequency to be calculated, the more accurate it is expected to 
be. Theoretically, one may calculate as many frequencies and mode 
shapes as there are nodes in the system. However, from a practical 
point of view, the accuracy of calculations drops very rapidly - as 
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shown in the above example. To increase accuracy, three general 
rules may be cited: 

1. Apply the boundary conditions correctly. Application of 
boundary conditions is not always straightforward. An ex¬ 
ample of ambiguity in this application is presence of fric¬ 
tion. 

2. Use much more nodes in the system than the required 
number of frequencies. Generally speaking, there should be 
enough nodes to describe the shape of the mode shape ac¬ 
curately. 

3. The lumped mass matrix gives very a good solution if only 
the first mode shape and frequency needs to be calculated. 
The advantage of consistent formulation is that stress and 
deformations are calculated with more precision. 



a) Good Approximations 



b) Poor Approximation 

Figure 7.6. Approximations for Different Mode Shapes 


Finite Element Formulation of Heat Conduction 


For the sake of completeness, let us review the finite element 
formulation of heat conduction in a bar with two different cross 
sectional areas clamp ed at one end and at a temperature of 7’ ] 


as 


shown in i Figure 7,7.1 Also, heat is being generated at each of the 


bar as well as where the two cross sections meet. Our goal is to 
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calculate the temperature distribution in this bar particularly 
where the two cross sections meet and at the free end. In this par¬ 
ticular situation, we assume that the problem is steady state and 
that not other mode of heat transfer exists. 



Figure 7.7. Heat Conduction In A Bar With Different Cross Sectional Areas 

Consider a generic line element of length L and assume that 
temperature varies linearly from one end to the other - each end is 
called a node. Then, a simple finite element model of this system is 
as in figure 7.8. 



Figure 7.8. A Finite Element Model of the Bar 


For a one dimensional heat conduction problem, we have 

0~(7 hot -r cold ). 

Thus, we utilize this equation and balance the heat equation (i.e. 
heat in less heat out equal to zero) for each node as follows. To 
maintain consistency, assume that the temperature at each node i 
is lower than its surrounding. 

At node 1: 
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{T 2 -T x ) + Q x = 0 


(7.9) 


K\*\ 

k 

For node 2: 

!h^L { T 3 -T 2 ) + ^(T ] -T 2 ) + Q 2 =0 (7.10) 

L 2 M 

At node 3: 


■(T 2 -T 2 ) + Q 3 =0 


Now, we may assemble the matrix equation: 


( X )] " ( T )l 0 [T 

-(f)! -(f ) 2 ^ 

0 " (J T )2 ( T >2 7,3 


(7.11) 


(7.12) 


One may easily recognize the similarity of this equation with that of 
equation 7.5. 

Some Basic Definitions 

So far there has been mention of terms such as nodes and 
elements but they have not been defined. 

Node: A node is a location in the model where variables, 

such as displacements, temperatures, etc. are calculated. A node 
contains degrees of freedom. 

Element: A building block for the model, it contains and dic¬ 
tates the relationship between nodes. 

Element Connectivity: It is a fist of nodes that make up an 
element. Furthermore, it is element connectivity that allows infor¬ 
mation to be shared in between elements 

Higher Order Elements: Recall that the finite element ap¬ 
proximations makes use of shape functions (also called trial func¬ 
tions or interpolation polynomial). Should this shape function be of 
the order of two or more, the element is known as a "higher order 
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element." In higher order elements, some interior nodes are intro¬ 
duced in addition to the corner nodes in order to match the num¬ 
ber of nodal degrees of freedom with the number of generalized co¬ 
ordinates in the shape function. 

Isoparametric elements: One of the strong points of the Finite 
Element Methods is in its ability to model curved geometries - such 
as round boundaries. If the same shape function is used to define 
both the field variables as well as the geometry, then the formula¬ 
tion is said to be isoparametric. 

Subparametric elements: if the interpolations function used to 
define the field variables has a lower order than that of the geome¬ 
try, then the formulation is said to be subparametric. 

Superparametric elements: if the interpolations function used 
to define the field variables has a higher order than that of the ge¬ 
ometry, then the formulation is said to be superparametric. 

The Finite Element Analysis Procedure 

To take advantage of any FEA approach, one must develop a 
basic understanding of the physical and engineering issues; a ba¬ 
sic understanding of the fundamental concepts of finite element 
method; and, a knowledge of the capabilities and limitations of the 
approach used. 

The six steps of finite element analysis 

Step 1. From the physical to the FEA representation: In this 
step, one has to decide what the true concerns are; i.e., what is it 
that needs to be determined; overall deflections, localized stresses, 
etc. By knowing what variables are needed, appropriate elements 
can be chosen which lead to an overall knowledge of how to make 
the transition from the physical to the FEA model. 

Step 2. Discretization: Once the FEA model is known, it is 
then subdivided into a number of elements. The analyst must take 
into account changes in geometry, material properties, loading, etc. 
Furthermore, questions on the size and the number of elements 
used as well as simplifications afforded by the physical configura¬ 
tion of the body and loading must be addressed. 

Step 3. Application of loads: In general, this step is simple 
unless loads depend on geometry or other variables. 

Step 4. Application of Boundary Conditions: constraints must 
be applied where physical boundaries are constrained or when 
modeling simplifications - such as symmetry - are used. 
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Step 5. Assembly of element equations and the solution 
phase. In today's sophisticated FEA programs, this step is done 
automatically. 

Step 6. Review of the solution and validation: This is the time 
that the analyst would study the field variables. Care must be ex¬ 
ercised here not to treat the results as absolute. The results must 
be validated in order to gain confidence in the solution. 

Steps 1 through 4 are generally called the preprocessing 
phase, step 5 is the solution phase, and step 6 is the post¬ 
processing phase. 
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8 

Design and Analysis for Mechanically 
Reliable Systems 


Reliability concerns those characteristics which are dependent 
of time - Stability, Failure, Mean Time To Failure/Repair, etc. It 
may also defined as the probability of a product or device per¬ 
forming adequately for the period of time intended under the oper¬ 
ating conditions encountered. 

However, if harm or injury occurs due to poor reliability, the 
manufacturer or distributor may have a responsibility to compen¬ 
sate for these losses and/or injuries. The general liability law holds 
that all those who introduce a defective product into the market 
are generally liable for the product and harm or injury caused. 

There is another definition: that which is to one’s disadvan¬ 
tage. This is not a legal definition but in many circumstances, it 
may mean loss of the customer’s good faith and business. So, it 
may be said that if a product has a low reliability, that the dis¬ 
tributor and the manufactures are still liable - even if no harm or 
injury has occurred. As Taguchi and Clausing (17) have explained: 

“When a product fails, you must replace it or fix it. In 
either case, you must track it, transport it, and apologize for 
it. Losses will be much greater than the costs of manufac¬ 
ture, and none of this expense will necessarily recoup the 
loss to your reputation.” 

In general there are two aspects of reliability that concerns an 
electromechanical engineer. These two aspects are mechanical and 
electrical. We will talk about each to some extent. For now, We 
need to remember that mechanical reliability is greatly impacted by 
temperature and vibration that in general may lead to fatigue and 
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failure. Electrical reliability is influenced by temperature as well. 
Other factors impacting this aspect of reliability include current 
and voltage levels as well as electromagnetic or radio frequency in¬ 
terference. 

It should be mentioned that there are other factors impacting 
reliability. These include shelf life, corrosion, chemical reactions, 
solderability, moisture, aging, etc. A discussion of these factors, 
however, is beyond the scope of this book. 

Since reliability is in a way a study of failure, it is important to 
first discuss failure. Failure is defined as the inability of the system 
to meet its design objectives. It may be that the system was poorly 
designed and never met its objectives. Or, it initially met its design 
objectives but after some time, it failed. In these two scenarios, 
clearly certain factors were overlooked! 

There are three types of mechanical failures; namely, failure 
by elastic deflection, by extensive yielding or by fracture. A discus¬ 
sion of these subject generally starts with stress analysis followed 
by the assumptions made. 


Stress Analysis 

Stresses are internal distributed forces, which are caused by 
external applied loads. Strains are changes in the form under the 
same loads. Consider a rod of length L and diameter A. One may 
intuitively recognize that the displacement of the end of this rod 
depends directly on the magnitude of the applied force - very simi¬ 
lar to the force-deflection relationship of a spring-mass system: 

MM 

F 



Figure 8.1. Force-Deflection Relationship 


Now consider what happens inside of the rod in Figure 8.2: 
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Figure 8.2. Internal forces 


The concentrated load is (internally) developed over the area of 
the cross section. Thus, one may express this distributed force as 
follows. 


F 

<j - — 

A 

Similarly, a distributed (average) displacement may also be calcu¬ 
lated. 


A 

£ - — 

L 

It turns out that a and s have a 
similar relationship to a spring-mass 
force-deflection curve. The slope of this 
line (E) is called Tensile Modulus, 

Young’s Modulus, or Modulus of Elastic¬ 
ity. 

Consider another scenario. A block under a shear force will 
also deflect. In shear the force deflection relationships are defined 
as follows. 

F = tA , F ^ 


where A is the area, r is the shear stress. 
t = Gy , 

where G is shear modulus and y is shear strain. 

In general both normal as well as shear stresses develop in 
solids under a general loading. For example, a cantilever beam un¬ 
der a simple load at the free end, exhibits both normal as well as 
shear stresses. 
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Figure 8.3. Normal and Shear Forces Developed in a Cantilever Beam 


Considering that deformations depend on the general state of 
stress as a minimum both constants E and G must be known. This 
is only true if the material is isotropic, i.e., the material behavior is 
independent of the direction of applied loads. 

General Stress-Strain Relationship 

A small cube of a material under general loading exhibits the 
following stresses. 



Figure 8.4. General State of Stress 


The strains are related to stresses through the following ma¬ 
trix relationship 

{e} = [S]{o} 

Copyright © 2003 by Marcel Dekker,-Inc,-All.Rights -Reserved, •-- 




The matrix [ S ] has 36 constants (unknowns); however, it can 
be shown that these 36 unknowns may be reduced depending on 
material behavior. They are as follows: 

Anisotropic Behavior 

Materials such as bone exhibit different behavior depending 
on the loading direction and location. In this case there are 21 in¬ 
dependent constants (unknowns) in the [ S ] matrix. This matrix is 
symmetric. 


' £ x 


$l\ $12 ^13 ^14 ^15 ^16 



e y 


$22 $23 $24 $25 $26 



< 


$33 $34 $35 S 36 



ty* 


^44 $45 ^45 


x yz 

y xz 


Symmetric S 55 


T xz 

Yxy 


$66_ 


X xy 


Orthotropic Behavior 

Materials such as wood or composite lamina have different 
properties in two principle directions. These materials may be 
modeled with 9 independent constants (unknowns) in the [ S ] 
matrix. This matrix is symmetric. 


' £ x 


^11 

^12 

$\3 

0 

0 

0 ' 


'<*x 

£ y 


$\2 

$22 
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0 

0 

0 
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0 
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x yz 
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0 

0 

0 

0 

$55 

0 


Z xz 

Yxy 


0 

0 

0 

0 

0 

$66 _ 


r xy 


Isotropic Behavior 

Most engineering material properties are independent of di¬ 
rection. These materials may be modeled with only 2 independent 
constants (unknowns) in the [ S ] matrix. This matrix is symmetric. 
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Note that this relationship only requires the knowledge of the 
following are 3 constants; namely, E , G , and v. E is the modulus 
of elasticity and G is the shear modulus, v is Poisson's ratio de¬ 
fined as the ratio of the transverse strain in the j direction to strain 
in the i direction when the stress is in the i direction. All three 
variables may be determined experimentally, however, the following 
relationship exists among them. 

G = ——— and (-1<v<1/2) 

2(1+ v) 


The Stress-Strain Curve 

The choice of material is an important factor in the overall 
system design both from a budget point of view as well as struc¬ 
tural integrity. As it was pointed out, different materials behave 
differently particularly if their behavior depends on the direction of 
application of the load as in anisotropic or ortho tropic materials. 
This rule applies even to isotropic materials. If various specimens 
from different isotropic materials are loaded to complete fracture, 
three different behaviors are observed: 


1 . 

2 . 

3. 


Brittle Materials - The stress- strain curve is completely lin¬ 
ear as depicted in Ifigure 8.57 1 

Ductile Material with Strain Hardening - The stress-strain 
curve is nonlinear but the initial portion of the curve is lin¬ 
ear as shown in |figure 8.6 

Ductile Material without Strain Hardening - The 
stress-strain curve is nonlinear but the initial Portion of the 


curve may or may not linear (figure 8.7), 
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Rupture 


Stress 



Strain 


Figure 8.5. The Behavior of Brittle Materials 


Ultimate Stress 



Figure 8.6. Ductile Material with Strain Hardening Behavior 



Figure 8.7. Ductile Material Without Strain Hardening Behavior 

Examples of these types of materials in the electronics pack¬ 
aging are as follows. Most ceramics and certain type of printed cir- 
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cuit boards may be considered as brittle. Most metallic cases and 
chassis are ductile but exhibit strain hardening. Solder as well as 
most plastics are ductile but do not have any strain hardening. 

When dealing with ductile materials, there are several regions 
of interest. Knowledge of these regions enables us to verify whether 
a system has failed or has the potential of failure. These regions 
are defined as follows. 

1. Proportional Limit: The point on the stress-strain curve 
where the curve begins to deviate from a straight line. 

2. Elastic Limit: Maximum stress, to which a specimen may 
be subjected when upon removal of the load, no permanent 
deformation is caused. 

3. Yield Point: A point on the curve where there is a sudden 
increase in strain without a corresponding increase in 
stress. 

4. Yield Strength: The maximum stress that can be applied 
without permanent deformation of the test specimen. This 
is the value for which there is an elastic limit or at 2% de¬ 
formation. 

5. Ultimate Strength: Also called tensile strength, is the 
maximum stress value obtained on a stress-strain curve. 

Simplification or Engineering Assumptions 

It is now appropriate to review a set of simplifying assump¬ 
tions that may be applicable in about 85% of all stress analysis 
problems. The first set of such assumptions is based on the idea 
that material behaves linearly and proportional limit is not 
reached. The second set of assumptions is based on simplifying a 
three-dimensional geometry into a one- or two-dimensional case. 

Linear Elasticity 

If the loads are such that the stresses with the material never 
reach proportional limit and the following conditions are met, it is 
said that the system follows the theory of Linear Elasticity or 
Strength of Materials. The theory of Linear Elasticity dictates that 

> Deformations, strains and rotations remain small. 

> Stiffness through the model does not change. 

> Boundary conditions remain the same, e.g., loading direc¬ 
tion with deformation does not change. 

> Material remains in the linear elastic range. 
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This is a potential pit fall for a great deal of designers in using 
Finite Element Analysis packages somewhat as a black box. Once, 
any of these assumptions are violated, the solutions obtained must 
be questioned and other options exercised. 

Geometric Simplifications 

The general equations of elasticity are 3-D and complicated. 
However, under certain conditions, it is possible to reduce the 
three-dimensional stress-strain relationship to two dimensions, 
thus simplify the solution procedure greatly. These are 

> Beam, Plate and Shell Theories 

> Plane Stress 

> Plane Strain 

> Axi-symmetric 

1. Beams, Plates and Shells are areas which have been stud¬ 
ied for a long time and a well established set of formulae for calcu¬ 
lating their deflections, stresses and natural frequencies under a 
variety of loading and boundary conditions exists. For more infor¬ 
mation on these topics see (11). 

2. Plane Stress- In general, a structure is under conditions of 
plane stress if the stresses that develop along one of the three or¬ 
thogonal axis are so small compared to the rest of the stresses that 
they can be assumed to be zero. This generally occurs if the thick¬ 
ness of the problem is much smaller than other dimensions of the 
problem and all applied loads are in the plane of the problem, then 
a z = 0 and v yz = t x , = 0 . 


X 

Figure 8.8. Plane Stress Condition 

2. Plane Strain- In contrast to plane stress, if the thickness of 
the problem is much larger than other dimensions of the problem 
and all loads are applied uniformly in that dimension, then condi¬ 
tions of Plane Strain exist. An example of a plane strain problem is 
a cylinder under uniform pressure. Another example would be the 
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strip foundation of a building. All cross sections along the X-axis 
of the cylinder and the strip foundation that are away from the 
boundaries are under the same loading conditions and, therefore, 
experience the same deformation. As a result there will be an out- 
of-plane stress. In general a structure is under the plane strain 
condition if the strains that develop along one of the three orthogo¬ 
nal axes, say the x-axis are zero. The longitudinal stress along the 
X-axis is not zero and can be evaluated in terms of ct z and cr y , while 
the shear stress along the X-axis is zero. 



Figure 8.9. A Cylinder under Uniform Pressure; Plane Strain Conditions 


3. Axi-symmetric- If the problem is a body of revolution and 
the loading and boundary conditions are also symmetric around 
the same axis, then it is possible to solve the equations in 
2-dimensions 



Figure 8.10. Axi-symmetric Conditions 

In finite element analysis models that intersect the center line 
(the Z-axis) require special treatments. This type of models cause 
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computational problems and the nodes on the axis of revolution 
must have the radial translational degrees of freedom constrained. 

Stress Concentration 

In certain classes of problems, the stress field is uniform over 
the domain with an exception of isolated regions where it becomes 
complicated. An example of this problem is a simple cantilever 
beam loaded with a uniform load along its axis but with a notch 
near its end. Most engineers recognize this notch as a stress con¬ 
centration point and expect that it would generate higher stress 
levels. A notch is a common stress concentration point. Other 
common causes of stress concentration are as follows: 

> Abrupt changes in section geometry such as the bottom of 
a tooth on a gear. 

> Pressure at the point of application of the external forces. 

> Discontinuities in the material itself, such as non-metallic 
inclusions in steel. 

> Initial stresses in a member that result from overstraining 
and cold working. 

> Cracks that exist in members caused by part handling or 
manufacturing. 

In the theory of Strength of Material, a stress factor is pro¬ 
vided depending on the geometry of the point but further discus¬ 
sions are beyond the scope of this book and the reader is referred 
to other sources such as reference (26). It may be added, however, 
in practice, a stress concentration factor for circular holes in plates 
is approximately given by (18): 


x + 0.3 

where x is the ratio of width of the strip to the diameter of the hole. 

Failure 

Now that a discussion of Stress Field has been given, we may 
turn our attention to the concept of Failure. In developing a part or 
assembly, the designer must determine possible modes of failure of 
the system and then establish suitable criteria that accurately pre¬ 
dict the various modes of failure. Some categories of failure are as 
follows. 
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Failure by Elastic Deflection 

Failure by elastic deflection means that the deflections are too 
large to be acceptable but stress values are not important as de¬ 
sign criteria. For example, consider a skyscraper that can sway 20 
feet on the top of the building. Compared to the overall length of 
the building, 20 feet is quite negligible, however, imagine yourself 
on that last stoiy as the building begins to move. Not too pleasur¬ 
able! So the displacements while not causing any harm are con¬ 
sidered unbearable by the end user. Another example of undesir¬ 
able deflections is in vibration of system components. When the 
amplitude of vibration is large enough that parts collide the system 
has failed. Other examples involve structures such as beams or 
shells which may buckle under compressive loads. Under these 
circumstances, the structure regains its original shape once the 
loads are removed. 

Failure by Extensive Yielding 

Under such circumstances, the system exhibits permanent 
deformations. Even plastic flow may take place. These are most 
significant with regards to simple structural members such as axi¬ 
ally loaded members, beams, torsion members, columns, or possi¬ 
bly thin sheets or plates subject to in-plane forces at ordinary tem¬ 
peratures or at elevated temperatures. Thus stress values play an 
important role as a design criterion. 

Another phenomenon that must not escape the designer is 
creep, particularly if plastic parts are designed to carry loads. 
Creep is the action by which strains are developed without an in¬ 
crease in the loading conditions. 

Failure by Fracture 

There are three types of failure by fractures: Sudden with no 
evidence of plastic flow, fracture of cracked or flawed members, 
and progressive fracture better known as fatigue. 

Let us review fatigue briefly. The internal stress state of a vi¬ 
bratory system may have a microscopic impact on the load carry¬ 
ing capacity of the system. This causes the system to fracture at 
stress levels far below the yield values. Fatigue strength of a mate¬ 
rial may be altered by such factors as frequency of cycling, cold 
working of material temperature, corrosion, residual stresses, sur¬ 
face finish and mean stress. The relationship between the stress 
levels and the number of cycles to failure is generally depicted in S- 
N curves. 

The S - N curve usually reported for a given metal is often 
taken to represent a 50% probability of failure. The values above 
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this line present a higher probability of failure and values below 
this line presents a lower probability of failure. Furthermore, if 
failure occurs below 100,000 cycles, the fatigue is said to be low 
cycle. Above this value fatigue is considered to be high cycle. See 
Figure 8.11. 



Failure Criteria 

Based on the foregoing discussion, a failure criterion must be 
set for each mode of failure. 

Failure by Elastic Deflection 

The failure criterion is set based on the knowledge of the 
maximum deflections allowed and making precautions not to ex¬ 
ceed these values. Another criterion that is often neglected is 
bucking. A simple test is to identify the components which are un¬ 
der compressive loads and ensuring that these components have 
not exceeded their load carrying capacities. For example, for a sim¬ 
ple slender column, the maximum load is 

P 

r max *“* ^ 

L 2 

Where E is young’s modulus of elasticity, / is second mo¬ 
ment of inertia and L is the column’s length. 
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Failure by Extensive Yielding 

The criterion used to identify yielding depends on the material 
for a large extent. These are numerated below. 

1. Maximum Principal Stress (Rankine’s Criterion) is a good 
criterion for brittle fractures. The part fails when maximum 
principal stress has reached yielding (or rupture). 

2. Maximum Shear Stress (Tresca’s Criterion) is a good crite¬ 
rion for ductile materials. The part fails when maximum 
shear stress has reached yielding. 

3. Maximum Strain (St. Venant's Criterion) is also a good cri¬ 
terion for ductile materials. The part fails when maximum 
strain has reached yielding. It gives slightly more reliable 
results than Maximum Principal Stresses. 

4. Von Mises failure criterion is based on Strain Energy Den¬ 
sity of distortion being equal to energy of distortion at yield. 
This is a good criterion to evaluate stresses for fatigue. 

5. Octahedral Shear Stress is the same approach as the strain 
energy but the formulation is based on shear energy. Gen¬ 
erally speaking Tresca and Von Mises give the best results 


Failure by Fatigue 

The number of cycles to failure in a fatigue environment is es¬ 
timated from the materials S - N curve. This curve reported for a 
given metal is often taken to represent a 50% probability of failure. 
This means that to develop a better or worse confidence the value 
of stress value must be adjusted. For example, consider figure 


50% of the samples under testing will fail after lxlO 11 cycles. Now 
suppose that we need to be confident that only 10% of the compo¬ 
nents will fail under these conditions. Therefore, we need to make 
sure that the stress levels will not exceed 13,000 psi. Another sce¬ 
nario is that higher part failures may not present a problem. In this 
case, the stress levels may be increased to as much as 21,000 psi 
leading to a 90% part failure by the time lxlO 11 cycles has been 
reached. 

A good rule of thumb is to keep the maximum stresses in a 
vibratory system below 25 to 30 percent of yield stress to have an 
indefinite fatigue life. 




8.11. Under a stress level of 17,000 psi, there is a probability that, 
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Life Expectancy 

Miner’s Index is used to predict the vibration lifetime of a 
system that experiences different levels of frequencies (and corre¬ 
sponding stress levels) 






Nj is the cycles to failure at stress level i and n i is the actual 

cycles of vibration at that level. Theoretically, the value of R should 
be equal to 1; however, as Steinberg (15) points out this value is 
closer to 0.75 for electronics equipment which are employed in 
random vibration environments. 

Example 1 

A certain material undergoes a stress level of 10,000 psi at a 
frequency of 100 Hz. Considering that maximum number of cycles 
it can endure at this stress level is 2.6xl0 10 , determine the time to 
failure. 

Frequency is defined as number of cycles per unit of time 
(usually seconds). Thus to calculate the number of cycles for a 
given frequency and time, one may multiply frequency and time to 
obtain the number of cycles. 


R = 


n _ number of cycles 


N 


N 


(time)(frequency) _ tf 

n ”77 


0.75 = ■ 


f(100) 
10 


f = 1.95x10“ 


2.6x10 


That is 1.95 x 10 8 seconds or 54,167 hours. 

Example 2 

Suppose that the values in the previous example are la levels. 
How would the life expectancy change under random vibration? 
Assume that the following relationship exists between the stress 
levels and the number of cycles to failure: 


Ni = 


7.253x10 

11.36 


55 


a 
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It is noteworthy to consider that a single degree of freedom 
system (e.g. a spring-mass system) will only resonate at its only 
natural frequency. However, under random vibration, while the 
resonant frequency remains the same, the amplitude of vibration 
changes bringing about changes in stress levels. For this problem, 
we have 

7 253 x 10^ 

At IS: N, = —-—- = 2.63xl0 10 , n, = 100(time)(68%) = 68t 

(10,000) 1136 

/"»55 

At 2S N 2 = —- X - - -, - r = l.OOxlO 7 n 2 = 100(time)(27%) = lit 

(20,000) 1 36 

—i q 1 A 55 

At 3S JV 3 = - = 1.00x10 s «3 = 100(time)(4%) = At 

(30,000) 1136 

Now apply Miner’s rule: 


-^L + —+ = 0.75 

yv, n 2 n 3 


6 St 


21t 


2.63x10 


10 


■ + ■ 


At 


= 0.75 


r = 17563 


1.00x10 1.00x10“ 


That is 17,563 seconds or 4.88 hours. If the system were be¬ 
ing tested in a laboratoiy, this means that it would survive the 
shake table for 4.88 hours. 

Thermal Stresses and Strains 

As electronics equipment are operated, the internal tempera¬ 
ture rises to a steady state value. Once the equipment is shut 
down, the temperature is lowered to that of the environment. 

There is also true temperature transience as electronics 
equipment is operated. Suppose that the electronics is used for 
number crunching. As the CPU is engaged in this activity, its 
power consumption increases and there is a corresponding in¬ 
crease in temperature. As this activity reduces for I/O activities or 
once the calculations are done, the power consumption is reduced, 
thereby reducing the temperature levels. 

Materials generally expand (or shrink) as the temperature in¬ 
creases (or decreases). As components begin to expand at different 
rates, they “push” against each other leading to what is generally 
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known as thermal stresses. Thermomechanical management, thus, 
involves the impact of temperature change on material behavior 
and its internal state of stress and strains. 

In electronics enclosure designs, one must be aware of the 
thermal load variations that the electronics components may un¬ 
dergo and the impact it may have on the overall system. 

Thermal Strains and Deflections 

In a uniform temperature field, the behavior of a uniform 
component depends on the geometry. In general, three conditions 
may exist: 

1. No restrictions: The component is free to deform. While 
there are deformations, no significant stress state is 
formed. 

2. Constraints: The component may be constrained minimally. 
In this case, only deflections and deformations take place 
but again no significant stress state is formed. 

3. Properly Constrained: As the material expands, there is not 
enough “room to move.” Therefore, stresses and strains de¬ 
velop. 

However, if the temperature field is not uniform and a tem¬ 
perature gradient exists, then there is the possibility that various 
segments of the specimen deform at different rates caused by non- 
uniform temperature distribution. 

There is yet a third condition. In the last two cases, it was as¬ 
sumed that the material is uniform. Many components in elec¬ 
tronics equipment are nonuniform; each segment having a differ¬ 
ent rate of thermal expansion. Furthermore, many engineering 
plastics used today, exhibit a dependency on the magnitude of 
temperature. This combined with a nonuniform temperature dis¬ 
tribution can potentially develop into severe stress gradients. 

Basic Equation 

In a linear static problem, temperature changes affect the 
strains in the following way 


£ total £ mechanical £ thermal 


where £ thermal = a AT , a is the coefficient of linear thermal ex¬ 
pansion and AT is the temperature change from a stress free state. 
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This equation looks deceptively simple. A treatment this sub¬ 
ject in any detail is beyond the scope of this book. The interested 
reader is referred to Thermal Stress and Strain in Microelectronics 
Packaging (19). However, FEA must be employed for practical 
problems. 

In electronics packaging, the impact of thermal stress and 
strains is seen primarily at chip, component, and board levels. 

Die Attachments 

Residual thermal stresses are introduced in the cooling step of 
the bonding process due to thermal expansion and mismatch 
among the die, the bonding material and the package. These 
stresses may cause the die to crack. 

In this regard, voids existing in the bonding layers lead to 
stress concentration areas. These voids are generated due to a va¬ 
riety of reasons including trapped gas, liquid or other impurities. 
Other stress concentration areas may be caused by layer separa¬ 
tion due to improper bonding, fatigue, creep or rupture. In addition 
to stress concentration, voids may also increase the chip operating 
temperature and cause hot spots. 

Methods of improving Die-Attach quality include bonding in a 
pure environment, application of pressure for good contact, and 
back-grinding and wafer thinning of GaAs devices 

IC Devices 

There are two types of plastic packages. Thin elongated pack¬ 
ages with large chips known as thin small-outline packages 
(TSOP), and high-lead-count large square packages with relatively 
small chips, known as plastic quad flat packages (PQFP). Suhir 
(20) provides means of calculating the maximum warpage in these 
package types. 

PCB Warpage 

The sources of thermal stress may be numerated as follows. 

1. Mismatch of global thermal coefficient (TCE) between major 
components 

2. Local TCE mismatch between subcomponents such as lead 
and solder 

3. Lead stiffness 

4. Thermal Gradients in the System (i.e. Non-uniform Tem¬ 
perature Distributions) 

5. Unbalanced component lay-up 

6. Nonplanar boards 
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The following are a few simple rules to minimize printed wire 
board’s warpage. 


1) Geometry Symmetry 

a) Asymmetric lay-ups create layer warpage. Use symmetric 
configurations about the board mid plane. 

2) CTE match 

a) Consider materials that have similar CTEs. 

3) Thickness Tolerance 

a) PWB warpage is very sensitive to layer thickness toler¬ 
ance; control it rigorously. 


Some Tips 

Studying the impact of temperature on a piece of equipment 
can be as crude as blowing hot air on various segments of a work¬ 
ing PCB using an air gun. However, for a more detailed study: 

1. Consider the effects of temperature on properties. Do mate¬ 
rial properties vary greatly with temperature? 

2. Conduct a heat transfer analysis to develop a better under¬ 
standing of temperature variations. 

3. Where a range of properties are given, use both ends of the 
spectrum. 

4. Use the temperature field to assess a need for a thermal 
stress analysis. 
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9 

Electrical Reliability 


We, as engineers, program and/or QA managers, must con¬ 
duct reliability calculations to develop an understanding of failure 
issues as a basis for the first year and extended warranty sched¬ 
ules or as a basis for maintenance and repair scheduling. To do 
this, we need to determine the reliability of components individu¬ 
ally and the assembly as a whole. For this purpose, we need to cal¬ 
culate Mean Time To Failure (MTTF) or Mean Time Between Failures 
(MTBF). 

Electrical reliability calculations are simple but they can be 
quite sophisticated if statistical theories are employed. Herein, the 
most basic approaches will be presented. To do this, we need to 
become familiarized with some terms. For a more detailed discus¬ 
sion of this topic the reader is encouraged to read AT&T Reliability 
Manual (21). 

Cumulative Distributions Function ( F{t )) is the probability of 
system’s first failure before time t, 

0 < F(t) < 1 

Survivor Function is the probability of surviving to time t 
(without failure). 

S(t) = l-F(t) 


Hazard Rate (/.(t)) is the instantaneous rate of failure of a 
population that have surviv ed to time t . A typical hazard rate 
curve is shown in 


Figure 9.1. 
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Figure 9.1. Hazard Rate vs. Operating Time 


There axe three regions in this curve: 

1. Infant Mortality - The initially high but rapidly decreasing 
hazard rate corresponding to inherently defective parts. 

2. Steady State - The constant or slowly changing hazard rate. 

3. Wear Out - It generally occurs in mechanical electronic 
parts or when degradation exists. 

Hazard rate is commonly expressed in units of FITs. One FIT 
is equal lO 9 per hour. In other words there is a lO 9 probability of 
failure in the next hour. Military standards are different in that 
they use Percent Failing Per 1000 Hours of Operation. 

Note that these three variables, i.e., F(t), S(t) and Aft) are re¬ 
lated and if one is known, the other two may be calculated. Fur¬ 
thermore, time (t) is generally expressed in terms of hours. 

First-Year Failures 

AT&T Reliability manual (21) defines the probability of 12 
month survival as: 

S(t) = e~^ 

There are 8760 hours in a year. So this equation may be alterna¬ 
tively rewritten as 

S(8760) - e~ 8760A 

Mean Time To Failure is defined as MTTF — — . 

A 

Probability of Failures is F{t) = 1 - S(t ) 
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Example 1 

A device has a 250 FIT hazard rate, determine the MTTF and 
the survival rate for the first year. 

A = 250xl(T 9 = 2.5 xlO -7 

MTTF = — =- l - - 

* 2.5xl0“ 7 

MTTF = 4,000,000 hours 
5(0 = e~ M 

5(8760) = e ~ S760{2 - 5 xl0_7) 

5 = 0.998 
F = 1-5 = 0.002 

Therefore for this device the survival rate is only 99.8% in the first 
year and only 0.2% failure. 

Example 2 

A device has a MTTF of 250,000 hours. What is the first year 
failure rate? 


MTTF 250000 
A = 4x10 -6 or 4000 FITs 

5(8760) = e-^ 876 °) - e - 4 ^0- 6 (8760) 

5 = 0.965 
F = 1-0.965 
F = 0.0344 or 3.44% 

In other words, a MTTF of 250,000 hours represents a 3.44% first 
year failure! 

Reliability Models 

The model employed in the previous two examples is called 
the Exponential Distribution model. Along with this model, there are 
two other distribution models that are commonly used. The first is 
called the Lognormal Distribution model and is generally used for 
accelerated testing of semiconductors. The second is called the 
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Weibull Distribution model and is widely used particularly for infant 
mortality calculations. In this model, the hazard rate varies as a 
function of usage. 

Weibull Distribution 


A(t) - A a t a 2 0 > 0, 0 < a < 1, t> 0 

- A t l ~ a 

S(t) - exp(—-) 

1-a 

Note A a is a scale parameter and is the hazard rate at one-hour de¬ 
vice usage. A long-term hazard rate is assumed at 10,000 hours. 
Furthermore, Weibull reduces to the Exponential Distribution 
model for a = 0. 

Example 

A device has a 250 FIT long term hazard rate, using the 
Weibull model determine the survival rate in the first six months. 
Assume that a=0.75 

A = 250x1 CT 9 - 2.5xl0~ 7 

This is the long-term hazard rate. To recover the initial rate, use 
the following equation. 


A{t) = Aj~ a 

2.5 xl0~ 7 = A c (10000-° 75 ) 
A a = 0.25 x 10 -4 


- ) / ,_a 

S(t) = exp (—^—-) 


£(43 80) = exp| 
S = 0.9992 


1 - a 

^-(0.25xl0“ 4 x4380 (1 ”°- 75) ) 


V 


1-0.75 


J 


F{t) = \-S{t ) 
F - 0.0008 
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Therefore for this device the failure rate in the first six months is 
only 0.08 %. 

Temperature Effects 

Temperature has an adverse effect on reliability. Reference 
temperature is 40° C and the corresponding hazard rate is denoted 
by subscript r. If the operating temperature is different than 40° C, 
then 


X u ~ AX r where 


4=expl 



f 




_l_ 

T 

1 r 



E a is the activation energy. K b is the Boltzman constant. Tem¬ 
peratures are in absolute values. 

Accelerated Testing 

A practical application of this phenomenon is in accelerated 
testing. Many devices can not be tested under their operating con¬ 
ditions to failure because of time scales involved. For example, 
some medical devices are required by the FDA to have a life ex¬ 
pectancy of several decades. The only practical means of demon¬ 
strating this life expectancy is through accelerated testing. 


AF = M7TFnomial - exp^- J-!— 

MTTF acce ] K b L Formal ^accel. . 

To find the Find MTTF under normal operating conditions; 
first, find the MTTF at two different accelerated temperatures. This 
allows the activation energy to be calculated. From these values the 
MTTF at normal conditions may be found. 

Electrical Stress Effects 

Other factors such as electrical stress may have an adverse 
effect on the device reliability as well. A component is generally 
rated either for a specific voltage, current, power, etc. Reliability 
studies have shown that once a particular percentage of this rating 
is exceeded, the hazard rate is increased. In MIL-HDBK-217F (22), 
this percentage is assumed to be 25% and is denoted by Po. 
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For electrical stress, MIL-HDBK-217F (22) provides the fol¬ 
lowing relationship. 

X u - A E X r where 
A e =exp[w(P 1 -P 0 )] 

Po is 25%. 

Pi and m are provided by MIL-HDBK-217F and Bellcore stan¬ 
dards (23). Per Bellcore guidelines, each component is rated at a 
particular electrical stress level. At that prescribed value Ae is 
equal to 1. Now, should the component be operated at stress levels 
below the set value its reliability does not change; however, should 
the operating stress goes beyond this limit, its reliability decreases. 

System Failure Rate 

To calculate a system reliability or hazard rate, we need to 
have a knowledge of the failure rate associated with the compo¬ 
nents and then a knowledge of the system itself. 

Component Failure Rate 

The component failure rate is simply a product of a hazard 
rate multiplied by various stress factors; namely electrical, tem¬ 
perature and quality - another factor undefined as yet. Quality 
factor takes into account manufacturing issues. Generally this 
factor is equal to 3 for commercial items, 2 or 0.9 for parts which 
are to be used in military applications. 

A SS i " k G i n Q i n S i n T i 

X is the failure rate and n denotes a stress multiplier. G, de¬ 
notes generic values, Q l is the quality factor as defined previously, 
Sj is electrical stress factor =1 for 50% stress, and 7} is the thermal 
stress factor = 1 for 40° C, 

System Failure Rate 

A system may be comprised of a group of subsystems; each 
made up of a many circuit boards. The circuit boards in a subsys¬ 
tem and indeed, the subsystems within a system may be arranged 
in series, parallel or a complex arrangement. A discussion of par- 
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allel and complex arrangements is beyond the scope of this work 
and the reader is referred to other sources such as (21). 

For a series network, the reliability of the system may be cal¬ 
culated as follows: 

n 

J-ss = n eY, n i^sSj > 

i=\ 

Where n E is the environment factor and N, is the number of each 
component in the system. 

This approach provides the strictest sense of reliability calcu¬ 
lation in which the failure of any component will flag the entire 
system as failed. 

To calculate the hazard rate for a system, tabulate the follow¬ 
ing items for each component: 

1. generic hazard rates 

2. thermal stress factor 

3. electric stress factor 

4. quality factor 

The steady state hazard rate is obtained by multiplying the 
above numbers. To obtain the failure rate of the entire system, add 
component hazard rates and multiply by an environmental factor. 
The values of various stress and environmental factors are avail¬ 
able from various sources such as Bellcore or other standards. 

Example 

Determine the failure rate and mean time to failure of an 
electronic sensor. This system is designed to be used only one hour 
a day. Next evaluate the first year repair volume should there be a 
production level of 12345 units. 

X = 2.28xl0 -5 

5(0 = e’** 

5(8760) = e -8760(2.28xl(T 5 ) 

5 = 0.819 

F = 1-5 = 0.181 
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This figure indicates that an 18% first year failure rate may exist. 
This would be true only of the system was operated continuously. 
However, the duty cycle is only one hour a day. Thus a new ad¬ 
justed MTTF must be found based on the duty cycle. In this par¬ 
ticular example each hour of operation is equivalent to one day. 
Thus, the adjusted MTTF must be calculated by multiplying the 
calculated MTTF by 24. 


Table 9.1. Hazard Rate Calculation Data 


Item 

Qty 

Description 

generic 

qual¬ 

ity 

elec¬ 

tric 

ther¬ 

mal 

Multiply 

All 

1 

13 

O.luF Cap. 

1.00E-09 

2 

1 

0.9 

2.34E-08 

2 

1 

IOuF Cap. 

1.00E-09 

2 

1 

0.9 

1.80E-09 

3 

2 

22pF Cap. 

1.00E-09 

2 

1 

0.9 

3.60E-09 

4 

1 

Diode 

3.00E-09 

3 

1 

0.8 

7.20E-09 

5 

1 

Network 

Res. 

5.00E-10 

3 

1 

1 

1.50E-09 

6 

4 

1.0M Res. 

1.00E-09 

3 

1 

0.6 

7.20E-09 

7 

1 

39 K Res. 

1.00E-09 

3 

1 

0.6 

1.80E-09 

8 

2 

470 K Res. 

1.00E-09 

3 

1 

0.6 

3.60E-09 

9 

3 

Potenti¬ 

ometer 

1.70E-07 

3 

1 

09 

1.38E-06 

10 

2 

10 K Res. 

1.00E-09 

3 

1 

0.6 

3.60E-09 

11 

1 

IC- 

Microship 

1.00E-05 

3 

1 

07 

2.10E-05 

12 

1 

LCD 

3.00E-09 

3 

1 

0.4 

3.60E-09 

13 

1 

Sensor 

2.50E-08 

2 

1 

1 

5.00E-08 

14 

1 

Battery 

1.00E-07 

3 

1 

1 

3.00E-07 

15 

3 

Connector 
(per pin) 

2.00E-10 

3 

1 

0.6 

1.08E-09 

16 

3 

Switches 

1.00E-08 

3 

1 

0.6 

5.40E-08 


Total =2.28E-05 


Environmental Factor = 1.00K+00 
System Failure Rate =2.28E-05 
Mean Time To Failure =43784.03 hours 


/Vr7YF a dj ustc d — 24 /V/77'F ca ] cu]ate d 
A/7TF adj usted =1050816.72 hours 

1 1 1 
adjusted M777 r adjusted 1050816.72 

^adjusted ~ 9.52xl0“ 7 

5(0 = 

5(8760) = g _876 °0 ) - 52xl0_7 ) 
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S = 0.992 
F = 1-5 = 0.008 


This indicates that only 0.8% of production or 99 units will likely 
fail in the first year. 
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10 

Some Analysis Tips in Using Finite 
Element Methods 


In this chapter, some analysis tips are presented. 

Plastics 

When using FEA to conduct an analysis of plastics, follow the 
general rules. 

1) Make sure to use appropriate elements. 

i) Generally small deflections /rotations assumption is over¬ 
looked. 

2) Make sure that “enough” elements have been used in the 
critical areas of the model. 

3) Many Plastics exhibit creep, i.e., prolonged application of 
force changes properties. 

Range of Material Properties 

Material properties often are given as a range. Or in other 
cases, component masses or dimensions may vary. A combination 
of these variations may bring about changes in the system’s re¬ 
sponse to external conditions that may be catastrophic. 

Therefore, it is prudent to conduct the analysis for these vari¬ 
ous combinations to ensure the integrity of the final production. 

CAD to FEA Considerations 

One must differentiate between the function of CAD programs 
and FEA Software. CAD packages are design tools developed to give 
a “touchy-feely” taste to the designer. FEA software is designed to 
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make calculations based on exact data. Thus, the two have very 
different roles. 


CAD to FEA, DO’s and DO NOT’s 

The following is a list of do’s and do not’s in CAD to FEA File 
translations. 

1. Do not just bring your CAD over to FEA. 

2. Create a file to be modified for FEA. 

3. Do you really need to simulate the entire model? 

4. Do you really need the degree of detail in regions of rela¬ 
tively no importance? 

5. Take out the extra’s before translating. 

Criteria for Choosing an Engineering Software 

Choosing engineering software for today is also choosing soft¬ 
ware for tomorrow: and like toothpaste out the tube, it is a difficult 
decision to reverse if it is wrong. There are many choices. How can 
one make an intelligent estimate and choose the right package for 
his/her application when such a large number of products are 
available? What are some of the related issues once the code is 
chosen and paid for? 

With the advent of the new generation of computers and the 
speed with which they operate, it is only natural for this tool to be¬ 
come common place in the world of engineering design and analy¬ 
sis. Up to about ten years ago, most engineers wrote their own spe¬ 
cialized programs to do certain calculations. General purpose 
commercial engineering programs were not commonly used. They 
were usually utilized by highly educated engineers in Research and 
Development areas for highly complicated tasks. The most effective 
means for marketing of these programs at that time was word of 
mouth! 

Engineering applications and computers have come a long 
way. On the one hand, personal computers and low-end worksta¬ 
tions have made it easy for design and manufacturing engineers to 
access these sophisticated programs more readily; and, on the 
other hand, the software vendors have recognized the potential of 
this market and have attempted to make their programs more 
user-friendly. In 1995, the engineering software industry was esti¬ 
mated to be around 50 billion dollars. Many new companies have 
entered this market and right now there are hundreds of programs 
available for the engineering community. 
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The following is a methodology for making software choices - 
take the points that work for you. These are summarized in six 
steps at the end of this chapter. 

What types of Engineering programs are there? 

By and large, there are four types of engineering software 
available today. These general classification are computer aided 
design and manufacturing (CAD/CAM), finite element analysis of 
stress and heat conduction problems (FEA), computational fluid 
dynamics (CFD), and finally specialized areas such as electromag¬ 
netic or combustion. I will examine each category briefly. 

Computer Aided Design (CAD) and Computer Aided 
Manufacturing (CAM) 

These software programs aid the engineer in developing drafts 
and blueprints. They are usually designed to produced three-view 
and isometric representations of objects with the given dimensions 
and special considerations. These drawings can then be used in 
manufacturing the piece. The same CAD programs may be used in 
Construction and Architecture in order to convey a sense of the 
building to be designed. Traditionally, these packages are for 
graphics alone and do not offer any analysis capabilities. More re¬ 
cently, however, these programs offer solid modeling capabilities 
which can be used to transfer model information to FEA packages 
for analysis. 

Solid Mechanics and Stress Analysis - FEA Software 

These programs actually do the various and necessary calcu¬ 
lations to evaluate whether or not a structural design would 
maintain its integrity through its life cycle. Historically, this is 
where engineering applications of computers and programming 
began on a wide scale. The matrix analysis techniques of structural 
engineering opened the ways to the applications of finite difference 
and finite element techniques. In return, these techniques allowed 
engineering simplifications of realistic systems and before long, 
FEA became a common place technique. Soon after, it was realized 
that these very techniques are also applicable in other areas of en¬ 
gineering and physics. Nowadays, FEA applications in solid me¬ 
chanics range anywhere from linear static to transient nonlinear 
analyses, to thermal stress and conduction calculations, and from 
simple impact problems to crash analysis. 
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Fluid Mechanics and Heat Transfer - CFD Software 

Once enough experience was gained in solid mechanics, at¬ 
tention was focused on fluid flow, heat and mass transfer type of 
problems. The CFD software may employ one of several analysis 
techniques available; namely, control volume techniques, finite 
difference techniques, and finite element methods. To this day, 
each formulation has shown its own unique characteristics and to 
my knowledge, no one is shown to be superior to others in general. 
There is a range of problems that are solved using CFD. Examples 
may be cited from simple incompressible viscous flow through 
pipes and pressure drop calculations to heating, ventilating and air 
conditioning to complicated supersonic flows. 

General Physics - CFD and/or FEA 

As analysis becomes a more integral segment of industrial de¬ 
sign cycle and as the computing power of modem computers in¬ 
creases almost exponentially, engineering simplifications and as¬ 
sumptions begin to take new forms. While in the past, coupling of 
various disciplines of physics were ignored, they can now be con¬ 
sidered and their effects be taken into account with the advent of 
new computational techniques and engineering software programs. 
By no means is the application of numerical methods limited to 
obvious applications of solid mechanics, fluid mechanics, or heat 
transfer. Various disciplines such as Magnetics, Chemical Reac¬ 
tions, Combustion, Molding, etc. have found their ways into vari¬ 
ous computer programs as well. 

There is a price to pay with these engineering advances. True, 
many complicated calculations may be done easily now; however, 
as the physics becomes more difficult, the numerical treatment 
becomes more delicate. For example, CFD problems sire by far 
much more difficult to solve than their FEA counterpart. It is esti¬ 
mated that 75% of all solid mechanics, i.e., stress analysis prob¬ 
lems and only 10-15% of CFD problems are linear. In a linear 
problem the mathematical equations do not pose any difficulty and 
they can be solved in a straightforward manner. However, if a 
problem is nonlinear, special considerations need to be made so 
that these mathematical equations can be solved. There are other 
added complications. As the physics becomes more complicated 
smaller elements are needed in order to capture and resolve 
changes in the primary variables such as velocity components, 
turbulence energies, or magnetic fluxes. These naturally result in 
much larger problem sizes requiring larger computer disk spaces 
and CPU time allocations. 
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Which software should I choose? 

First determine your general area of application, i.e., solid 
modeling, solid mechanics, fluid dynamics, and so on. Then, to an¬ 
swer this question adequately, one has to research the following six 
areas: applications, the completeness of the package, manuals, 
customer support, software quality assurance, and, finally, inter¬ 
facing to other programs. 

Application 

A client once asked me to help him choose the proper package 
for his application. My first thought was that he should go with 
package A because it was the leader in the industry in stress 
analysis - with a yearly lease price of $10,000 - as well as package 
B, another industry leader in fluid flow analysis - with a price tag 
of $17,000. Upon examining his application, I realized that he 
needed a nonlinear stress analysis solver with occasional applica¬ 
tions of fluid mechanics and heat transfer. I suggested that based 
on his applications, he should choose package C, because this 
company offered their product on a modular basis and he only had 
to pay for what he was going to use. True, companies A and B were 
industry leaders, but the user would end up paying for many fea¬ 
tures - such as crack tip propagation, stamping, nonlinear plastic 
flow, and so on - that would not have been used at all. The neces¬ 
sary modules of Software C cost this client about $7,000. 

When choosing a software package, it is important to know 
the extent of your applications and the physics involved. This, I 
believe, is one of the most important factors for a decision. And it is 
important to remember that some software vendors sell their pack¬ 
age as a whole and others in modules. 

Application, however, is not the only issue to be considered. 

Complete Package 

It is a normal assumption to think that the program I am 
about to buy or lease comes all in one package! That is not so. 
There are many “solvers” out there that can solve some veiy diffi¬ 
cult problems, however, they do not have any modeling capabili¬ 
ties, nor any tools to look at the results once the analysis is done. 
These programs generally depend on information provided by a 
third party software in order to work - without which, they are use¬ 
less. One has to be mindful of this fact particularly when compar¬ 
ing prices. 
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User Interface, Manuals, and Training 

As general-purpose commercial engineering codes are be¬ 
coming more and more available to nonexpert users, user inter¬ 
faces play more important roles, and code developers dedicate 
larger portions of their resources to user-friendly interfaces. Four 
levels of user interface exist: 

1. Graphical User Interface (GUI) where information is input 
through the usage of some picking device such as the 
mouse. The more sophisticated the GUI, the less usage of 
the keyboard 

2. Textual User Interface where the information in input pri¬ 
marily through the keyboard. 

3. Manuals, training courses and notes such as tutorials. 

4. Customer support. 

These factors are intertwined and need to be looked at closely. 
The following are only a sample of questions to be entertained. The 
question of user support will be treated separately. 

1. How user-friendly is the GUI (text base interfaces are rap¬ 
idly becoming a past trend)? How well is the information 
and logical flow of modeling laid out on the screen? How 
much access to manuals does one need to start using the 
program? 

2. How user-friendly are the manuals? Are they written as a 
reference or as tutorials? Would the vendor send copies for 
evaluations to your site? 

3. Are there specific manuals to teach the new user how to get 
started and come up to speed in a short period of time? 

4. How inexpensive and accessible are the training classes? 
One point to keep in mind is that some "Introductory 
Classes" can be so simple that one would wonder "why pay 
several hundreds of dollars just to sit in the class!" Look at 
the course syllabus before actually registering. Another is¬ 
sue is that some companies do not allow for cancellation 
once the money is paid! Or they have a surcharge for re¬ 
scheduling. 

Customer Support 

Support is an issue that is difficult to evaluate before buying 
the code. Yet, it is an extremely important issue because users - 
even pros - come across obstacles at times that they can not re- 
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solve. These obstacles range from simple negligence on the part of 
the user or not knowing portions of the code to discovering an un¬ 
detected bug. In a way, support engineers are the back bone of any 
software (or even hardware) company. The questions to research in 
this area are: 

1. The actual number of support engineers - are the code de¬ 
velopers also the ones who support it? 

2. Their level of competence - to be a good support engineer 
not only requires the ability to "see" the user's model and 
what causes the obstacle, but also the technical compe¬ 
tence to be knowledgeable in all of the disciplines that the 
code has to offer. Many times a user calls for support with 
an issue that has to do with the physical nature of the 
simulation and not the software. The support engineer 
must be knowledgeable of the assumptions and limitations 
of the code. 

3. Friendliness - rapport which eases communications. 

4. Timeliness - how quickly do they respond to phone calls, 
faxes and e-mails? Is there any type of procedure that over¬ 
sees these activities and their timely responses? 

The best technique to find the answers to these questions is to 
talk to other clients and users of that software and network with 
them. 


Quality Assurance (QA) 

QA seems to be the “buzz” word these days. Many companies 
pride themselves for having a QA department or a set of QA proce¬ 
dures in place. Some software companies have even taken the ini¬ 
tiative to become ISO certified. It may be that for the consumer, it is 
more important to talk about Software Quality Assurance (SQA) 
than QA. SQA depends on maintenance, error reports and bug fixes, 
user interfaces, and Finally developing new features. Software Per¬ 
formance Testing (SPT) may be a better term for what is generally 
known as QA. 

In general, the user's encounter with SQA is mainly through 
maintenance, error reports, and bug fixes. Furthermore, the major¬ 
ity of users make the assumption that SPT has been successful 
and never ask for verification of the code prior to purchase. 
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Maintenance 

From a user's point of view, maintenance means upward and 
platform compatibility, and life cycle. It implies that once a feature 
is developed, it works equally the same across any computer plat¬ 
form and would produce identical results. Moreover, updates 
should support the capabilities available in previous versions. Life 
cycle, here, refers to the number of updates that are supported in 
one major release. As an example, suppose that version 1.0 of a 
program has just been released, that is VI.0. The number 1 refers 
to a major release and the number 0 refers to the updates. Say, 
after six months, new updates - either features or bug corrections - 
are released. This new program may be called VI. 1 because noth¬ 
ing major has been done to the program and all changes are con¬ 
sidered to be minor. The next release might be VI.2, VI.3, and so 
on. Then, the software company may decide to make major 
changes and come out with version V2.0. Usually, the input files 
used for VI.x versions will not be accepted by V2.0 major release. 
The time span from V1.0 to V2.0 is called a life cycle. 

Error Reports and Bug Fixes 

Each software developer must have a comprehensive error re¬ 
porting system with different classes of error identifications. This is 
a measure of how well the program developers and SQA keeps 
track of bugs. If such a system is not in place, the same errors are 
likely to happen time and again. These errors and their implica¬ 
tions must be reported to the users after a formal investigation and 
verifications. And for specific cases, solutions and/or work around 
should be offered. Ask about it, if the marketing people have not 
clarified this point. 

Software Performance Testing 

From a user's point of view SPT is a set of problems, which 
should check every aspect of the program for the following: 

1. Backward compatibility. That is to say that the new version 
of the code should easily accept the previously defined in¬ 
put files and produce the same results - with the exception 
of error corrections. 

2. If the program allows for two ways of solving the same 
problem - such as different techniques for solving equations 
- they produce the same results. 

3. The mathematical limits and assumptions of various capa¬ 
bilities are tested. 
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Validation 

At times, validation problems are presented as QA results 
taken for SPP. Validation issues are different than software testing 
in that they try to show that the program produces results that are 
the same as those obtained by other means. There are three types 
of validation studies: 

1. Comparison with Exact Results. Unfortunately, exact re¬ 
sults are only available for a very small group of simple 
problems. Usually one-dimensional problems belong to this 
category. 

2. Comparison with Experimental Results. Traditionally, this 
has been the most reliable and most effective approach to 
validating any theory. But it, too, has its drawbacks for use 
as a technique for validating numerical works. Errors are 
almost always expected and most numerical techniques 
show "good" agreement with the experiments and not nec¬ 
essarily an "exact" agreement. One has to exercise care not 
to denounce any code or algorithm on this basis because 
many mathematical models have been developed to resem¬ 
ble the physical model under specific conditions. If these 
models are utilized in a way contrary to their basic as¬ 
sumptions, the results will not be dependable. 

3. Comparison with Other Numerical Solutions. This is 
probably one of the most common methods of validating 
any one code on different types of physics. Somehow it is 
believed that if the code compares with some other code, it 
renders the code validated! In reality, it only proves that 
there are no programming errors. 

In the final analysis, the one who truly validates the code is 
the user. He or she validates the code by utilizing the program in 
his/her real life engineering problems time and again. And he/she 
trusts that it produces reliable results for use in his analysis. 

Interfacing to Other Programs 

In today's engineering environment, and with the recognition 
of concurrent engineering, it is essential that various engineering 
software - even competing ones - have the ability to exchange data. 
It is quite common that a designer would create a model using a 
CAD program and then send it to the analyst for further work. 
Generally, neutral files such as IGES (Initial Graphics Exchange 
Specification) would allow for such information to be transferred. 
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However, other companies have collaborated to make this transfer 
much easier. It is important to know of such capabilities, if any. 
The following is a set of probing questions to ask: 

1. In case files are imported from CAD programs: Does the 
translator recognize both the wire frame and solid models? 

2. Does the user have a choice to indicate which entities are 
translated? This is an important option particularly when 
CAD models are exported to FEA and CFD programs. CAD 
models usually contain symbols (such as dimension lines 
as well as dimensions) that if brought into the analysis 
package introduce meshing difficulties. These "extra" have 
to be deleted before any finite element discretization can 
take place. 

3. Does the data flow in one direction or is it both ways? In 
other words, if software A is capable of recognizing files 
from other programs, how many other programs would rec¬ 
ognize files produced by software A? 

Conclusions 

It is not easy to decide which software to buy. However, one 
has to review three levels. First, one has to define the scope of the 
work to be done by the software, second, the duration of the work 
and finally the extent of usage. Considering these three and com¬ 
bining that with the available budget and knowing what to look for 
in each product, making the decision may become easier. So, here 
is the 6-step formula 

1. Determine the extent of your application. 

2. Decide whether you need a complete package or your appli¬ 
cation can be satisfies with certain modules. 

3. Research the type of training that the vendor provides to 
help you get started. Next, try to develop a feel for the clar¬ 
ity of the manuals. Then, find out how easily the program 
can be used. 

4. Contact other users and ask them how professional the 
customer support team is. 

5. Research the QA techniques and procedures used by the 
vendor and the extent of the tested capabilities. Does QA 
monitor even the manuals and the training courses? 

6. Finally, find out how easily, this software interfaces with 
other programs. Are there specific interfaces or should gen¬ 
eral neutral files (such as IGES) be used. 
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Design Considerations in an Avionics 
Electronic Package 


In this chapter, the mechanics and thermal design procedure 
and calculations (based on both analytical means and finite ele¬ 
ment methods) of an airborne ARINC 746 Cabin Telecommunica¬ 
tions Unit (CTU) are discussed along with the relevant FAR’s (24), 
as well as ARINC standards (25). Finally, reliability and failure 
rates of components and their effects on the performance of the 
system are studied. In designing avionics electronics packages, 
four areas of engineering need to come together and work concur¬ 
rently as a team. These are electronic, electro-mechanical packag¬ 
ing, thermal and finally vibration. Often times, thermal and vibra¬ 
tion analyses are done only as an after-thought through experi¬ 
mental means to ensure Federal Aviation Authority (FAA) Federal 
Aeronautics Regulations (FAR) and/or ARINC standards. 

Equipment design is influenced by its operational environ¬ 
ment and its impact on susceptible components. For example, if a 
system is operating in sun light, it can experience as much as 350 
watts of heat loading per square foot of exposure. It may also be 
that the same system is experiencing random vibrations due to a 
variety of sources adding to the severity of stress field especially at 
critical frequencies. The same considerations must be made for 
EMI and other environmental issues. In short, a system can not be 
designed adequately by ignoring the operational environment. In 
the avionics industry, regulatory agencies and standards organiza¬ 
tion impose rules such as Federal Airworthiness Regulation 


This segment was completed with the help and input of Mr. Robert E 
Walter. 
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(FAR)(24) and ARINC(25) as well as Radio Technical Commission 
For Aeronautics (RTCA)(26) have developed sets of test procedures 
to enable design engineers to test and verify their designs and pre¬ 
vent any safety-threatening failures. Much too often, however, 
thermal, mechanical and reliability evaluations have been left up to 
these tests and as the last item on the design agenda with no 
analysis taking place. This practice relies heavily on the designer’s 
previous experience and lack of any substantial changes to new 
systems. Otherwise, time- and budget-consuming design cycles 
must be undertaken before the product may be released to the 
market. 

There is another drawback to testing without analysis. These 
tests only reveal whether a system fails or not. It does not give any 
substantial information in regards to mean-time-before-failure 
(MTBF) reliability, mean-time-to-repair (MTTR) and so forth. These 
issues are important when one considers business issues such as 
scheduling repairs, product warranties, or consumer satisfaction 
with the product. 

Herein, analysis of a hypothetical commercial aviation elec¬ 
tronic package is discussed. The particular system is assumed to 
be an ARINC 746 (25) Cabin Telecommunication Unit (CTU). 

DESIGN PARAMETERS 

The following is a set of design parameters developed based on 
various regulations and standards. These areas cover operational 
characteristics, electrical and mechanical design. It should be 
noted that the design of electronics, while not discussed here, is 
within the scope of operation with its constraints set forth within 
the electrical design parameters. Packaging aspects are within the 
realm of mechanical design parameters. 

Operational Characteristics 

The model Cabin Telecommunication Unit (CTU) is similar to 
a public branch exchange (PABX) computer. It generates com¬ 
mands for the establishment and termination of air/ground, 
air/SATCOM and seat/seat communications. Furthermore, the 
CTU has interfaces to the aircraft’s networks and to the radio tele¬ 
phone control computer. It is also linked to embedded processors 
in the radio allowing user’s to obtain air-to-ground modem links. 
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Reference Documents 


To gain FAA certification and industry acceptance for this de¬ 
vice, compliance with the following reference documents is man¬ 
datory: 


• ARINC 746 - Cabin Communication System (CCS) 

• ARINC 600 - Air Transport Avionics Equipment Interfaces 

• RCTA/DO-160C - Environmental Conditions and Test Pro¬ 
cedures for Airborne Equipment 

• FAR Part 25 - Federal Aeronautics Regulations 

Electrical Design Specifications 

The specifications of the electrical design must be as follows: 

• The AC input is to be protected using a single pole circuit 
breaker. 

• The DC outputs of the power supply are to be protected 
using fuses. The DC output return is not permitted to float 
above ground potential. The high voltage DC side of the 
power supply is to be tied to a power interrupt capacitor. No 
batteries are to be used for power interrupts unless the bat¬ 
tery meets FAR Part 25, or a special containment is used. 

• The CTU must incorporate an internal chassis/safety 
ground independent from the return line per ARINC 600 and 
746 techniques for grounding and bonding. 

• The CTU must incorporate a single point ground as de¬ 
scribed in ARINC 600. All PCB’s should have an internal 
ground plane. This ground plane must be tied to ground at a 
single point. PCB standoff locations are not to be tied to 
chassis ground, and thereby defeating single point ground¬ 
ing. 

• In case of loss of power, the power supply and capacitor 
combination shall be capable of maintaining the output load 
for 2 seconds, a time in excess of RTCA/DO-160C and 
ARINC 600 specification. The loss of power is defined as out¬ 
put drops of 3% below nominal values. 

• The CTU must comply with RTCA/DO-160C for EMI/RFI. 
Experience has shown that all wires are to be shielded 
and/or twisted. All inputs/outputs should either be isolated 
or filtered. ARINC 600 connectors should be used. 
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Mechanical Design Specifications 

The specifications of the mechanical design must be as fol¬ 
lows: 


• The unit must weight less than 20 pounds. 

• ARINC 600 standards dictates a size limit of 6 MCU\ 7.50” 
(width) x 12.76” (length) x 7.64”(height) (WxLx H). 

• It is imperative that the CTU does not create any fire or 
smoke under any circumstances considered to be a hazard 
to flight operations. It is strongly recommended that all wires 
conform to MIL-W-22759, MIL-W16878E and cables to MIL- 
C-27500. It is against FAA rules to use any cables jacketed 
with PVC. All PCB’s are to be FR4-G10 or better. 

• The packaging sheet metal must be aluminum. No fractures 
or tears are permitted on formed sheet metal parts. The as¬ 
sembly must retain its mass during sudden deceleration - 12 
g’s per FAR part 25 and RTCA/DO-160C crash safety regu¬ 
lations. 

• All screws, washers, and nuts must conform to military 
standards. All fasteners must have a self locking device 
such as nylon or deformed threads. All PEM style inserts 
must be self locking style and only installed after irradiating 
and plating. 

• The temperature range is from -15°C to 55°C and altitude of 
zero to 15,000 feet. No temperature derating is permitted. 

• The CTU will be tray mounted. Air will be moved through 
the tray metering plate and enter the CTU from the bottom 
and exit from the top. The bottom and top surfaces will have 
several hundred 4 mm diameter holes. Eveiy effort must be 
made to eliminate other air leakage points. ARINC 600 dic¬ 
tates that an easily removable internal air filter be installed 
that filters particles down to 4 microns. Airflow impedance is 
critical because vital aircraft systems must not starve for 
ample air flow. ARINC 600 dictates an air weight flow of 220 
kg/Hr/KW at 40°C and 1013.25 mbar. No exterior side sur¬ 
face is to exceed an average of 60°C and no single point may 


* MCU is defined as follows: “The Modular Concept Unit (MCU) is the basic 
building block module for use in commercial airplane avionics system de¬ 
sign. This specification provides for the standard interfaces between the 
MCU-type Line Replaceable Units (LRUs) and the electrical wiring, envi¬ 
ronmental control systems, and supporting structures.” 
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exceed 65°C at an inlet temperature of 55°C (ground) and an 
altitude of 15,000 feet 

Electrical and Thermal Parameters 

The following parameters in the electrical system have been 
specified: 

• The power supply in this unit must have 130 watts of out¬ 
put with a minimum power factor of 0.90 and a minimum 
efficiency of 80%. 

• The nominal input voltage of 115 volts AC at 1.6 amps and 
will have an input of 184 VA. This gives a system thermal 
dissipation of between 165.6 and 184 watts. 

ANALYSIS 

The scope of this chapter does not allow an electrical and/or 
electronic analysis of the CTU. The mechanical analysis focuses on 
the following three areas per ARINC 600 requirements: 

• Thermal analysis to ensure that average surface tempera¬ 
tures do not exceed 60°C and that no localized temperature 
exceeds 65° C. 

• A static analysis will be performed to investigate a 12 g’s per 
axis load per FAR Part 25 regulations and RTCA/DO-160C 
procedures. This analysis will ensure that the ultimate 
stresses are not exceeded and is designed to ensure that the 
unit will remain in tack. Functionality of the unit after test¬ 
ing is not an issue. 

• The last analysis is to test for random vibration per 
RTCA/DO-160C requirements. 

Reliability and Mean-Time-Between-Failures (MTBF) calcula¬ 
tions are yet an other aspect of analysis which will be dealt with 
briefly. Depending on the field application of the unit, either MIL- 
HDBK-217F (22) or Bellcore Standard (23) TR-332, Issue 5 may be 
used. In this application, Bellcore Standard TR-332, Issue 5 will 
be used to derive an estimate of the lifetime of the electronics. 
Other methods may be employed to calculated the mechanical (fa¬ 
tigue) lifetime of the packaging itself. 
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Thermal Analysis 

The CTU has the following components and maximum poten¬ 
tial heat dissipation distributions: 


1. PCB 1- uP Board, 26 watts 

2. PCB 2- Memory board, 26 watts 

3. PCB 3 - Network board, 19.5 watts 

4. PCB 4 - Network board, 19.5 watts 

5. PCB 5 - Power Supply, 54 watts 

6. PCB 6 - El Network Board, 19.5 watts, quantity 2 


This totals 184 watts. ARINC 600 requires that no more than 220 
kg/kW/hr flow be used. The task, therefore, is to find out whether 
the system is within ARINC’s airflow budget. 


_ Kg AnA W \kW An An Kg 
220—xl84 — x-= 40.48- s 


kWhr 


1 1000IT 


hr 


Kg 

In other words, our system may not use more than 40.48—of 

hr 

available cooling air. To determine whether this is a sufficient 
amount, recall equation 5.9: 


m = ——— where m is mass flow rate 
Cp^T 

This equation may be rewritten by substituting the value of 
C p =0.2789 for air and maintaining units of Watts for Q and 

Celsius for temperature: 

Rewrite this equation and solve for AT : 

A7 ,_ 3.58560 
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40.48 


In other words, with the available flow our temperature rise 
will be 16.3°C. ARINC requires a temperature rise not more than 
10°C. This flags a problem. 

The first reaction maybe to anticipate changing the electronic 
design. Since this analysis is taking place after the design cycle is 
completed, any changes in the electronics will be extremely expen¬ 
sive. An alternative is to study the standards further and investi¬ 
gate any other possible solutions. 

ARINC procedures were developed for standard avionics rack 
systems- implemented in the electronics bay of the airplane. The 
rack system provides coolant air to all electronic equipment in¬ 
cluding the vital systems intended for survival of the flight. How¬ 
ever, it is possible to develop a private rack which delivers coolant 
air independent of the standard system. The source of coolant air 
for this “private rack” may be the cabin air circulating through the 
system. This is a way to eliminate any regulatory or airframe 
manufacturer restrictions on the availability of coolant air flow. 

Based on equation 11.1, the auxiliary source must supply a 
65.88 kg/hr flow rate for a 10°C inlet to outlet temperature rise. 
The step to calculate the required flow rate as well hot spot tem¬ 
perature were outlined in Chapter 5. There, it was determined that 
the maximum component temperature will be 69°C with a flow rate 
of 31.7 CFM. To provide this flow, a 400 Hz fan will be utilized 
which will deliver a minimum of 36 CFM. 


Load Carrying and Vibration Analysis 

To study the characteristics of the system under various me¬ 
chanical loads, a finite elemen t model utilizi ng 4446 shell elements 
was developed as depicted in Figure ll.lj A static analysis was 
conducted where 12g loads where applied along the three main 
directions. It was assumed that the front and back panels did not 
contribute to the load carrying capa bility of t he sys tem. Von Mises 
stresses as shown in |Figures 11.2| through |l 1.4 were below any 
yield values. The system will be able to withstand the required 12g 
static loading. 

Furthermore, the natural frequencies and mode shapes of the 
system were calculated, the system was subjected to a random vi¬ 
bration excitations in three axes per RTCA/DO-160C proce dures. 
The input level was 0.01 g 2 /hz for 10 to 2000 hertz range. 


Table 
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11.1 


gives th e first five natural frequencies of the system and Fig- 


lure 11.5~l and 11.6| present the first mode shape as well as the over¬ 
all 3S random stresses generated in the CTU. The stresses gener¬ 
ated by random vibrations may be used to estimate a mechanical 
life time for the PCB’s and the box. 



Figure 11.1. The Finite Element Mesh Depicting 4446 Shell Elements 



Figure 11.2. Stresses Under 12g loading in X direction 
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Figure 11.3. Stresses Under 12 g loading in Y direction 



Figure 11.4. Stresses Under 12 g loading in Z direction 
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Figure 11.5. The First Natural Mode of Vibration (Exaggerated) 



Output Sot R«pon» data fan DYNRE3 
Contour: Plate TopVonM*« Sbe» 

Figure 11.6. The 3S Stress Levels Produced By Random Vibrations 


Table 1. The First Five Natural Frequencies _ 

Mode 1 2 3 4 5 

Frequency 185.0 214.2 218.3 220.5 223.6 

The Lifetime may be estimated using the following relation- 
ship(14): 
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T = 


B 


( 11 . 6 ) 


sec . 


4 fa b 0.271* 



b 


+ 0.0428 


Where /is the response frequency. B and b depend on the material 
used: 

For G10 Glass Epoxy: b = 11.36, B = 7.253 x 10 55 psi/cycle 

For Aluminum: b = 9.10, B = 2.307 x 10 47 psi/cycle 

The Lifetime of this unit is essentially infinite. Per RTCA/DO- 
160C standards, any laboratory testing does not exceed more than 
2 hours. This length of time is only sufficient to screen for systems 
to fail in their infancy and unfortunately may not be indicative of 
field behavior and reliability. 

Reliability and MTBF Calculations 

The reliability calculations for this system were based upon 
Bellcore TR-332(23), Issue 5, rather than MIL-HDBK-217F(22). The 
reason for this is that the system is defined as telephone equip¬ 
ment rather than avionics. TR-332 also allows for greater use of 
commercial grade parts and assigns high quality factors to military 
standard components. The following equation was used to generate 
the system failure rate: 


(11.7) 


A - A G KQ7t S 7t T n E 


where: 

A = Failure rate (FIT’s) or failure per 10 9 hours 

A g = Parts generic failure rate or failures per 10 9 hours 

Uq = Quality factor 

jt s = Electrical stress factor 

; rt T = Temperature stress factor 

n E - Environmental factor (a factor of 6 was used) 

First Year Electronic Failures 

It is not practical to provide a bill of material (BO M) for this 
CTU here. However, for a given BOM, a table similar to Table 9.1 


Copyright © 2003 by Marcel Dekker, Inc. All Rights Reserved. 






may be constructed. The failure rate for this system has been cal¬ 
culated as: 

A, = 692 FIT’s and may be rewritten as A, = 69 x 10‘ 9 hours 

The AT&T Reliability Manual (21) defines the system survival rate 
of the first 12 months (8760 hours) per the following equation: 

S t = e~ M = 0.994 

In other words, 99.4% will be in service after 1 year without re¬ 
pairs. 
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